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Abstract 

A global tinie-discretized scheme for the Navier-Stokes equation 
system in its Leray projection form is defined. It is shown that the 
scheme converges to a bounded global classical solution for data which 
have polynomial decay at infinity. Furthermore, the algorithm pro- 
posed is extended to the situation of initial-boundary value problems. 
Algorithms constructed in a different context (cf. [H [TUl [3112]) may 
be used within the proposed scheme in order to compute the solution 
of Lcray's form of the Navier-Stokes system. The method can be ex- 
tended to incompressible Navier-Stokes equation systems on manifolds, 
which is shown in a subsequent paper. 

2000 Mathematics Subject Classification. 35K40, 35Q30. 

1 Introduction 

In its Leray projection form the incompressible Navier-Stokes equation is 
a semilinear partial integro-differential equation system, where the integral 
term is quadratic with respect to the gradient of velocity. This integral term 
requires a careful treatment in order to control the growth of the solution 
(whatever scheme you propose). This is the main difference to the multi- 
variate Burgers equation. Indeed, for the multivariate Burgers equation, i.e., 
the Cauchy problem 



(1) 

. u(0,.) =h, 

on [0,cx3) X M" (where u is some strictly positive constant, i.e., > 0), and 
1 < i < n) global solutions may be constructed via the a priori estimate 

maxsup |Mj(t, x)| < maxsup |/ij(x)|, (2) 
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which may be obtained from estimates of the form 

d_ 

dt 



^\\u{t,.)\\H^ <\\u{t,.)\\Hs+.Y. \\D''u,D^Uj\\L2-2\\Vu\\ls (3) 



«J l"l + l/3|<s 

for some positive s G M (this is the standard notation for Sobolev spaces 
which will be defined below). This type of estimate is valid also for 
initial- value problems on [0, oo) x Q with periodic boundary conditions which 
are defined on the n-torus T". From this point of view the incompressible 
Navier-Stokes equation system is a multivariate Burgers equation system 
(of Cauchy type) with an external force, and with the negative gradient 
pressure as source terms. In addition an incompressibility condition for the 
velocity is satisfied, i.e., the divergence of the velocity is zero. Indeed, the 
Navier-Stokes equation system is defined via three equations, the nonlinear 
diffusion equation 

— -zyAv + (v V)v = -Vp + fe^ t>0, xgW, (4) 
the incompressibility condition 

V • V = 0, t > 0, X E M", (5) 
and the initial condition 

v(0,x) = h(x), xGR". (6) 

The Leray projection for this equation system is obtained from a Poisson 
equation for the pressure, i.e., we have 

-Ap = Elk=i (ail-.) (4^0 ' 

with the solution 

P(,., = -/^^K,.(.-.)E(||£)('.!')*- («) 



Here, 



^ In if n = 2, 
Kn{x) := { " (9) 
1^ \ Ixp-", if n > 3 

is the Poisson kernel. Here, |.| denotes the Euclidean norm and a;„ denotes 
the area of the unit n-sphere. Since the function x — )■ \x\^ is locally integrable 
if and only if /i > —n we observe that for n > 3 the functions 
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are integrable around x = for 1 < I < n (you may observe this explicitly by 
writing the derivative of the kernel in (jlOp in polar coordinates). Note 
that for all x £ M"' we have a natural upper bound 



d 

dxi 



but the derivatives of the velocity in the expression for the pressure in ([5]) 
should have some decay at spatial infinity in order that the Leray projec- 
tion form makes sense. Later, we shall observe that a solution v can be 
constructed which decays at spatial infinity such that the integral on the 
right side of ([8]) exists globally. More precisely, we shall see that for initial 
data h = {hi, • • • , /i„)^ which are in Sobolev spaces [H^ (M")]" for s large 
enough (i.e., for hi £ H" (M") for 1 < i < n) we have a certain decay at 
spatial infinity of the solution such that 



9f dvj 
. ^dxj dxk 



E 



(«■!/) 



dy < CO. 



(12) 



Note that a global classical solution justifies that we speak of 'the' solu- 
tion. Hence, it makes sense to define the Navier-Stokes Cauchy problem for 
divergence free velocity fields by 



at ^ z^r- 



=1 dxj 



j,k=i 



dxj dxf^ 



{t,y)dy, 



(13) 



I v(0,.) 



for 1 < i < n, and where external forces fgx are set to zero for simplicity. In 
order to have polynomial decay at spatial infinity of the data h we take a 
standard assumption that we have h € [ff'^]" for all s € M. This assumption 
may be weakened a bit (this may depend on the regularity which you want 
to achieve- cf. below), but from an algorithmic or physical perspective it is 
satisfying. Note that here and in the following we write 



dvk dvj 
dxj dxj. 



, , Sufc , , dvj , , 



dxi 



dxi 



(14) 



for the sake of brevity (similar for sums of functions etc.). We look for 
classical solutions in the space of divergence free vector fields, i.e., in the 
space 



V G 



C^^'^([0,oo) X M") 



div V = 



}■ 



(15) 



where the solution v should be bounded and regular (i.e. derivatives should 
exists in a classical sense such that uniqueness is guaranteed), and where 
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some decay at spatial infinity guarantees that ()12p is satisfied. Here, the 
function space CjJ'^ ([0, oo) x M") is the space of scalar functions which have 
bounded time derivatives up to first order and bounded spatial deriva- 



tives up to second order, and the function space 



C^'([0,oo) 



is 

the space of vector- valued functions v := (ui,-- - ,Vn)'^ with components 
Vi E C(|'^([0,oo) X M") for all 1 < i < n. Here, and in the following we 
write vector- valued functions in boldface letters. From this point of view 
it seems natural to measure the construction of the Navier-Stokes solution 
V = {vi,--- ,Vn) in a |.|i^2-iiorm of functions with globally bounded time 
derivatives up to first order and spatial derivatives up to second order (for 
each component function Vi ) plus an integral norm related to (jl2p . Note that 
the norm |.|i^2 does not lead to a Banach space, i.e., the limit of successive 
approximations with finite |.|i^2 norms will not have a finite |.|i^2-iiorm. 

The incompressible Navier-Stokes equation cannot be solved by a simple 
global fixed point iteration. Similar as in the case of the multidimensional 
Burgers equation (cf . ) we choose a time-discretized scheme and construct 
fixed points which are local in time. However, compared to the multidimen- 
sional Burgers equation controlling the growth of the integral term in (|13|) 
is an additional difficulty. For Z > 1 let 

v^''= ,vP'^Y : [1-1,1] xM"^M" (16) 

be a solution of ()13p on the domain [Z — 1,Z] x M" in transformed time 
coordinates r with 

t = PIT (17) 

and with initial data v''''(Z — 1,.) = v''''~^(/ — 1,.) being the final data 
of the previous time step number I — 1, where v'''^(0, .) = h(, .). Note that 

vP''(/-i,.) = v(e!; =1 Pmj ■ j for / > 1. Here, the idea is that choices pi < 1 
small enough lead to a local converging iteration scheme in transformed 
coordinates r on domains [/ — 1, Z] x while on the other side the numbers 
pi may be chosen large enough such that the sum of time-step size pi diverges. 
The latter condition implies that the scheme is global in time (the discussion 
here is preliminary; cf. below for a more extensive discussion of this time 
discretization). It seems difficult to control the growth of the solutions vf'' 
directly if the time step size series {pi)i is large enough such that 

N 



^ Pit CO as A^too- (18) 



1=1 



However, the requirement (jlSp for a global scheme is weaker than the re- 
quirement of a scheme where the time step size is bounded from below in- 
dependently of the time step number I. Indeed the requirement (fT8]) allows 
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us to choose time step sizes pi of order 

Pi-\ (19) 

in order to define a global scheme. This does not mean that we shall con- 
struct a solution which has a linear bound with respect to time. Indeed 
we shall construct a bounded solution v where all components Vi are gob- 
ally bounded with respect to the |.|i^2 norm. However, choosing time step 
sizes of order (fT9|) is useful in order to control the integral magnitude (fT2]) . 
However, linear growth of the solution with respect to a supremum norm up 
to a certain order of derivatives is not sufficient in order to define a global 
scheme. 

Therefore we shall construct a bounded solution of an equivalent prob- 
lem with solution v'"''' defined recursively via a series w'^'P'K Here 'bounded' 
means bounded with respect to the supremum norm. For the integral mag- 
nitude [12] it suffices to have linear growth in time. The choice of the time 
step size [19] will ensure the local convergence of the scheme. Here for each 
I > 1 the functions v^''''', v'''', r' are all defined on the domain — 1, Z] x M". 
For T = I — 1 we have v^''''' = v^'^''~^ for all I > 1. The series (r')^ is de- 
signed in order to control the source terms of the equations for v'''''''. For 
integers / > 1 (time-steps) we shall construct a series of real numbers (pi) 
and a family of recursively defined functions (r-) , Z € N, 1 < i < n, where 

rl:[l- 1,1] X M" ^ M, (20) 

and consider the equation system for the functions 

vl'''' = vf'^+rl (21) 

It is not sufficient to construct a global upper bound for -|- r of course. 
We shall construct a global upper bound for r = (r');>i = {r[, ■ ■ ■ ,r\^)i>i 
and a global upper bound for v*"'^ and shall conclude that v = v'' = v''''' — r 
has itself a global upper bound. For / > 1 the functions r\, 1 < i < n (which 
determine a global function r) satisfy themselves for each time step number 
I a multivariate partial differential equation system which is determined 
dynamically within the solution scheme for v^'^ and r. In order to explain 
this in more detail we first consider the equations for the functional series 
v''''''', I > 1 for a rather arbitrary class of functions r- G C^'^ ([Z — 1, Z] x M"), 
i.e., for all 1 < z < n the functions {t,x) — )• r'(r, rr) are assumed to be 
bounded functions with a bounded time derivative ( in a weak sense at the 
integer values I > 1 and in a classical sense elsewhere) and bounded spatial 
derivatives up to second order on the domain (1 — 1,1] x M."". Furthermore, 
the functions are defined recursively with respect to the time step number 
I. For r = (ri, • • • , r„)^ : [0, oo) x M" M" we have 

ri{T, x) = rl{T, x) iff (r, x) e [1-1,1] x M". (22) 
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We construct global functions v'"''' : [0, oo) x M" — > R" and r : [0, oo) x M" 
time step by time step. First we prove that for some bounded function r the 
function {t,x) — >• v''(t, x) = v^'''(r, x) is bounded and a solution of a equa- 
tion system which is equivalent to the Navier-Stokes equation system. Then 
we conclude that there is a solution of the Navier-Stokes equation system 
V = V -|- r which is bounded and in C^'^ ([0,oo) x M") (the space of func- 
tions which have classical time derivatives of first order and classical spatial 
derivatives of second order). Bounded classical solutions with polynomial 
decay at infinity lead to further regularity and then to uniqueness. Indeed, 
it is well-known that energy estimates imply that the Navier Stokes equation 
'determines' v among the smooth solutions. Next for each time step number 
I the restriction v^''''^ of the function v^''' to the domain [/ — 1,Z] X R", i.e., 
the function 



^iPi' Pi' I I 



(23) 



satisfies 



, ^r'^'P'\l - 1,.) = v'^'''''-i(/ - 1,.), 



(24) 



where 



I v^n / dv'-''''^ sr^n r,p,l dr\ 



-pi 



Note that the functional L^'\r^; v^'P'^) is affine with respect to s v''''''' and its 
first derivatives. Sometimes it is useful to consider the linear part Lf''''^(r'; v'"''''') 
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of the functional L^' (r'; v^'^'^). Therefore we write 
= :Sf' (r')+Lf'°(r';v-W). 

(26) 

Note that the operators L^' and L^' ' are determined by the relation ()23p 
and the Navier-Stokes equation system. Furthermore, in (j24p we have 
v^'/''0(0,.) = h in case Z = 1. Note that Lf''°(r'; v'"''''') is a linear differ- 
ential operator with respect to the function v^'''''. Note that all terms of 
the functional L^'' have a factor pi. The numbers pi measure the time step 
size in original coordinates and will be chosen such that at each time step 
we can construct a local time solution of the incompressible Navier-Stokes 
equation by an iteration procedure where we use polynomial decay at spa- 
tial infinity of the iterative approximations of the solution. Furthermore, 
the functions r' are chosen such that at each time step / that we can control 
the source terms on the right side of (|24|) depending on the computation of 
the data from the last time step. This is a crucial idea of our method. More 
specifically 

a) the choice of pi ensures the time-local convergence of an iteration 
scheme which computes approximations v'''^''^'' for the local solution 
^'r,p,i Qf ii^Q incompressible Navier-Stokes equation at each time step 
I, i.e., the solution on the domain — x M" with respect to trans- 
formed time coordinates r. For each time step number I > 1 the local 
iteration determines the function v''''''' as a limit of the functional se- 
ries v'"'^''^'', fc > 0. For given r and each A; > the function v'''^''^'' 
satisfies 

n,7yV" ^ ""i I n, V" r,p,k-l,l dv^'"' 

—g^ Pl^ 2^j=l—Q^, <-Pl2^j=l^j dxj 

= Lf''(r';v^'/''*^-i'0 + 

v'-.P.fc.«(/_ 1,.) = v'''^''^i(Z- 1,.). 

(27) 

Here, in case of fc = (resp. A; — 1 = — 1) we define 

^r,p,-i,Z(^^ rr) = ^^'■'P'^-\l - 1, x) for all (r, x) e [I - 1, /] x M". (28) 
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We shall determine a the time step size 

C 

Pl = T 



(29) 



for some constant C > where all local schemes for v'"'^''^'^ converge 
while Ylff^i t CO as M t oo. The latter property of the time step 
sizes makes the scheme global. There is a difference to the previous 
scheme for the multivariate Burgers equation which has a time step 
size of order l/l (cf. [^). In the present scheme we use the time 
step size ()29p in order to keep the additional integral term of the 
Navier Stokes equation ('additional' compared to the multidimensional 
Burgers equation) under control. Note that it is only for analytical 
reasons that we use a decreasing time step size of this form in order 
to make the scheme global (because the sum of time step sizes goes to 
infinity as the number of time steps goes to infinity). A scheme with 
increasing time step size is preferable from a numerical/computational 
point of view. Once we have proved that the solution function v is 
bounded itself we can alter the scheme for numerical purposes to a 
scheme with increasing time step size. 

/3) In order to determine the function for each / recursively we assume 
that — 1,.) and r'~^(/ — 1,.) have been computed (in case 

I = 1 the function v^'^{l — 1, .) equals the initial data function h which 
is given, and will be set to zero). Given this information of the 
previous time step we first determine a function (pl which is constructed 
in order to control the growth of the functions r-. Furthermore, the 
functions (f>[ and r' are constructed in order to control the growth of 
the local Navier-Stokes solution wp'''' and r- at time step 1. At each 
time step we control the growth of the local solution v''^'' as a limit 
of an iteration procedure of successive approximations 

k 



where we control the growth of the successive approximations of this 
series, and then show that this control is preserved in the limit k oo. 
The advantage of this approach is that all members of the series 
in (j30p are solutions of linear equations. At time step / > 1 we 
compute the function r' in terms of the data of the previous time 
step first. Having determined this function r' appropriately we show 
that the growth of the function v^'^''''' is controlled. Then we shall 
show that for small time step size pi the sum with the summands 
§-^r,p,k,i _ ^^r,p,m,i _ ^r,p,m-i,i^ q£ (j5U|) is small enough such that some 

growth property established for v'"'''''^'' is preserved for all approxima- 
tions v^'P'^''- and in the limit for v'"'''''. At the same time all time step 




(30) 



m=l 
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size numbers pi have a lower bound of form ()29p . Hence, the scheme 
is global. We exploit that we have some freedom in the choice of the 
functions . The construction is such that the global function r which 
equals the local function r' on each domain [Z — 1, x M" is bounded. 
Each r' is itself a solution of a linear equation with with certain source 
terms ip^ which serve as 'consumption terms' in certain regions of the 
domain [/ — 1, /] x M" and control the growth of the functions for v^''''^'' 



and r . For each 1 < i < n the function cj)^ is determined in terms of 
the information gained at the previous time step, i.e., in terms of the 



functions x — t- r' ^{l — l,x) and x 



^{l — l,x). For simplicity of 



notation we denote these functions by rl~^{l — 1, .) and — 1, .) 

sometimes. The functions are constructed grosso modo as follows 
(details will be given in the proof of the main theorem below). For 
each 1 < i < n and at each time step I the function (p[ is constructed as 



a sum 



+ 



in two steps. The function (j)^'^ is constructed in 
terms of v'''^'^~^(/ — 1, .) in order to control the growth of the function 



Similarly, the function is constructed in terms of r and the 
function v'''''''~^(/ — 1, .) in order to control the growth of the function 
v'^''P. Assume that for 1 < i < n we have 



r,p,l—l 



(/-l,.)|o= sup ^■''■'-^(/-1,.)| <C 



(31) 



for some C > 0. We consider two sets where the modulus of these 
data exceeds a certain level. Let 



D 



X\V, 



{l-l,x) £ 



c 



c 



and let 



D 



vA-l 



X\V, 



^' 2 



(32) 



(33) 



Then there are several ways to define functions (j)^''' on domains (/ — 
1, Z] X which will serve as source terms and 'act against the tendency 
of becoming large' which is encoded in the sets ^ ^ and ■ ^ to 
some extent. Maybe the most simple form is to define these function 
(/>^'' in a form with no real time dependence. (This has some slight dis- 
advantages because we have to do some extra work concerning joint 
regularity of the solution functions v'"''''^ at the discrete times I, but 
we can deal with these problems without too much difficulties (cf. dis- 
cussion in the following remark below). Hence, we shall define (details 
are given in the proof of the main theorem) 



: {1-1,1] X 



[T,X 



(34) 
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where 



1,2 



is determined such that 

-1 if X G D 



1 if X G 



vl-l 



and such that 



sup 

x6K" 



!V,S,l 



ix)\ < 1. 



(35) 



(36) 



This is only grosso modo what is going on. The uniform upper bound 
for the functions (p[ with respect to the |.|i^2-iiorm ensures that the 
distance 



dist D 



D 



> c 



(37) 



for some constant c > which is independent of the time step number 
Furthermore, we shah extend the definition of the functions (/)[ to 
the whole domain [/ 



1,/] X M" such that 

0^(r,.)iL2 <c:(a + /c,; 



(38) 



uniformly in the time variable r G [/ — 1, /], where denotes the usual 
Hilbert space of squared integrable functions. Indeed, we shall have 
related = (R")-estimates for the functions and v^'^'', where 
denotes the Hilbert space of functions with derivatives up to 
second orderNote that the estimate in (j38p is linear with respect to 



the time step number I. Therefore our choice of /)/ ~ } is appropriate. 



The estimate (j38p is closely related to the estimation of the term 



_d_ 



Knix 



j,k=i 



dvk dvj 
dxj dxk 



{t, y)dy 



(39) 



in the Leray projection of the Navier stokes equation. Indeed, esti- 
mates of the form ()38p for imply that we have a similar linear upper 
bound for (|39p in case n = 3. In the case of ?i > 3 our method may be 
adapted by use of certain Banach space norms i?**'^. 



Next for the construction of 



assume that for 1 < i < n we have 



sup |r^i(/-l,.)| 



(40) 



where will be a constant which depends only on the initial data 
dimension and viscosity. We shall construct r' such that |r'|o < 
with respect to the supremum norm |.|o, and for all / > 1, and all 
1 < i < n. Let 



D 



{x|rK/-l,x)G [f } 



(41) 
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and let 



D 



r,l-l 



91' 

2 



: [1-1,1] X 



where € C^'^ is determined such that 



(42) 



(43) 



rA-l 



and such that 



sup 



-1 if X G D. ,^ 
1 if X G 

6r''(x)| < 1. 



(44) 



(45) 



Furthermore the extension of (/)[ to the wholedomain is such that the 
norm has a hnear bound with respect to the time step number / 

Remark 1.1. The definition of the source term functions (/)[ is such that 
we may have bounded jumps at the time discretization points I of the 
scheme. However, in the representation of the functions ■yp''''^'' and 
the bounded source terms become integrated over time starting from 
T = I — 1 at each time step /. This will ensure that the construction 
of the global solution function v^'^f and are uniformly bounded con- 
tinuous over time. We shall see that this is sufficient for us. However 
it is also possible to construct global functions (pi which equal (/>^ on 
each domain (1 — 1,1] x M" which are differentiable from the beginning 
defining 



,v,l, 



[T,X 



1,1] X ^ 
:= sin^ (7r(T 



and 



: [1-1,1] X 



(46) 



(47) 



'(r, x) := sin^ {it{t — {I 

Then we get global bounded solution functions v'^'P and r which are 
differentiable across time points t = I. However the growth control 
of the local functions v'^''''^ and becomes a little bit cumbersome. 
Hence we proceed with the definition which leads to global solutions 
which are bounded continuous at time points r = / first. 
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We continue to outline the main ideas within item The functions 
serve as source terms or 'consumption terms' in the following parabolic 
equation for ( which is a linearising of a time-local Navier-Stokes 
type equation with source term (f)\ among other source terms), where 
for 1 < i < n we have 

-PI kr. {-tM^ - y)) ^,m=i [^^'-^) H - 1, v)dy + 

(48) 

and where 

tP-1,0/ l^l rol-l\ — I V^n p,l~l di\'''''~^ , v-^n r,p,l—ldr'~^ 

Li (r' ^ V-'^'' ') = +pi J2j=i r]' -^Q^ + PI Ei=i vf' 

(49) 

Note the signs on the right side of (|48p . Note that dependence on the 
function V'^ involves only the function x — )• v^'''''~^(Z — 1,.) which 
is known at time step /. Assume that this equation has a classical 
solution = {r\, ■ ■ ■ , r^) in C^''^ on the domain [/ — 1, Z] x M" (we shall 
show existence of such classical solutions for time-local Navier-Stokes 
type systems below). Plugging (jlH]) in (f27|) for A; = (note that for 
k = Q y'.'P^^~^'^ = = v^j''''''~^) we get the equation 

{n r,p.O,l o2 ^'iPiOJ , r,p,0,Z 

^ — /'i'^ Ei=i + /'i Ei=i ^/'^ - 1, ■) = + 
^r,p,0,«(; - 1,.) = v-^'P'^-^l - 1,.), 

(50) 

where we shall point out that 5' is small compared to (pl in the relevant 
regions Df^-, D'/^-, Df^+, and (provided that pi is small 

enough) . 

The choice of (pi is such that in these relevant regions the behavior of 
the solutions v'''''''''' and r' is controlled. Moreover, we know that the 
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function 

^r,p,l ^ ^r,p,0,l ^ ^^r,p,k,l ^ (51) 
k>l 

is small where the sum of functions 

§^r,p,k,l ^ ^r,p,k,l _ ^r,p,k-l,l^ ^52) 

is a 'perturbation' of the function v*"'^'*^'' if pi is small We may choose 
pi in such a way that a certain control of the behavior established 
for v''''''^'^ is also valid for all functions v^'^''^'', A: > 1 and for the 
limit v'''^'^ Proceeding with time step numbers I in such a way while 
the sum "^"^i pi of the time step numbers pi goes to infinity ensures 
boundedness of the function v^'f. However, boundedness of the latter 
function does not guarantee boundedness of the difference v — r = 
v''— r. We have to ensure that r is bounded, i.e., there is a global bound 
of all functions r' independently of the time step number I. We shall 
see later in detail that the construction outlined can be implemented 
such that the functions v^'^'' and r' are uniformly bounded, i.e., there 
is a number C > such that < C and |r'| < C for all / > 1 

and 1 < i < n. The constant C > can be computed a priori and 
depends only on the dimension n, the viscosity > 0, and the initial 
data h. Moreover, the sequence defines a globally bounded function 
r : [0, oo) X M" which is (weakly) differentiable with respect to time 
and has bounded spatial derivatives in a classical sense. The functions 
v'^'/''^ : [/ - 1,/] X R" ^ M", I >l which solve the equations ([Ml) for 
each/ with v/'''^''(/ - 1, .) = v/'''^''"i(/ - 1, .) for / > 2 and v^''^''(/- 1, .) = 
h(/ — 1, .) for / = 1 define a global function 

v^'^P ■=vP + r: [0, oo) X ^ M", (53) 

which equals v''''''' on each subdomain [/ — 1,/] x M", / > 1. Since r 
is bounded this function v^'*" satisfies a system of equations equivalent 
to the Navier-Stokes equation system and a bounded global solution 
of the Navier-Stokes equation system is given by 

v(t,.) :=v^'''(T,.)-r(T,.). (54) 

The regularity v G [C^'^ ([0, oo) x M")]" is easily obtained if v'^''' and 
r are globally Holder continuous and have a certain decay at infinity. 
This leads to uniqueness. 

However, before we get involved with the construction of the sequence 
(pi) and (r^) in detail let us consider ()13p again. From the point of view 
that the incompressible Navier-Stokes equation is an extension of the multi- 
variate Burgers equation a solution of the equation is an element of the 
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class of divergence-free solutions of a certain family of multivariate Burgers 
equations where for each function f we add a certain source term to ([1]) 



involving functions f G 



. More precisely, consider the 



family (v'^)^ where for each fixed f the function v-^ satisfies the Cauchy 
problem 

at ^ 2^j=i 2-.j=i - 

{£-Knix - y)) E-.=i (S^it) (''y)dy, ^''^ 

^rf{0,.) = h. 

Then a fixed point of the map (defined on regular divergence free vector 
fields with values in divergence free vector fields depending on time) 

Floe : f — v-^, where divf = 0, and divv-^ = (56) 

is a solution of the incompressible Navier-Stokes system. Let us denote a 
fixed point of the latter map by 

v*=v^*. (57) 

Formally, such a fixed point has a representation in terms of the fundamental 
solution r* of the scalar equation 

where we denote v* = {v\, - ■ ■ , I'j^)^, i.e., we have the formal representation 
{l<i<n) 

v*{t,x) = hi{y)T*{t,x,0,y)dy+ 

/o /r" (£-^n{y-z)) Elk=i (Stfe) {s,z)r*{t,x,s,y)dzdyds. 



(59) 



Note that for all 1 < i < n we have the same fundamental solution T* 
of a scalar equation involving first order coefficients Vj,l < j < n which 
are the same for all 1 < i < n. Well, if we construct such fixed points in 
a time-discretized scheme, then we may use the a priori estimates of the 
multivariate Burgers equation of type and a representation of type 
at each time step in order analyze the growth of the solution. Note that we 
cannot control the integral terms involving the Poisson kernel in the Leray 
projection form this way. Nevertheless the representation ()59p gives us a first 
hint how the integral term can be controlled in a time-discretized scheme 



14 



using our considerations above. There is a mutual dependence of the choice 
of the functions r • and the choice of the numbers pi which ensure the local 
convergence in time. The choice of the latter numbers depends on the local 
scheme. It is not necessary to start each time step with the solution of the 
corresponding multidimensional Burgers equation. Indeed, an alternative 
way of constructing a divergence-free fixed point 

: f ^ (60) 

via a time-discretized scheme may be the following: locally in time the family 
(v-^)j. of vector- valued functions v-^ = (^v(,--- ,vCj may satisfy for each 
f = (/i, • • • , /n) the equation 



f dv( _ sr^n d'^vl sr^n r dv{ 
dt ^ ^j=l ^ ^j=l dxi 



3 



(61) 



I v/(0,.) = h. 

Here, we search for a fixed point solution f* with Fgiob(f*) = v-'^*, where 

divf* = divv-^* = 0. (62) 
Let us consider this map more closely. Start with the equation 

— z^Av-f + (f • V)vf = - V/ t>0, X G M", (63) 

at 

for some scalar function , and where 

fejuG C^^'2([0,oo) X M") "|divu = 0, u(0,.) =h}. (64) 

In the following we denote spatial derivatives in the form fij := and 

fij^k '■= Qx Qxfc ^^'^^ ^ usual for example in the literature on the theory of 
general relativity. Sometimes we feel free to denote time derivatives in the 
form /j^i = For the divergence divv-^ we have 

— div v-^ + uA div v-^ + ^ fj ^ divv^ = - ^ fijvl- - Ap^, (65) 

j ^ i,j=l 

where divv''^(0, .) = divh = 0. Now let F'^ be the fundamental solution of 

8 

_r/-z.Ar^ + ^/^T^ = 0. (66) 
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Then the solution to equation (|65p with zero initial data has the represen- 
tation 

divv^(t,x) = / / ( - V fi^.^-ApA {s,y)Tfit,x,s,y)dyds. (67) 

Jo Jr" y J 

Now let v-'^ be a solution of ()6ip . Then we have the representation 
Vi{t,x) = /jg„ hi{y)Tf{t,x,0,y)dy 

+ /o/m"/r" {£-^n{x-z)) T.lk=i (isjfe) {s,z)Tf{t,x,s,y)dydzds 

(68) 

Hence if we solve (|61|) for v-^ in the form (|68|) , then this is the same is solving 
(1631) in the form 



v{{t, x) = hi{y)Tf{t, X, 0, y)dy 
+ /o/m" (^^') {s,y)Tf{t,x,s,y)dyd6 



(69) 



/(t, x) = - / (i^„(x - y)) V ( 1^^) (t, y)dy, (70) 



along with 



and such that 

divv-^ = (71) 

is ensured. Hence, if we find a fixed point v* of the map f — t- v-^, then this 
function satisfies a) the equation ([UT]) , and b) the equation (jU (with f^x ^ 0) 
along with 

r " / f)v* dv* \ 

p{t,x) = - (Kn{x-y))^ (-^-^]{t,y)dy, (72) 

and divv* = 0, and v*(0, .) = h(.). The latter ansatz has the advantage 
that we can easily preserve the incompressibility condition at each step of 
approximation of the local solution, i.e., the solution at each time step. Note 
that we speak of 'the' local solution here in the sense of our construction, 
i.e., the approximating equations involved in each time step have a unique 
solution. Moreover, we shall prove that the local solution is a classical solu- 
tion, and then it can be shown that the classical solution is indeed unique. 
In any case (whether we start with the multivariate Burgers equation or 
with a linearized Navier-Stokes equation), it seems that such a fixed point 
cannot be obtained by a global iteration scheme. 
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Recall that it is indeed sufficient to solve ()13p (this is well-known). A 
regular solution v of equation (jlSp satisfies 



v(0,.)=h, 

along with p of form ()72|) . The associated equation for the divergence div v 
has zero initial data and a source term 



E 



dvj dvi 



■ ■ -I L 7 



dxi dx 

which becomes zero, since 



A.P, (74) 



for regular v with the indicated decay at infinity. Here, A^; denotes the 
Laplacian (where derivatives are with respect to the variables Xj), and where 
some regularity of the function v is sufficient. Note that higher order dif- 
ferentiability is easily obtained if a solution in C^'^ has been obtained (cf. 
Section 4). 

Back to the time discretization, the time step sizes are given by a series of 
positive real numbers {pi)i£N where N denotes the natural numbers (starting 
with / = 1) and the size will be chosen in such a way that the iteration at 
each time step / converges and such that we have global convergence, i.e., 
the time step sizes are large enough. In order to ensure global convergence 
it is sufficient to have a time step size of order pi = j- for a constant C > 
which depends only on the viscosity z^, the dimension n and the initial data h. 
Prima facie it seems that our recursive construction of a bounded function 
r which equals the functions r ■ locally on — 1 , /] x M" leads us to the 
conclusion that a uniform lower bound for the time step size pi may be 
possible (note that the numbers pi are time step sizes from the point of view 
of the original time coordinates). However, we need to have a bound for the 
integral magnitude in ()12p too, and this leads us to our choice of a decreasing 
time step size. From a numerical point of view we hope that the opposite may 
be possible, i.e., that smoothing effects may allow to increase the time step 
size pi as time goes by, i.e., as the time step number / increases. Once we have 
proved that there is a bounded classical solution we may discard the control 
function r and consider this possibility. We shall ensure that the global 
solution (which equals v^'P''- on each domain [/ — 1, /] x M") is bounded. 
Since r is bounded, this implies that we have a bounded global solution 

= v^'^ + r (here, v'' equals v'''' on each domain [I — 1, x M"). Note that 
v(t, .) = v''(r, .) where the former function is a solution of the Navier Stokes 
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equation. Note that we construct the function r time step by time step 
together with the function v'^'''. Furthermore, note that = v^''', where 
the first function refers to the equivalent system in original time coordinate 
t. Accordingly v** satisfies the equation 

PI {£-Kn{x - y)) E • ,=1 (US) it, y)dy + n,u (^6) 

. v'-(0,.)=h(0,.), 
where 

(77) 

{l-Knix-y)) E-,=i (Stfe) (t,2/)d2/ + LO(r;vO. 

Our solution scheme computes v'' and r (in original coordinates) simultane- 
ously, and a global solution of the 

1. Navier 



-Stokes equation system is then obtained by addition v = v'^ ' + r (we use 
the same symbol r for simplicity). Note also that r(t, .) = r(T, .), hence v is 
sounded since is bounded and r is bounded. Note that on the domain 
Ylm=i P'm.i X]!n=i /'"i X 1^" the solution of the Navier-Stokes equation v 
equals v'''' on the domain [/ — 1, x M". As already remarked dependence of 
the numbers pi on the time step number I is for numerical purposes because 
the scheme with flexible time step may take advantage of the smoothing 
effect of the scalar densities involved in our scheme (implying larger time 
step sizes as time goes by (cf. remarks in section 5)). 

Remark 1.2. Note that in [B] we used a scheme with time step sizes 

PI = (78) 
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in order to solve the multidimensional Burgers equation, where C* > is a 
constant which does not depend on the time step number / (it depends only 
on the data h, and the dimension n , and the viscosity constant v > and 
will be determined below). The scheme considered there is not sufficient in 
order to prove global existence of solutions for the incompressible Navier- 
Stokes equation since we have an additional source term which is quadratic 
with respect to the gradient of the velocity. The introduction of the functions 
r • for I > 1 and 1 < i < n is the crucial difference which makes it possible to 
control the source terms in the Leray projection form of the incompressible 
Navier-Stokes equation. 

Let us have a closer look at the time transformations t — )■ r, and the 
introduction of the function r. We consider the Cauchy problem (|13p on the 
domain D = [0, oo) x M"" where our interest is in the case n > 3. We shall 
solve the Cauchy problem in subsequent time steps / > 1 on the domains 

A = m-i,Ti] xM«, (79) 

where Tq = and Ti = Ti^i + pi for / > 1. Instead of considering a 
scheme with small time step size in original coordinates we may consider a 
equidistant scheme in transformed coordinates with time step size 1, i.e., 

r^ti{T)=piT if r G [1-1,1]. (80) 

The transformed domains are denoted by DJ = [1 — 1,1] x R". The original 
time coordinate t has an index here in (jSOp in order to indicate that we are 
actually considering infinitely many different time transformations (one for 
each time step number /). However, in order to keep notation simple we 
drop such indices in the following. The index / makes clear which domain 
we consider. On each domain DJ we have a transformed Cauchy problem 

-w- - z2j=i ^ + Pi Ej=i as- = 

pi /r" {£-Kn{x - y)) Zlk=i (^^y)dy, 

, v^''(/-l,.) =^rP'^-\l-l,.). 

(81) 

Hence the initial data of the Cauchy problem for v^'' on the domain DJ 
(the Ith. time step) are the final data of the Cauchy problem for v^''~^ on 
the domain DJ_^. Note that v^'^(0, .) = h. However, as indicated above we 
introduce another sequence of real functions (r') and consider for each time 
step / > 1 and on the domain DJ the function 

(r, x) v'"'''''(r, x) := v''''(r, x) + r^. (82) 
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If for time step number I the function v'''' solves ()8ip . then the function 
■^r,p,l gQ^ves the equation 

dv^'''''' d'^vl''''^ sr^n r,p,l dv^'''-'' r P,h I r o l\ , 



dT 

J 



y'''P'\1-1,.) =v^'''''-1(/-1,.). 

(83) 

The sequence (r') will be chosen (recursively with respect to I) in such a way 
that the global growth is controlled, i.e., that the components of v''''''' are 
bounded independent of I. Furthermore, since r is bounded this implies a 
bound on v''''. The sequence (pi) is chosen in such a way that for each time 
step number / the local scheme for v'"'^'' converges. At each time step I we 
approximate v''''''' iteratively by functions v'''''''^'' for A; > 0. Indeed, at each 
time step I we shall construct a solution v^''''' of (|83p in form of a functional 
series 

k>l 

where for each /c > 1 we shall have a contraction of the successive approxi- 
mations 

k 

^r,p,k,l = ^r,p,0,l ^ ^^r,p,m,l^ (gg) 

m=l 

where 

^§^r,p,k,l^ _ ^^r,p,k,l _ ^r-,p,fc-l,i| ^ ^^^r,p,k-l,l _ ^r,p,k-2,l^^ j^gg-j 

Here the symbol |.| represents some appropriate norm and we have a con- 
traction constant c € (0, 1) which turns out to be independent of the time 
step number I. In general we shall prove a contraction property of type ()86p 
for specific series {^P'^'^)k defined as in (iii) below with a specific first element 
v^'O;' for each time step I > 1 (alternatively, we could use the methods in (i) 
or (ii)) below to set up a local scheme). We have (at least) three possibilities 
for the local iteration (given (r') which is computed first at each step). 

(i) We start for each I the iteration with the corresponding multivariate 
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Burgers equation, i.e. for each / the function v^'^''''' solves the equation 



dv. 



r,p,0,l 



E^2 
7=1 + Pi z:.,=i 



r,p,0,l 



dxi 



Lf (r'(Z-l,.),v^'A'-i(/-l,.))+r^,,(/-l,.)+ 



PiIm 



a 

dxi 



Knix-y))E 



dxi 



dxk 



{I - l,y)dy, 



^ ^^^Pfi^\l - 1,.) = w'-'P'^-^{l - 1,.). 



(87) 

We shall see below in which sense we have time derivatives at the 
points where r = / — 1 for some / > 1. For the latter equation we have 
the a priori estimates and hence a global regular solution which 
is bounded by the respective initial data. The initial data v'''^''~^(Z — 
1, .) are the final data of the previous time step which are the result 
of an iteration which we have to define next. We define a series of 
multivariate Burgers equations for k> 1 where w^'P'^'^ solves 



E^9 r.p.k.l 



n r,p,k,l dv. 



En r 
i=l 



r,p,k,l 



dx-j 



dxi. 



(t, y)dy, 



.r,p,k,l 



{l-l,.)=yr''P''-\l-l,.). 



Then we may prove - with an appropriate choice of pi , r/-that locally 
this scheme leads to a fixed point. 

(ii) We may linearize for each time step the equation ()83|) and define 
a sequence of function v^'P'^^^ which are solutions of integro-partial- 
differential equations, i.e., we start with the solution v'''''''^'' of 



dv 



r,p,0,l 



d^V 



r,p,0,l 



W Pi" Z-j=l Sxj + Pi ^j=l 

= Lf (r',v-''''0'') 

+Ar + Pi /r" imi^nix - y)) Elk=l 



'■'j dxi 



dxi 



dxi. 



(/ - 1,.) = v'-'P'l-^il - 1,.). 



(r, y)dy, 
(89) 
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Then for A; > 1 we define recursively functions v'''*^'' to be solutions of 



r,p,k,l 



, = 1 -^ZT- + Pi E, = l 



r,p,k-l,l dv. 



r,p,k,l 



dxi 



^ ^r,p,k,i(^l _!,.) = ^■>P^i-^(l - 1,.). 



' (t, y)dy, 



(90) 

Note that this is an iteration of global equations, i.e., linear partial 
integro-differential equations. 

(iii) Alternatively, we may define a sequence of functions v'^'''''^'' which solve 
linear parabolic partial differential equations. We may choose an iter- 
ation starting with v'^'''''^'' which solves 

/ i~i r.p.O.l r^o r, 0,0,1 , ^ f) r, 0,0,1 



j=l —^5— + Pi Lj=l ' 



+Pi /r" (£-^n{x - y)) Elk=i 



(91) 



(t, y)dy 



Then for A: > 1 we define recursively functions w'^'P'^'^ to be solutions 
of 



( d^-'P'^^l a^^'"-'^'' , ^„ r,p,k~l,ldvl''''''-^'' 

Pl^ l^j = l + Pi — 



dxi 



+^lr + PI /r" {£-^nix - y)) Elk=l 

^ ^r'-'P'^^^l - I, .) = v^'^''-i(Z-l,.). 



{t, y)dy, 



(92) 



Let us discuss the different advantages / disadvantages of the three approaches 
for the local solution at each time step / > 1. 

ad(i) The solution of the multivariate Burgers equation is closer to the solu- 
tion of the Navier-Stokes equation than the linear approximations in 
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(ii) and (iii). However, at each substep k of the /th time step the func- 
tion v''''^'^ may not be divergence-free (a divergence free vector field is 
obtained in the hmit). We have a priori estimates (essentiahy a max- 
imum principle) for the multivariate Burgers equation. On the other 
hand from a numerical point of view linearized equation are preferable 
in iteration schemes. 

ad(ii) This iteration has the advantage that we can ensure that at each time 
step I > 1 each approximation we can ensure that v''''^'' is divergence 
free, i.e, divv^''^'' = for all A; > 0. However, at each substep k of a 
time step I we have to solve linear partial integro-differential equations 
which may be complicated from a numerical point of view. 

ad(iii) In this case we have to solve local scalar linear parabolic equations in 
order to determine v''''^'' at each substep fc of a time step From a nu- 
merical point of view this is interesting since we have good approxima- 
tions of solutions for the involved parabolic equations where the second 
order terms form a Laplacian (cf. ([ID]), [1]). On the other hand, the 
local iteration does not take place in a space of divergence-free vector 
fields in general, i.e. a divergence free vector field is obtained in the 
limit. 

In this paper we choose the third alternative approach for the local so- 
lutions and construct for each time step I a fixed point in some appropriate 
function space of maps on domains Df. This allows us to define the most 
efficient algorithm among the three alternatives in section 4. Compared to 
standard discretization schemes all three versions of our analytical scheme 
have the advantage that there is no spatial discretization. Next let us con- 
sider the local iteration in the form of the third item (iii) above in more 
detail. At each time step / we consider maps of the form 

{ ^yr'-'-f'P'^ = F[{f), (93) 

where f = (/i, • • • , /n)"^, and where v^'f'P'^ = ^vl'-^'^''' , • • • , Vn'^'^'''^ satisfies 
the equation 

w/m" [£-Kn{x-y)) E-.=i (S^St) ir,y)dy + rl^, (94) 

, v^'^'/'''(/ - 1,.) =v^'''''-i(/-l,.). 

Here we assume that we have solved the Navier-Stokes equation in trans- 
formed coordinates for r < / — 1 and the initial data at / — 1 are given by 
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the solution at this time. The domain of the map F-[ is 

Dpr := {f e [d'^ ([(/ - 1), /] X M-)] " |f(Z - 1, .) = V'P'^-Hl - 1, .)} ■ 

(95) 

In case / = 1 we have v'''-^'^'^(0, .) = h, of course. With the appropriate choice 
of pi and r • a local scheme for v''''''' may be defined in terms of a functional 
series (v'"'^'^''')^ with lim^i-oo ■v^'P'^^'- = v'"'''''. We start the iteration determing 
for 1 < i < n the functions (f>[ and then the functions r' (solving a certain 
equation as scetched in /3 above). Then for each A; > we have 



' dv'^'''''''^ v^n d'^v'-''''^'^ , sr^n r,p,k-l,ldv. 



r,p,k,l 



t v''.p.'(Z - 1,.) =v'-'P'^-\l-l,.), 

(96) 

along with v*"'^' = v''''''' Let us look at the difference of two successive 
approximations v''''''^"^'' and v*"'^'^'' where we consider fixed functions f 
and g instead of ■v^'^p^^~^'^ etc. as coefficient functions. Comparing v^''''' 
with v^'^'' leads us to an expression for the difference which we denote by 

This function satisfies the equation 



dv. 



+Pl /m^ - y)) Elk=i (§|St) 



(98) 



-PI /m" {£-Kn{x - y)) Elk=i (§|St) y)^^' 

^ ^^rj,g,p,l^l - 1,.) = 0. 

Here we have 

(99) 

j^P,l,f,0 ||j.i^^r,/,p,/^ _ i^P,l,9,0 (^^l ^^r,g,P,l-^ 
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Note that the source term ^ does not appear on the right side of this 
equation for the difference, furthermore all the terms involving only rl and 
the derivatives of (because they cancel out). Hence the first term v^'^'*^'' 
of the functional series v'''''''''' should contain the essential information con- 
cerning the growth of the solution at time step / and the higher order terms 
of the local iteration. In the first iterative time step of an iteration scheme 
related to ([98]) we shall have fj = u^''''^'' and gj = y'^^'P'^'^^ such that the first 
term on the right side related to pSj) becomes 

n rj r,g,p,l n n r,p,0,l 

j=l J j=l J 

In [6] we considered a bound of u^''^'' and vf'j'^ of form C*Ci where C; 
depends linearly on / and C* is a constant independent of the time step 
number I. Linear dependence of Ci of both terms with respect to I is sufficient 
in order to make our scheme global as the sum of time step sizes ^i^i Pi 
is unbounded. However, in this paper we shall construct a uniform bound 
CiC* for v^'j'^''' which is independent of the time step number I. Let us 
stick to this second term of right side of (|98|) for a moment. In the fc-th 
iteration step the contribution to the first term on the right side related to 
becomes 



En 
7 = 1 



r,p,k,l r,p,k-l,l\ dv- 

v/ — v/ 



/ 0. X 

Pl2^j=iy"j M dxj ^ l^ni=l dxj^'^i I' 

such that we may estimate these terms by using a contraction property of 
the differences 5v^'P'^'^ which we shall observe for the specific series 

T^r,p,k,l _ ^r,p,0,l _|_ Y^k ^ ^^r,p,m,l _ 

(102) 

^r,p,0,l _j_ ^^r,p,m,l ^r,p,m— 1,(^ 

Now consider again the equation (|98p . Let us consider r' = for 1 < i < n, 
since we have no source terms in the higher correction terms for Sv^'P'^'^ 
(in the proof of the main theorem below we shall see the additional terms 
j^p,ij,gfi ^j,^ ^^r,p,k,i'^ with f = wP'^''' and g = v''''^'' do not alter the reasoning 
of time local convergence essentially) . Note that for r • = we denote 

^P,k,l = ^0,p,k,l ^ ^r,p,k,l^ (^]^Q3^ 

and similar for all components v^'^'K Since f G Di^j we have that for each 
1 <i <n the fundamental solution of the scalar equation 



i=i i=i 
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exists. Similar, since g G we have that for each 1 < i < n the funda- 
mental solution r' of the scalar equation 



exists. Then formally we may represent the solution of the equation (j94p 
in terms of the this fundamental solution (recall that we consider r = for 
simplicity). Since w^'P'^ and \3,p,i have the same initial data we have for 
1 <i <n (recall that r = this time) 



- lo /iRn Ej=i ifj - 9j) {s, ?/)%-(s, y)r^(T, x; s, y)dyds+ 

Jo P<- /m" /r" ^'^n,i{z — y) X 

Ej,fc=i Uk,jfj,k) {s,y) - YTj,k=i i9k,j9j,k) is,y)jT^j{T,x; s, z)dydzds, 



(106) 

where we denote 

fk, = || etc. (107) 

for simplicity. (Note that in ()106p we may express the first spatial derivatives 
of the function in terms of an integral involving the respective first spatial 
derivatives of the fundamental solution F^). In order to construct a fixed 
point we shall use some decay at spatial infinity. Indeed we need a certain 
decay at spatial infinity of the approximating functions v''''^'' of v-^''''^ in 
order to estimate the limit of the functional series (v'''*^'')^ (with respect to 
some appropriate norm, for example a Sobolev norm I-Ih" for s > ^ + a. 
In order to construct an iteration scheme for the higher order corrections at 
each time step I, i.e., the functions Sw'^'P'^'^ for A; > 1, we observe that 

E",fc=i ifk,jfj,k) {s,y) - Ej,fc=i i9k,j9j,k) is,y) = 
Ei,fc=i ifk,jfj,k) {s,y) - Ei,fc=i ifkd9j,k) {s,y)+ 
E",fc=i i9k,jfj,k) {s, y) - TJlk=i {9k,j9j,k) (s, y) = (;lo8) 
Ei,fc=i {fk,j{s, y) + gk,j{s, y))) x 



E",fc=i {fjAs,y) - 9j,k{s,y))^ . 



In order to deal with the problem of constructing a fixed point in a func- 
tion space of infinite domain we use the standard assumption concerning 
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decay of regular initial data at spatial infinity, i.e. we assume that the 
map X —7- h.{x) = {hi{x),--- ,hn{x)) is a given function with components 
hi in (M") n for all 1 < i < n and s large enough. Here, we write 
= (M"). Now three different procedures are possible in order to con- 
struct a fixed point for each time step / in an appropriate function space, 
i.e. locally in time. The first method a) below makes use of embedding the- 
orems for convergence of local schemes. It is remarkable that Hilbert space 
theory suffices in order to deal with the most interesting case of dimension 
n = 3. The following method can be generalized in order to include arbi- 
trary dimensions by considering embedding theorems for spaces of Holder 
type (so called Zygmund spaces C|). These spaces coincide with classical 
Holder spaces for noninteger values s. We have the standard embedding 
theorem 



Theorem 1.3. For all s G M" and p G (l,oo) 



(109) 



for r 



s — - 



Here the spaces i?*'^ are defined similarly as in the case p = 2 and 
coincide with the spaces H^'P for integers k which have weak derivatives up 
to order k in L^. We do not repeat the definitions here (which can be looked 
up in standard textbooks) since this of marginal importance for us, i.e., in 
order to remark that the following method a) can be realized for arbitrary 
dimension n > 1. We shall see later that the integral magnitude in (|12|) is 
finite. 

a) The following method can applied in the case of arbitrary dimension 
(cf. preceding remarks), but for simplicity we consider the case n = 3. 



This is a construction in 



1, 



where 



Hf := {/ G a {[I -l,l]x K") |/(t, .) G (M") , V t G - 1, /]} . 

(110) 

This is the space of vector-valued functions h = (hi,-- - ,/in)"^ with 
hi G C^'^, i.e., hi is in C^'^ with bounded derivatives of first order 
with respect to time and second order with respect to space. For 
the series v''''^'' defined in {iii) above we can establish that v''''^'' G 



cl'^{[i-i,i]xm-) 



such that for some < c < 1 we have 



r,p,k,l r,p,k—l,l\n 



ll,2> 



where 



i=l 



,j,k\o 

j=l j,k=l 



(111) 

(112) 
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and where |.|o denotes the supremums norm. Then we can show that 
i^[Hfr^.U[Hf]\ (113) 

Then our scheme leads to a series (v''''''^)^ G C]''^ {[I -1,1] x M") n 
\Hl^^ " with a hmit 

Since n = 3 the functions v^'^ are Holder continuous with respect to 
the spatial variable (uniformly in time r) and the fundamental solution 

r^i'' of 

^-^''^E^ + Z'^E^ a^-^f (r'V) = (115) 

exists. Then we can use the representation 

vl''''\r, T, x) = hi{y)r[;^{T, x; 0, y)dy+ 

lo Pi /r" /r" ^ri,i{z — y) X 

E -,^=1 ) y)r:/(T, x; 5, z)dy(izds (116) 

~ !o /m" ^1t(S' y)r[;V> a;; s, z)dydzds, 

+ /o" /m" ''°(*'' v)(s, y)r[;V> 2;; s, z)dydzds, 

where the existence of the first spatial derivatives on right side of 
()116p may be established directly. Existence of second order spatial 
derivatives and first order time derivatives may be derived using partial 
integration. 

Remark 1.4. In a subsequent paper concerning Navier-Stokes equa- 
tions on manifolds we shall show that from a closer look at (jllip we 
may even conclude that for some c S (0, 1) we have a time- local con- 
traction with respect to Holder norms, i.e., 

I r,p,fc+l,/ _ p,k,l\n I ^ \ r,p,k,l _ p,k-l,l\n /i i rj\ 

|V ^ la/2,2+Ql - V lo/2,2+a- K^'-') 

For compact mainfolds this implies that we have a local contraction 
on a Banach space, and this simplifies the analysis. Note, however, 
that the usual the introduction of another time transformation which 
introduces a potential term of a specific sign is critical. A naive use 
of this Schauder type estimates within our scheme does not lead to 
a global scheme. Here we use standard notation for Holder norms, 
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where the first subscript a/2 refers to the modulus of Holder con- 
tinuity with respect to the time variable and the second subscript 
refers to the modulus of Holder continuity of the (second derivatives) 
of the functions v'^'f'^'^ with respect to the spatial variables. More 
precisely, define the Euclidean distance in M""*"^ between the points 
yi = {ti,xi),y2 = {t2,X2) by 



e{zi, Z2) = - t2| + \xi - X2I. (118) 
If It; is a function in a domain D C M""'"^ we denote for a G (0, 1) 

\w{yi) - w{y2)\ 

[W\a/2,o,,D = sup ^ . (119) 

vi^y2;yi,y2(iD e° (2/1,^2) 

Next define 

\wU/2,a;D = klo,D + Na/2,a;D-, (120) 

Similarly, we define 

n 

W\a/2,l+a,D = X] Kj]a/2,a,D , (121) 

and with the notation wt := ^ we have 

n n 

[w]l+a/2,2+a,D = '^[wt\a/2,a,D + [w,j,k\a/2,a,D- (122) 

j=l j,k=l 

This leads to the notation of more regular HoDer spaces 

\w\a/2,l+a:D ■= klo,l;D + [w]a/2,l+a,D (123) 

and 

\w\l+a/2,2+a;D ■= \w\l,2;D + [w]l+a/2,2+a,D , (124) 

where we use the notation 



i;D := \w\o-D + X \wx,\o-D, and (125) 



1=1 



IMOAD ■= \ w\0;D + X \Wxi\0;D + + \Wx,Xj\o;D, (126) 

1=1 *J=1 



and 



\w\l,2;D ■= \w\o-D +'^\Wx,\0;D + \wt\o;D + \WxiXj\o;D- (127) 

i=l ij=l 



29 



Note that the latter norms do not define Banach spaces. However, be- 
low we shall use the fact that a functional series which is uniformly and 
absolutely bounded with uniformly and absolutely bounded derivatives 
can be differentiated term by term. For this matter (|127p is useful. If 
the domain D is determined from the context we shall suppress it in 
notation, especially if D is of the form [S, T] x . 

b) Another possibility (a more elementary way) is the following. We use 
the relation (jllip and apply it iteratively starting with the function 
vP'^-i and get a series (v*"'^'*-''')^ G C^'^ with 

k 

^p,k,l ^ ^r,p,l-l ^ §^r,p,m,l ^ f(jl,2 _ ^] ^ M")j " , (128) 

m=0 

where 

§^r,p,m,l _ ^r,p,m,l _ ^r,p,m-l,l ^ (129) 

Next assume that for all x G M" and for multiindices a with |a| < 2 
and k < 5 we have constants Cak such that 

For each substep k of each timestep / we shall show that for |a| < 2 
the approximations v'"''''*^'' satisfy 

^Qyr,p,k,l^t,x)\ < ^^JJ^^ if l«l < 2 . (131) 

Then the decay relation for the solution helps us to control convergence 
of our scheme by consideration of convergence of a functional series in 
transformed spatial coordinates. In transformed spatial coordinates 
we may use 

Proposition 1.5. Let K be a compact domain and assume that a 
functional series {hk)k '■ Cb{K) — )• M is given. Assume that for all 
X £ K the series ^^k&N^kixo) converges, and assume that for 1 < 
j < n the series ^^m^kj converges uniformly in K. Then the series 
X^feeN ^fc converges uniformly to a function h* which is differentiahle 
and such that for 1 < j < n we have 

hAx) = ^hk,,{x). (132) 

km 

The transformation uses the polynomial decay of order 5 described 
above. Application of this lemma then leads to a classical solution 
(locally in time) on a compact space Ki (the image of a coordinate 
transformation of DJ) and then to a classical solution on Df. 
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c) The third possibihty is to do the contraction estimate in a stronger 
Holder norm, i.e. 

\ r,p,k+l,l _ r,p,k,l\n \ ^ \ r,p,k,l _ r,p,k-l,l\n 

1^ ^ ll+a/2,2+al - ^ ll+a/2,2+Q- V'''^'^) 

and then establish decay at spatial infinity as in b). Then in trans- 
formed coordinates for the equivalent series (w'''^''^'')^. we have mv^'P'^^^ g 

Ci+a/2,2+a i^i), and w'''^'' E Ci+Q,/2,2+Q (Ki) as k-[ oo. Note that the 
latter space is a Banach space. 

We may summarize our construction as follows. A solution of the Navier 
Stokes equation system which is given by a global fixed point 

^^^rf = F{^), (134) 

where f = (/i,--- , fn)'^ is defined on [0, oo) x M", and V'^ satisfies the 
equation 

f dv{ _ d'^v{ r dv{ _ 

dt ^ ^j=l 'dx[ + ^j=l ■'J dxj - 

/m" {£-K4--y)) Elk=i {^^) ir,y)dy, 

I v/(0,.) = h(0,.). 

(135) 

This fixed point may be denoted by 

V* :=F(v*) = v^*. (136) 

It seems difficult to construct this fixed point directly. Therefore, we con- 
struct a fixed point of an equivalent problem which has the solution = 
V + r in successive time steps / on the domains [I — 1,1] x R", where 
r = (ri, • • • ,r„) is a bounded function. We solve for v'' in time steps with 

functions v'''^'' = • • • , Vn'''] , via a functional series 



^r,p,l ^ ^r,p,0,l ^ Svr,p,k,l_ (-^37^ 
k>l 

We establish a contraction property for the specific functional series v^'P^^^^ 
and prove convergence to v^'''''. The function r equals r' on each domain 
(/ — 1,/] X M", and is defined recursively in such a way that the solution 
is bounded and r bounded globally in time. Here, solves at each 
time step a linearized equation of Navier-Stokes type with 'consumption' 
source term. Note that this is crucial for global convergence. The local 
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contraction property for the higher order corrections fohow from a priori 
estimates of fundamental solutions of approximative subproblems. We shall 
use the Levy expansion of scalar linear fundamental solutions. We may also 
use the following standard result. 

Remark 1.6. Note that the application of standard Schauder estimates (as 
will be used in the second part of this article) requires a certain sign of the 
coefficient of the potential term. This sign is not given for the successive 
problems for v^'P'\ I > 1. However, since the functions r' are bounded, we 
can always consider equivalent problems related to 

ul'''''\T,x) = e-^-vl'''\r,x) (138) 

with A > large enough have coefficients for potential terms with the same 
sign as the time derivative. 

This is also a difference to our construction in [6] for the multidimen- 
sional Burgers equation. If the transformation (|138|) and standard Schauder 
estimates are used in order to establish a local contraction property of our 
scheme, then it is not sufficient to prove linear growth of the solution glob- 
ally. However, since we prove he boundedness of the solution we can even 
use standard Schauder estimates at each time step in order to prove local 
convergence. For the sake of completeness let us refer to this standard results 
which may be found in |llj . 

Let L be a scalar parabolic partial differential operator defined by 

Lf{t,x) = §if{t,x) - Elj=iaijit,x)g^f{t,x) 

(139) 

+ Er=i Ut, ^)irj{t, x) + c{t, x)f{t, x) 
for ah / G {Qt) for arbitrary T G (0,oo). 

Remark 1.7. Note that any potential term c/ has the same (positive ) sign 
as the time derivative. 

Let us assume that the spatial part of the operator L is uniformly elliptic 
with ellipticity constant A G (0, oo) (ellipticity from below) and A G (A, oo) 
(ellipticity from above), i.e. for all 

V C G \ {0} , ^ G M" : A|e|' <Y.a,j{t,x)i,ij < A|eP (140) 

Furthermore we assume that the coefficients have bounded Holder norms, 
i.e. 

I«ijl<5/2,5,QT ^ \bi\s/2,S,QT ^ \c\s/2,6,Qt ^ ^ (141) 
for some constant K > {). 
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Theorem 1.8. Assume that {l-jC^ o-nd jl-jl^ hold and that c < r for some 
constant r > ( note that minus sign in front of the second order coefficients 
in the definition of the parabolic operator L, cf. remark |j.7| jj. Let g G 
C2+'5(R") and f e C^/'^'\Qt). Then there exists a unique function u G 
(ji+&/2,2+S (^Qj,^ which solves the Cauchy problem 

Lu = f on Qt, 

(142) 

u(0, .) = g{.) onW. 
Moreover the estimate 

W\i+5/2,2+5,Qt ^ ^ (b|2+<5,M" + |/|5/2,<5,Qt) (143) 
holds for some constant C > depending only on n, A, A, 6, K. 

If we want to apply the latter estimate then it is not enough to show that 
the growth of v''''''' is at most linear. This is even more obvious as in the 
case where we deal with the weaker norm |.|i,2- However, using the stronger 
Holder norm leads to an additional difficulty with respect to the control of 
the integral magnitude ()12p . Instead, we shall show constructively that we 
have a classical bounded solution with an integral magnitude (jl2p which 
has linear growth with respect to time. This is achieved by a dynamical 
choice of the functions r-, i.e., a recursive construction with respect to the 
time step number I. In principal this allows us to apply the Schauder type 
estimates locally with respect to time. However, in this paper we use classical 
expansions of fundamental solutions. In case of the Navier-Stokes equation 
on compact manifolds (the flat n-torus is a simple example of a flat one) 
we do not have to deal with the integral norm seperately. However, some 
different difficulties arise if the coefficient of second order depend on the 
spatial variables. We shall treat this case in a subsequent paper. 

In order to prove the existence of a bounded solution v = [vi, ■ ■ ■ , u„) of 
the incompressible Navier Stokes equation in its Leray projection form we 
construct a bounded function r = (ri, • • • , r„) which is smooth with respect 
to the spatial variables and differentiable with respect to time except at 
a discrete set of points. Moreover the functions r — )• rj(r, x) are in H^''p 
for all 1 < i < n and p > 1, i.e. the first time derivative is in L^ for 
p > 1. The values where the time derivatives of r — )• ri{T,x) do not exist 
in a classical sense coincide with the time discretizaion of our scheme, i.e. 
with the values of natural numbers (r E N). The bounded function r is 
constructed recursively at each time step I. Hence r is determined by a 
sequence of functions = (r^,--- ,r^), where each is defined on the 
domain [I — 1,1] x R". Instead of solving the incompressible Navier-Stokes 
equation we solve a corresponding equation for the function 

(r,x) ^ v^''' := v'' + r, (144) 
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time step by time step, i.e., in the form v^'''' = v^' + r for / > 1 where 
each w^''P'^ is defined on the domain [/ — 1, Z] x R". We first determine from 
data of the previous time step, i.e., from v^'''''"^ (or from the initial data h 
if I = 1. This is done in such a way that the locally constructed solution for 
v''''''' stays bounded independently of the time step number /. 

This is the outline of this paper. In the next Section 2 we prove the main 
theorem, i.e., the existence of a bounded regular function r such that := 
V + r satisfies an equation equivalent to the incompressible Navier-Stokes 
equation in an (almost) classical sense. This leads to a classical solution v 
for the Navier-Stokes equation itself. Then in Section 3 we draw conclusions 
concerning regularity, uniqueness, and extensions to equations with external 
forces. Finally, in Section 4 we sketch an algorithm for the solution of 
the Navier-Stokes equation by iterative use of local analytic expansions of 
scalar parabolic equations. This algorithm simplifies the construction of the 
global solution via the function v^' since we can use the information that the 
solution is bounded. 



|a:h(x)| < rr-^k (145) 



|9r9^fex(t, x)\ < , for all (t, x) G [0, oo) x M", (146) 



2 Main theorem 

We shall assume that the initial data function x — )• h(x) = {hi{x), • • • , hn{x))^ 
with components hi in C°° (M") for all 1 < i < n satisfies 

Cak 

(iTR 

for all X G M", and for all multiindices a and and integers k with some 
constants Cak- 

Remark 2.1. If there is an external force term i^x included on the right then 
we may assume that 

Camk 

(1 + \x\+t) 

for any multiindices a, integers /c, and nonnegative integers m, and with 
some constants Camk- 

Essentially, the assumptions on the external forces ^ex and the initial data 
h mean that these functions are located in Sobolev spaces of arbitrary 
order s S M, i.e., for all s G M we have 

h G [H' (R")]" , (147) 

and for all t G [0, oo) 

fe..(t,.)G[i/^(R")r. (148) 
Here, for s G M is an ordinary Hilbert space with 

W = (M") = {/ G S' (M") 1/(0' (1 + < oo} , (149) 
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where / denotes the Fourier transform and S' (M") denotes the space of tem- 
pered distributions on M". In order to keep notation simple we consider the 
case fi = from now on (however, this is not essential and it will be clear 
that our construction can be generalized immediately to the case involving 
source terms satisfying ()148p or ()184p ). As indicated in the introduction 
our proof of a bounded regular solution of the incompressible Navier-Stokes 
equation consists of three main ideas: a) we introduce a time discretization 
and a series of linear time transformations t = piT such that time step size 
1 in r-coordinates is related to a small time step size in original coordinates 
and small coefficients of spatial derivatives in transformed time coordinates 
such that an iteration procedure leads to a local solution in time; b) in order 
to control the growth of the solution we introduce a locally regular func- 
tion r and solve an equivalent problem for the function v*" := v -|- r at each 
time step where r is itself a bounded function. The function r solves a lin- 
earized Navier-Stokes type equation, but with a 'consumption' source term 
which has been explained in the introduction to some extent. This source 
term consists of n components 0' = ,</>^) such that the functions 

■^r,p,i g^^^ J.I ijoTjnded independent of the time step number I. Here, the 
'consumption' function (j)^ and the function are defined recursively on the 
domains DJ of each time step number I and depend on the solution v'^'^"^ 
and r'~^ defined at the previous time step or dependent on the data h if 
Z = 1; c) At each time step we ensure that the solution has some decay at 
spatial infinity which is necessary in order to prove the existence of a time 
local fixed point. In the case n = 3 it is sufficient to have f[''''^(r, .) € 
for alU > 1 and r G [1-1,1]. 

Remark 2.2. In the introduction of this paper we mentioned that global 
regular solutions of the multivariate Burgers equation can be obtained by 
a priori estimates of the form ([3]). If we consider the method (i) described 
above and start at each time step with the solution of the associated Burgers 
equation with a certain source term then we need to prove that solutions of 
the associated Burgers equation satisfy a certain decay at spatial infinity (in 
case the initial condition at each time step satisfies a certain decay condition 
at spatial infinity). This is true in any case a), b), and c) of local fixed point 
construction above. For example, if we consider the method c) of the fixed 
point construction, then a decay of polynomial order 5 at spatial infinity is 
sufficient in order to prove the global convergence of our scheme to a classical 
solution. Hence, in order to apply method i) for the local iteration part it 
is sufficient to prove the following: 

Proposition 2.3. The Cauchy problem for the multivariate Burgers equa- 
tion with initial data h G H'^ for s G M has a classical solution u G 
[Ci+a/2,2+a ([0, oo) X M"))" such that u(t, .) satisfies il83\) for k < 5, i.e. 
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we have for all t G [0, oo) 

Cak 

(1+1^ 

More over the solution is bounded by the initial data, i.e., ^ holds. 



I^X*' •)! < , fork<5 and \a\ < 2. (150) 



The proof of this proposition essentiahy follows from our considerations 
in [6] and from estimates of type ([3|). The result in Proposition 1 can be 
sharpened in order to construct more regular solutions. In this Section we 
deal with the construction of classical solutions (which is essential). How- 
ever, since we use the alternative (iii) in this paper we do not use this 
proposition in the present proof. 

Next let us consider the construction of the function r more closely. 
Assume that we have solved for v^'P'^ ;= v^'™' + r"* for 1 < m < / — 1 
such that the solution for v^''' is constructed on the domain [0,1 — 1] x M". 
Especially the functions r™ for 1 < m < I — 1 have been constructed. 
Consider the restriction v'''''''~^ of v'"''' to the domain [I — 1, /] x M". Equation 
()24p may be written explicitly in the form 



En U-V-- sr^n r.p.tov.' n 

,=1 + PI E,=i = V'i, 

v'''/'''(/- 1,.) = v^'^''"i(/ - 1,.), 



(151) 



where 



(152) 



,r.p,l-l 



+PlE]=irf +P/E"=i^^i''''' (153) 



+Pi /m" {-ir^nix - y)) Elk=i 

In order to control the growth of the functions v^''''' at time step I we control 
1) the growth of v^'P'^^^ which satisfies a linearized equation, and 2) ensure 
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that the correction v''^''— v^'^''''' = Xlfc^i ^v^''''*^'' is smah enough by choosing 
pi appropriately. Let us look at the linear approximation v^'^''^'' of the 
function v^''''' first. The equation for v'"''''"'' is 

(154) 

+pi [-tM^ - y)) ^lk=i {^^'-^) - ^^y)dy, 

Note that all functions with index I — 1 are defined on the domain [/ — 
2, / — 1] X of the previous time step. Therefore these functions are alway 
evaluated at time / — 1 if they occur in the /th time step. Considering the 
linearized equation for v''''''*^'' has the advantage that we may represent the 
solution of ()154p in terms of the fundamental solution of the equation 

j=l 3 j=l J 

This representation involves a spatial integral with the initial data and a 
source term integral with respect to space and time involving the right side 
of the equation ()154p . The classical representation involves an integral with 
initial data which is solution of an equation where we can apply the max- 
imum principle. The second term can be controlled by choice of r' as in- 
dicated in the introduction. Next determining r' in a straightforward way 
such that the right side of p54|) becomes a function approximating the 
function -i/'' (|15ip would lead us to a Navier-Stokes type equation for . The 
equation is 

^ I ' ■ 



A I I V^"- I 9r'. 



-PI /m" {m-Kn{x - y)) Ei,fe=i (Sjfe) {r,y)dy 



(156) 



,1,0 



r'(/-l,.)=r'-i(/-l, 
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However, solving for r in this form would imply that we have transferred 
the original difficulties to control the growth of the solution to r'. Instead, 
we consider a different right side which approximates vl'*^ for Pi small. 
Here, 'approximation' is in the sense that ■0'''^ = 4>\+ {i^^i^ ~ 4>^i^ > and the 
difference becomes small if pi becomes small (cf. proof of main theorem be- 
low). Furthermore ?/'-''^ approximates Here, we let r\ solve a linearization 
of ()15ip . i.e. the function solves 

< +Lf''0(r'-i(/-l,.);v'^W-i(/_i,.)) 

-PI {-tMx - y)) Elm=i [^^'-^) (1 - ^,y)dy + 

, r'(Z-l,.)=r'-n/-l,.)- 

(157) 

Note that it is the same r' which solves this linearized equation which is the 
equation with a different source term Hence the difference of ^'"'^ and 
(f>[ may be obtained by subtracting equation (|157p from equation ()156p . We 
get 

-PI /r" {£-^n{x - y)) Elk=i (Stfe) 

+Pi /r" {£-^n{x - y)) Elk=i (%- 1) " ^'^^^2/ ^'^^^^ 
+Lf'''° (r'-r'-i,v-'^''-i) 

We shall indeed see that we can make the right side of ()158p small indepen- 
dently of I with the right choice of pi and r'. How small do we want to get 
it? Well, we want to have 4)\ close to il^^^ and ijj^^ close to -0' where we want 
to choose such that it controls both, v^''''' and r'. This is done as follows: 
first we choose based on the dynamic information of v'''^''"^ and r'~^. 
We do this as indicated in the introduction, and as we explain in a more 
elaborative way now. Further details are provided below in the proof of the 
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main theorem. For each I < i < n and at each time step / the function 
is constructed as a sum 4>[ = (j^l'^ + (j)/ where (jf'^ is designed in order to 
control the growth of v^''''"'' and v^'^'' and is designed in order to control 
the growth of For all 1 < i < n the function t/)^'' is constructed such 
that 0^'' gets the opposite sign of v['^'^ ^(Z — 1, .) whenever the modulus of 
the latter function exceeds a certain level. Assuming inductively that for 
1 < i < n we have 

r,pJ— 1 / 



sup 1,.)| < C (159) 



for some C > for all 1 < i < n we let D^j^^ ^ be the set of coordinates where 
^np,'-i exceeds the level C/2, i.e., we have = — 1, x) E 

and we let D^'^r^ be the set — G \—C, — ^] |, i.e., the set 

of arguments where t;'"'^''^! is below the level —C/2. Then we define 

(/.^'^ (/-I,/] xE"^M, <pf{T,x):=^l'''\x), (160) 

with a spatial function (jh'^'^ G Ch"^ using partions of unity such that 



6.' {x) 



-1 if X G 



1 ifxGD"''"', 



(161) 



and such that 



Later, in the proof of theorem 12.41 we shall see that 



sup \cl)r'\^)\ < 1- (162) 



10,2 



< Cs, (163) 



for some constant Cs independent of the time step number I. The construc- 
tion of (j)^'^ will ensure that the function ^r^''P'^ = v^'^ + r by the supremum 
of the initial data -v'^'f'''~^{l — 1, .) on the time interval [/ — 1, At the same 
time we have to ensure that the function r' is bounded. Otherwise we cannot 
ensure that v' := v^'^'' — is bounded. For the construction of we may 
assume from the previous time step that for 1 < i < n we have 

sup |r^H^-l,-)l <C°, \r'-\l-l,.)\i,2<Cr, (164) 

for some C > 0, and where C^,Cr are constants independent of the time 
step number I which will be determined below. Recall that we defined sets 
D^''^^ and D^_^~^ in order to define r' such that its growth is controlled. The 

definition of involves subsets C D^'^J^ and D'^^ ^ C -D^',- ^ sets 

such that { Dl\r^ U Dl\r^) n {0%^-^ U 0':^}-'] = 0. Defining of-^ = 
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1 y D'^'-, 1 latter relation is easily seen by rewriting these sets as 
subsets of 



and 



D 



r.O,«-l 



D:;y-':=\x\rlil-l,x)G 



C 



c 



x\r: 



{l-l,x)e [-C,- 



c 



x(^Dl 



Recall that we defined 



<^['':(Z-l,/]xM'' 



with a spatial variable function (/>[''^'' G CjJ'^ such that 



-1 ifxGZ)[;J'-\ 
1 if X G 

We shall see in the proof of our theorem 12.41 that 

sup \4>Y'\x)\ < 1, \fi'''\x)\o,2 < Cs, 



(165) 
(166) 
(167) 

(168) 
(169) 



where the constant Cg is as above. Having constructed = + (j)'/ for 
all 1 < i < n we solve the linear equation ()157p for Note hat all terms 
on the right side of (jl57p have the factor pi except for (/>-. Choosing pi small 
enough we can ensure that 4'\ dominates the source terms in the critical 
regions of arguments where the moduli of the function r'~^ exceed a certain 
level. This will allow us to control the growth of the function We shall 
see that classical a priori estimates ensure that 



kill,2 < Cr 



(170) 



for some constant Cr independent of the time step number I. Then we plug 
in this function r' into the right side of (|154p . and using (|158p we get 



En 
7 = 1 



dv- 



r,p,0,l 



4' = 4>l + {4' - 4 



(171) 



.[''''°''(/-l,.) = <'''''-^(/-l,.). 



Then we may represent the solution of ()17ip in terms of the fundamental 
solution of the equation (|155p and get 



vI''''^'\t,x) = ff,nvl'^'^ ^{l - l,y)T^{T,x;0,y)dy 
+ /r" ('^i + - is,y)r\T,x;s,y)dsdy. 



(172) 
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The first term of this representation can be estimated using the maximum 
principle for the corresponding equation without source term. The second 
term can be estimated (and becomes small) if piCr becomes small (indepen- 
dently of the number I). Finally we have to ensure that the 'higher order 
corrections' of the functions v^''''^'\ 1 < i < i.e., the functional series 
'^'kLi dv'^'P'^'^ which we have to add in order to compute the function v^'P''- 
from its linear approximation v'''^''^'', becomes small. 
We prove 

Theorem 2.4. Given any dimension n and a viscosity constant u > let 
h satisfy \14'T\) for any s G M. Then there is a global classical solution 

V G [Ci'2([0,oo) X M")]" (173) 

to the Navier-Stokes equation system ^B^, ^ which satisfies 

v^{t,.)eH^ (174) 

Remark 2.5. The proof below also shows that for all t > and all 1 < i < n. 

•)! - ^^^"^i^fc for A; < 5 and |a| < 2 (175) 

for some constants Cak > which depend on dimension n the viscosity con- 
stant 1/ and the initial data h. Moreover, below in Section 4 we shall use the 
proof of theorem 12. 41 in order to improve the scheme above for numerical pur- 
pose and in order to adapt it to initial-boundary value problems. Especially, 
the function r is not needed for the algorithm since the the boundedness of 
the solution v is shown. The main purpose of the function r is to show 
boundedness of the solution v. Once this is done it is not needed anymore 
and the algorithm simplifies. 

Proof. We do the proof in four steps. 

1) In a first step we prove existence of a bounded local classical solution 
v^'' assuming existence of a regular function ■vP'^~^{1 — 1,.), i.e., we 
prove existence of the local solution v^'' on [/ — 1, /] x M" for / — 1 > 
where bounded data wf' ~^{l — 1, .) G C| n are given. Here C| = 
C|(M") denotes the space of functions which have bounded derivatives 
up to second order. In this first step we set r = 0. Recall that 
^P,i — y^r,p,i _ ^ ^]^^^ ^p,/ _ ^o,p,i _ ^r,p,/ -j^ i^g^gg r = 0. Moreover, 

we show that the local solution has certain polynomial decay at spatial 
infinity (as indicated in the statement of theorem 12.41 above) if the 
data have this property. Especially, we show that for each component 
function v\'^'^ of the function v^'' we have 

<'^''(t,.) Gi/2 (176) 

for ah r G [/ - 1,/]. 
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In a second step we show that the local iteration of the first step of this 
proof can be extended in order to show existence of a local solution 
v'"'^'^ on — 1,/] X M" for a certain class of functions assuming 
existence of a regular function — 1, .) and a bounded regular 

function with a linearly bounded integral magnitude for r' (as in 
()12p ). i.e., we show existence of a bounded local solution v'^''''' on 
[/ - 1, X M" for I - 1 > where bounded data vl'^'^~^{l -1,.) eC^ 
and r^"^{l — 1, .) G are given, and where 

I E|(gfi)(*. i-n 

Moreover, we show that the local solution has polynomial decay stated 
in theorem 12.41 if the initial data of the respective time step have this 
property. 

In a third step we show that there exists a recursively defined bounded 
function r and a global iterative scheme of the solution v^'^ = v + r. 
Especially, for each time step number / the bounded function r' is 
defined via a linearized Navier-Stokes equation along with a source 
term = [(/)[,■■■ ,4'n) (among other source terms). These source 
terms (f)[ are defined in terms of the functions v^'^''~^(/ — 1,.) and 
r^~^{l — 1, .) which are known from the previous time step. They are 
defined in order to control the growth of the function v^''^. Here we use 
the fact that we have some freedom in order to define the functions 
rl- Indeed these functions are defined in terms of the 'consumption 
terms' 4>K We shall conclude that the functions v'"''''^ and r are globally 
bounded, i.e., there exists a constant C > independent of / such that 
for all 1 < i < n and / > 1 we have 

< C. (178) 

Moreover at each time step / > 1 we shall ensure that 

\rl\o<C, \vP\ = \vl'P''-rl\o<C. (179) 

The construction of the function r which equals at each time step 
I is such that a) step 2 above can be applied with pi and r' as they 
are specified in this step, and b) such that the functions r' and their 
spatial derivatives up to second order on (/ — 1,/] x are bounded 
and have a finite integral magnitude (|177|) . The integral magnitude 
increases linearly in time. Hence the choice ~ j in order to have 
convergence of the local iteration scheme while the scheme is global. 

In a fourth step we show the existence of a globally bounded classi- 
cal solution of the Navier Stokes equation, i.e., we show that = 
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+ r G 



+ r' E 



[Ci'2([0,oo 

y1,2 



C;'^([/-l,/]x 



and that for alH > 1 we have v'''^ = 
accordingly. This is not an im- 



mediate consequence of step 3) because for all x G M" the function 
T — >• r(r, x) is only weakly differentiable at the integer values / G N 
in general. Since the sequence (pi) is chosen such that J2iLi Pi t ^ 
as t oo we conclude that the global solution v in original time 
coordinates exists. 



2.1 step 1: proof of local existence of solutions at each time 
step 

The proof of step i) is in two substeps. First we prove a contraction property 
on the domain [1 — 1,1] x R" for the series {v^' ' )fe, < « < n with 

k 

m=l 

/I 12 

for a fixed time step / where we assume that v^' [I — 1,.) G Cj,' n 
at the /th time step. In a first substep we construct an iteration in 
Cj''^ {[l-l,l]xW) which is not an iteration in a Banach space. However, 
it is useful to have classical representations of approximating solutions in or- 
der to do some estimates. The limit function is Holder continuous if we have 
have a certain decay at infinity. Therefore, in a second substep of this first 
step we show that we have a certain decay at spatial infinity such that we 
can use one of the method a) (in case n = 3 with Hilbert space theory, i.e., 
showing convergence of the series v^'^' ' in the Hilbert space , or showing 
convergence of the series v^''''^''' in with respect to H^''^ Banach spaces in 
the general case ) or the even more elementary method b) (where n is an 
arbitrary dimension n) in order to prove the existence of a local fixed point. 
The method c) of the introduction is considered in a subsequent paper- it 
becomes important if one considers the Navier-Stokes equation on manifolds. 

Remark 2.6. Note that the argument of this first step is local. We establish 
some property for the function v'''^ assuming that some property holds for 
the function v'''^~^(l — 1, .) of the previous time step which serves as initial 
data for the Ith time step, i.e., v''''(Z — 1, .) = v''''~^(/ — 1, .). This means that 
the argument here is about local existence at an arbitrary time step I. The 
global existence (induction over /) is then considered with the introduction 
of the function r in the third step of this proof. 

For I > 1 assume that v''''~^(Z — 1, .) has been computed, i.e., the initial 
data of equation ()8ip have been computed, and assume that for all 1 < i < n 

(181) 
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for some constant Cg ^ , C[ ^ , Cg 2^ , C[ 2 which are given from the previous 
time step I — 1. In order to get estimates on Hilbert spaces we use an 
inductive assumption for the H^-novm, i.e. we assume that 



EE 

1=1 \a\<2 



Ivf't'H - l,x)\'dx < Ci,' + l)Ci~2\ (182) 



The reason for this bound -which is hnear in time- wih become apparent in 
step 3 of this proof. 

Furthermore, in order to apply the alternative method b) described in 
the introduction we assume that initial data function of the previous time 

step X — )• v^''~^(Z — 1, x) = ^t^i''~^(? — 1, x), • • • , Vn''~^il — 1, x)^ with com- 
ponents w^''""*^ which for all 1 < i < n satisfy 

\dy^'^'\l-l,x)\< (183) 

for all X G M", and for all multiindices |a| < 2 and and integers 1 < A; < 5 
with some constants Cak- 

Remark 2.7. If there is an external force term included on the right then 
we may assume that 

\drd^{eAt,x)\ < ^^_^f^jl^^k for ah (t,x) G [0,oo)x]R", (184) 

for any multiindices a, integers k, and nonnegative integers m, and with 
some constants Camk- 

In case I = 1 (the first time step) we have ^rP'^~^(l — 1, .) = h and the 
constants Cq, Cj*, Cg 2 are given in terms of upper bounds of the components 
of the functions hi and derivatives of the functions hi. The series ()180p can 
be generated by an iterative application of the map 

FrA^yrf'P'\ (185) 

starting with the value of Fi applied to the initial data of step /, i.e., starting 
with the function v'''*^'' = Fi {yP'^~^{l — 1, .)) . Here, — 1, .) represent 

the final data of the previous step or the data h in case 1 = 1. The function 
^f^,p,l ig determined by the solution of the equation 



,=1 + PI E,=i fj 



Pi /r" (£-Knix - y)) Elk=i (m^§k) (186) 
rf'P'^il - 1,.) = vP'^~^l - 1,.). 
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We want to show that the map Fi is a contraction on this series with respect 
to some appropriate norm, i.e., we want to show that 




IJv^.'^.'l < ij^v^.'^-i.'l (187) 

for some suitable norms |.| and /c > 1, and where we start the series with 
^pfl,l — ^p,l-i'^^ The norms are the |.|i^2 norm in order to have classical 
representations of approximating solutions in terms of certain fundamental 
solutions, and the |.|//2- and |.|//2,p-norms in order finish the estimate for 
the limit according to method a). Furthermore, we shall show that we have 
some decay at infinity. We shall also see that we have a linear growth of the 
integral magnitude for v'''^, i.e., we have a linear growth upper bound with 
respect to time of the magnitude 



{t,y)dy. (188) 



Since products of functions may be pointwise estimated by half the sum of 
the squares of the factor functions, the linear growth upper bound of (jl88p 
follows from the linear growth upper bound of the i/^-norm of the solution 
function. Note that the term (jl88p is related to the integral term of the Leray 
projction form of the Navier-Stokes equation naturally. Since ()188p can be 
estimated quite naturally in term of norms, it is quite natural to estimate 
the solution of the Navier-Stokes equation in the norm. Interestingly, 
in case of n = 3 this is also sufficient in order to have an embedding in 
Holder continuous spaces. Recall that in case / = 1 we have v''''"^ = h. 
In order to prove a contraction property of type ()187p we consider the map 

Fi{i) = v^''''' for general f G C^J'^ first, where v^'^'' satisfies the equation 

()186p . We assume that the data ?;f'' ^{l — 1,.) G CjJ'^ are bounded (this 
is certainly true in the case 1 = 1, where v^'' ^(i — 1, .) = h). Note again 
that in this first step we let r = 0. We consider the methods a) and b) of 

1 2 

the introduction and consider first the normed space C^^' (DJ) along with 

12 

the norm |.|i.2 of the introduction. Note that this normed space, i.e., C^^' 
equipped with the norm |.|i^2 is not a Banach space, but we can proceed 
according to the methods a) or b) of the introduction. We shall determine 

^p,i — ^yP'^ ^ . . . ^ as a limit of a functional series with members in C^'^, 

and with an additional decay at infinity, i.e., for all 1 < i < n we construct 
uf'' as a limit of the approximations f^''^'' 

C'/ {DJ) 3 vt' = + X: Svr\ (189) 
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where we prove in a second substep that all functions f^''^'' have a certain 
decay at infinity. Here, the function v^'*^'' = (^v^'^'^ , ■ ■ ■ ,Vn^'^^ solves 

-57 Pi^ 2^j=i Q^i + Pi 2^j=i 

= PI - y)) E (^^) ir,y)dy, 

^ vP^^'\l-l,.)=vP^'-\l-l,.). 

Furthermore, for A; > 1 we consider the difference functions 5^^'^'^ = v''''^'' — 
^^p,k-i,i -^hgj-e for all k we have v''''^'' = v'"'''''^''|r=o- Note that we have 6vP'^'^ 
recursively defined to be solutions of the equations 

- pi^ Ej=i -9^ + PI Ej=i ^s— = 



/m^ - ,)) E • .=1 (^^^^J (191) 

j^p,fc,;(; - 1,.) = 0. 

Here, for /c = 1 we denote v'^'^~'^'^ = t'^'""^'' := v'^'^~^ . Solving these equations 
recursively leads us successively to the functions 



k 

rP,k,l _ ,.p,0,/ 



v''''"'^ (192) 



m=l 

which approximate the local solution v''''. In order to apply the methods a) 
and b) of the introduction we show first that 

\S^r'''''\2 < \\S^rP'''-''^\l^2. (193) 

Then in order to apply method a) of the introduction we show that for all 
T £ [/ — 1, /] we have 

\S^r''^'Hr,.)\H^<\\S^'^'~'■{r,.)\H^. (194) 
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uniformly in r. We shall also consider how polynomial decay at infinity is 
preserved by the local scheme if the initial data from the previous time step 
have this property. In this work on the Navier-Stokes equation system we 
use methods a) and b) of the introduction in order to show the existence of 
time-local solutions for general dimension by elementary means. Recall that 
in case of method b) it is shown that the members of the functional series 
(v^'*^'^)fc satisfy a specific polynomial decay at infinity such that a spatial 
transformation leads us to a functional series in a Banach space, and this 
method applies for any dimension n. Hence we have two different proofs in 
the case n = 3 corresponding to the methods a) and b) of the introduction. 

Remark 2.8. Locally a sharper estimate with respect to the Holder space 
Ci+a/2,2+a (Df) IS possible. However, it is easier to control the growth with 
respect to the |.|i^2 norm. This control of the growth is essential for the 
global existence proof. Classical Schauder estimates introduce a potential 
term with a certain sign, and this would lead to a exponentially decreasing 
series of time step size numbers pi if applied naively. In this case the sum of 
the time step size numbers pi is finite and the scheme is not global in time. 
Hence, method c) of the introduction may be applied locally, but a globally 
some additional work is needed. 

Note that in case n = 3 it is essential to show in the second substep that 
v^''''\t, .) G (M") for alll < i < n, because the limit vf'\T, .) € (M") 
for all 1 < i < n (uniformly in r) implies that 

vP^^{T,.)ieC''cH^{M.^) (195) 

for some a > uniformly in r for all 1 < i < n. 

Remark 2.9. In order to extend the proof to general dimension n we may 
use the estimate For all s G M" and p G (1, oo) 

for r = s — ^. Hilbert space estimates for n > 3 are harder to obtain. The 

relation v'^'^'\t, .) G (M") for some s > m + ^ + a and some m > 0, and 
some a G (0, 1) for all 1 < i < n implies 

i;f'(T,.) G C i7"(M"). (197) 

Concerning method b) and c) of the introduction we note that the spatial 
transformation on a finite domain does not preserve the uniform ellipticity of 
the operator (represented by the positive constant v in the original equation). 
In this respect the method a) is a little bit more elementary than methods 
b) and c) of the introduction, because we do not need a spatial coordinate 
transformation in this case. In order to check the contraction property of 
the functional series (v''''^'')^>q we check the contraction property of the 
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functional ()185p on an appropriate functional subset of the function space 

" 12 1 " 

Cjj' (M") . First let us stick with the original spatial coordinates and con- 
sider the equation for ti/''''' more closely. Using relation (jlOSp from equation 
(fTnH]) we get 

= -Pi /iRn E"=i ifj - 9j) {s, y)-%^{s, y)^){T, x; s, y)dyds+ 

JLi Pi /r" /r" ^riA^ -y)^ (198) 

Yllk=i Uj,k - 9j,k) (s,y)) ^T^fiT,x;s,z)dydzds, 
where Tj- is the fundamental solution of the equation 

^ -AT^f + Pi^fj^ = (199) 



X 



on the domain DJ . 

Here, and in the following we write {fj — gj) {s, y) := fj{s, y) — gj{s, y), 
{fj,k - 9j,k) {s,y) = fj,k{s,y) - gj,k{s,y) etc. for simplicity of notation. We 
show that pi can be chosen such that the iterations 

^rK''~^) G 5*0, (200) 

1 2 

where 5o is a subset of C^^' , i.e., 

f,g G {f||/ik2 < 2C{7/&|/i(r,.)|H2 < 2C7'-i unif. Vr G 1,^]} =: Sq, 

(201) 

along with 

\v?''-\2<C[;l (202) 

The choice of the space 5*0 is related to the fact that we use classical repre- 
sentations of approximating solutions which may be estimated in terms of 
Gaussian estimates. These Gaussian estimates are convolutions which are 
estimated in terms applying Young's in equality (cf. below). 

Remark 2.10. You may ask why we consider the contraction in the space 

12 

C^' with norm |.|i^2 at all, since this is not a Banach space. The reason 
is that classical representations are very useful when we establish a local 
scheme. This will become apparent in step 3 of this proof. 
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The factor 2 in the definition of 5*0 is due to the fact that the first 
order coefficients vary with the the local iteration number k (see below). 
The requirement that hi S in ()20ip is due to the fact that we need 
to establish a contraction with respect to the norm in order to apply 
method a) of the introduction. In order to get a uniform bound of this 
first order coefficients which are independent of the iteration number k we 
establish an estimate of the form 

Er=ii^/''''-^r'ii,2<iif-gi?,2- (203) 

Recall from above that we equip the vector valued functions of the classical 

1 2 

function space C^' with the norm 

n 

1=1 

(we may also consider the maximum etc. instead of the sum, of course, 
but such considerations do matter only if we consider implementations of 
associated algorithms; the reference to the domain DJ is not denoted if it is 
clear from the context for sake of simplicity of notation) . 

Remark 2.11. Note that the estimate ()203p is one step of our construction of 
a solution in order to apply the methods of item a) and b) of the introduction. 
Locally in time there is a refined estimate for Holder spaces, i.e., we can get 
estimates of the form 

Etl b/'"'' - ^r\w2,2+a < C|f - g|^+,/2,2+„- (205) 

Here, the Holder norm is 

n 

|f|l+a/2,2+a = I /i I l+a/2,2+a ■ (206) 
i=l 

Then the method c) of the introduction can be used in order to construct 
a solution at each time step I. The application of the methods a) and b) is 
described below. 



In order to prove the contraction property (|2U3|) it is essential to estimate 
! second c 
to estimate 



the second order derivatives of the difference vj'^'^ — wf''''', i.e., it is essential 



J,P,l _.,9,P,l _ L/'P'^ _ y9,P,l\ (207) 



,f,P,i _ „,9,P,l _ I ^J'P'I- _ ,,5<P<' 

dxkdx 

Note that we consider the case = for all 1 < i < n in this first step of 
the proof. 
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We consider the two summands on the right side of (|198p separately. 
Accordingly, let us define 



V, 



= -Pi I»n Ej=i ifj - 9j) (s, (s, y)r^(r, x; s, y)dyds, 

and 

= S{i~i) Pi /r" /r" Kn^i{z — y)x 

(209) 

( (Ej,fc=i ifk,j{s,y) + gk,j{s,y))j X 

(E",fc=i {fj,kis,y) - gj,k{s,y))^ ^T^f{T,x;s,z)dydzds, 



a,p,i 



(208) 



Since f , g G 



we know that the classical fundamental solu- 



tion exists in its classical Levy expansion form. Since the initial data are 
C^'^ we know that is in C^'^ (Df) for all 1 < i < n. 

Especially, {t,x) — )• Ej=i (/j ~ S'j) {'^i^)'^fj'\'^i^) is Holder continuous 
in the spatial variables uniformly in r G [/ — 1,^]. 

Hence, classical Levy expansion theory of fundamental solutions (cf. [3]) 
shows that the function v('^'^'^{t, x) — vf'^'^{T, x) has spatial derivatives up 
to second order and that for spatial derivative of second order of the first 
summand we have the representation 



-Pi /ru Ej=i ifj - 9j) (s, y)%— (s, yW^f,j,ki'^, s, z)dydzds. 

(210) 

We shall see below how to rewrite this with only one spatial derivative of 
the adjoint of the fundamental solution F^. Next concerning the function 
yf'P'i'"^ _ y9.p.^2 £^j, fixed time parameter s the integral function part 

X ^ Jjgn Kn,i{x — y) X 

( (E",fe=i (/fe,i(s,2/) +5fc,j(s,y))) X (211) 

Ej,fc=i ifj,k{s,y) - ffj,fc(s,y))) )dy 
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looks formally like the solution of a Poisson equation 

n 

- Ap = ^ {fk,jfj,k - 9kj9j,k) ■ (212) 
j,k=i 

It is indeed a solution if the right side is in L^, if pointwise products of 
functions of type x — )• fk{s,x) and x — )• gk{s,x) are in H^'^. Again, since 

\fk,jhk\ir,x) < \ (|A,/(t,x) + \fj,k\\r,x)) (213) 

it is sufficient that we have /a:(t, ■),gk{T, ■) G (even /^(t, ■),gk{T, .) G if"^ 
is sufficient at this stage) . We shall see that the functions which are produced 
by iterations of the functional Fi starting with f^'' satisfy this requirement. 
Therefore in the following we may assume that this requirement is satisfied 
such that the function in ()21ip exists. We shall consider this in more detail 
later in this first step of the proof, and see that it follows from our polynomial 
decay condition on the initial data and the properties of the functional 

Fi. The following fact is rather classical. 

Lemma 2.12. Let n > 3, and let f ^ (M"). Let K he the fundamental 
solution for the Laplacian. Then 

{t, x) ^ uit, x) := / fiy)Kix - y)dy (214) 

is a well-defined locally integrable function which solves the Poisson equation 
on M" with data f, i.e., 

Am = /. (215) 

Remark 2.13. If n = 2 then we need the stronger assumption / f{z) log \z\dz < 
oo and the theorem is still true. 

Hence, if fj.,gj G C"^'^, then the function ()211|) is indeed the solution of 
(|212|) . Concerning the regularity of the solution u of the Poisson equation 
we have another classical result due to Holder himself. 

Lemma 2.14. Let n,k be some positive integers. Let Q C M" be an open 
set, and let f G C^+°'{Q). If u is a distributive solution of 

An = /, (216) 

i/ien u G ^2+°+'= (O). 

1 2 

If fj,9k G C^' for 1 < j,k < n, then the first order derivatives fj^i,gk,m 
are Holder continuous. This implies that the function of equation ()21ip is 
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in C^~^" (note the derivative of the kernel K). Hence, the integral 

f, 0,1,2/ \ q,p,L2/ \ 
= /(/-I) Pi /r" /r" Kn,i{z - y)x 

( ( \ ^^^^^ 

(E",fc=i {fj,k{s,y) - 9j,kis,y))^ ^r^j:j f^{T,x;s,z)dydzds 

exists and represents the second order derivatives of the second term (|2U9|) . 
The equations ()210p and ()217p contain second order derivatives of the fun- 
damental solution Fj. Note that the function 

z Kn,i{z -y)x 

( (Ei,fc=i (/fcj + (s, y)) X (218) 

(Ej,fc=i ifj,k - 9j,k) (s, y)) ) dy 

is Holder continuous, hence the function v^'^'j^'"^ — ffj'^^'^ is continuous. In 
order to do the contraction estimate it is useful to rewrite equation ()217p 
and equation (|210p . 

Note that the function ()218p is solution of a Poisson equation with a 
right side which is Holder uniformly in r. Hence, since the function (|218p 
is in C^^" with respect to the spatial variables and uniformly in r, we may 
rewrite (j217p shifting differentiation from the fundamental solution to the 
function ()218p (shifting one differentiation is enough). This way we have 
only first order derivatives of the fundamental solution in the representation 
of vl 'p^''^ — vf'pl'"^, and first order derivatives of the fundamental solution are 
integrable and can be handled more easily. Note that the derivative of the 
fundamental solution is with respect to the x-variables and the integrals in 
(|210p and ()217p are with respect to the conjugate y variables. In order to 
shift the differentiation one may use a relation of the fundamental solution 
and its adjoint, i.e., the relation 

r^(r,x;s,y) =r)*(s,y;r,x). (219) 

where F^* is the adjoint of Application of derivatives of this relation 
(|219p is exactly what we need. However, let us show this in detail in terms 
of the Levy expansion of the fundamental solution (as a by-product, this is 
also an independent proof of (|219p ). First we shift the derivative Ki of the 
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kernel K in (|217p . This is no problem since the inner integral with respect 
to z represents a convolution. Hence, we have 

= - Pi /r" /m" Kn{z -y)x 

( ( \ 

(Ej,fc=i - 9i.k) (s, 2/)) ) r'^^^^p(T, x; s, z)dydzds. 

Remark 2.15. Note that in this representation we have even only first deriva- 
tives of the difference {fj^kis, u) — gj,k{s, y))- This may be convenient in some 
circumstances (if we want a sharper estimate) but it is not necessary in case 

1 2 

of our modest goal, i.e. the contraction on a subset of C^' with respect to 
|.|i^2 norm. 

Now let us show explicitly how we shift one derivative of Tj^p First, 
recall from classical theory that the fundamental solution Tj in the Levy 
expansion form is given by 

r'f{T,x;s,y):=N\T,x;s,y)+ [ [ N\T,x;a,OH'^,C; s,y)dadC, (221) 

where ^ 

N\T,x;s,y)= \ expf- J^"^' V (222) 

and (p is a recursively defined function which is Holder continuous in x, i.e., 

oo 

</.(r, x; s,y)=J2 iLiN%{T, x; s, y), (223) 

m=l 

along with the recursion 

(LiiVOi(T, x; s, y) = LiN\t, x; s, y) 

(224) 

— ^ f 9N' 

— Pi 2^j=i Jj dxj ' 

{LN%+i{t, x) := £ {LN\t, X- a, i))^LN\a, s, y)dadi. 
First we look at ()220p . We may write 

f,p,l,2/ \ g,p,l,2/ \ 

(225) 

= - So Pi /r« '^'(■5, y)r/,^,p(T, x; s, z)dydzds 
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along with a Holder continuous function S 

S{s,y) := Kn{z - y)x 



2 (E",fc=i + 9k,j,i) {s, y)) X (226) 

(Ei,fc=i - ffi,fc) {s,y)) )dz, 

which is uniformly Holder continuous with respect to the time parameter s. 
Now we may write 

= - !o Pi /r- y)N!m,p{T, x; s, y)dsdy 

+ £ Ir^ S{s,y)N^^ p{T,x;a,C) (E^i(-^^^')g(c^. y)) S{s,y)dad^dydzds. 

(227) 

We may do the partial integral term by term and get exactly the (second 
order differential of the) Levy expansion for the adjoint. Well it does not 
matter so much that we get exactly the adjoint. What matters is that we get 
an absolutely convergent sum where only first order derivatives of the kernel 
A^; are involved, since such first order derivatives are locally integrable and 
can be estimated easily. We may denote the result of this partial integration 

from equation (|210p we get 



for the Levy expansion sum by T^l* (which happens to be the adjoint) . Hence, 



= -Pi lo /r« E"=i ifj - 9j) {s, yWif\s, ?/)r^,fc,^(r, x; s, y)dyds 

= Pi 1(1-1) /r" Ei=i ifj - 9j) (s, y)vl^[i{s, yW)l{s, y; r, x)dyds 

+Pi /(/-I) /r" Ej=i Uj,m - gj,m) (s, y)vif\s, y)T^f*k{s, y; r, x)dyds. 

(228) 

Next, we estimate the integral with integrand F^*^ by a constant Cyj, using 
estimates for local integrability of the first order derivatives of the fundamen- 
tal solution in Levy expansion form. Here, we may use standard estimates 
of the form 

< exp f-XfJ^^y^] , (229) 

for some constants C, A and for |x — y| > 1 in order to have integrability at 
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infinity, and then use 

(r — s)'^\x — y 



N^{r-s,x,y)\ < _ ,.+1-2. (230) 



for a G (0.5, 1) on a domain where |x — y| < 1 in order to have local 
integrability. We can formulate have both standard estimates independently 
of f G 5o. We get 

s "''*'=' - vscrrir« • ''''' 

where 



C = sup Cf, (232) 

f&So 



and 



A = inf A/, (233) 



in order to Furthermore, we get 



|r^^^(r,a;;s,y)| < _ ^^^^^ _ ^|„+i_2^: 



(234) 



for some G (0.5, 1) and with some constant . This is locally integrable 
with respect to time and space. This implies that there exists Cr > such 
that ^ 

sup / / sup \T^l*\dyds < Cr, (235) 

(r,x)6D[ JK" feSo 

and such that ^ 

sup / / sup \Th*-\dyds < Cr, (236) 
{t,x)£Di Ji-1 Jr" feSo 

where 1 < i < n. Hence, since g G Sq, for the second order derivatives of 
the first summand ()208p we get 

I f,p,l,l g,p,l,l\ ^ 
I i,k,m i,k,m — 

Pi /(l-i) /m" E"=i l/i - 9j\oW^''''^\o,2Cr 

+Pi /ign E"=i \ fj,m - 9j,m|o|v^'''''|o,iCr (237) 

<PiIli-i)k^\^-s\hW'''X,2Cr 

< piC[--,'C^\i - g\l,. 
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Note that in this estimate the second Cr comes from the estimate of |v^'''''|q 2- 
From equation (j220p we get 



_J,p,l,2 ,o,p,i,2| 



< /(/-I) Pi Kn,m{z y) X 

E",fc=il/fc,i,* + 5fcj,»|o) X (238) 



<PiCi^n24C{72'Cr|f-g|^,2- 

Here the constant Ck is from the spatial integral with integrand Kn^m- We 
may assume without loss of generality that Cr > 1. The constant Ck 
depends on dimension only and may be subsumed by a constant C*. In the 
following we shall use the latter constant generically if we want to subsume 
all constants depending only on dimension. We summarize that equation 
HMD and (I2H8D gives 

n 

E - <;--lo < PiC*nC'rC[-:,'\i - gll,. (239) 

j,m=l 

Similarly we get 

El=i H'!^' - < Pic:clc[-;,'\f - g|^,2- (240) 

and 

_ ^a,P,i\^ < p,QC2c{72V - g|^,2- (241) 

An estimate for the time derivatives of v- '''''' and v^''''' can be reduced to the 
estimates (|239p . (|240p . and (|24ip via the defining equation for 

^ = P.- E •=! -^-Pi EU = 0- (242) 

Remark 2.16. Alternatively we could use defining equations 



j=i - PI Ej=i fj^^+ 
/m" (irK^i^ - y)) E - .=1 (Hf it) ir,y)dy 



(243) 



for the functions vj'^'^ and analogous defining equations for the functions 
vf'K However this would enlarge the estimate constant. In this case we get 
another factor Ck and another factor C\ ^ ■ The method of growth control 
described in step 3 still applies. 
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Summarizing we get 

< piC:C^C[;,'\i - g\l, (244) 

< pic:c^c[;,'\{ - g\l,. 

for some generic constant C* > depending only on dimension n. We may 
choose 



PI < TTZT^OT- (245) 



1 

4C*CpC| 2^ 
Then for f,g £ Sq we have 

|v/W - v^^'^'^l^^s < i|f - g|?,2- (246) 

Note that 

f,gG5o^v^'^''G5o (247) 

for this choice of pi. Furthermore, note that the constant Cr depends on the 
class Sq. Hence, we do not have a contraction on the whole space but on a 

1 2 

subspace Sq of C^' . Since the definition of So contains a factor 2 and we 
start with < C[~2 the iterations are F;'^(v'''^~^), k > I belong to Sq 

indeed by choice of pi. Hence, the contraction condition leads to an iteration 

1 2 

in a subspace Sq of C^' . However, note that this is not a Banach space. In 
order to have the functional series in a Banach space or in order to have the 
limit of the functional series v'''' = lim/;i-oo v^''^'' known to be in some regular 
space, we need some decay at spatial infinity. Let us summarize first what 
we have achieved so far on the way of the construction of a local solution. 
Iteration of the map Fi starting with the function v''''"^ leads to a series 
of functions vP'^''- = F^~^^ (^v''''~^) where the components of the functions 
wP'^'^ satisfy 

k 

^P,k,l ^ ^PAI + ^ ^^P,m,l ^ ^1,2^ (248) 



and where the increments satisfy a contraction property 

(249) 



\S^P,k,l\n^ < piC*ClC[-2^\6vP'''-^'^\l2 



< l\5vP'''^^'^\l2- 

The approximations uf''^'' can be represented in terms of fundamental solu- 
tion r^, where for A: > 1 T^ is the fundamental solution of the equation 

fi-p,.Ari + «X;<--g=0, (250) 

j=l 
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and for A; = the function Fg is the fundamental solution of the equation 



dr 



j=i 3 



(251) 



The choice of pi ensures that the first order coefficients in (|'25Up and (j25ip 
have a uniform bound. This leads to the estimate of integrals involving 
integrands by a uniform constant Cr such that 



ll,2 = + 2^ C'V'^' |i,2 < 2Ci_2 ' 



(252) 



m=l 



and, hence, to a functional series in Sq. This leads to the contraction prop- 
erty of the series (v''''~^) with respect to the norm |.|i,2! and we have 
(|193p especially. 

Next in order to apply method a) of the introduction we show that the 
elements of the functional series f^''^''(r, .) are in uniformly with respect 
to r S — 1, Z] and for all 1 < i < n. This means that we extend the estimate 
(j249p above to the space of functions f with components 

/. G Hf := |/i G Cl'^{Dl)\h{T, .) G if^ unif. for all r G [Z - 1, /]} . (253) 

Here 'unif.' is an abbreviation for uniformly. This is a good space in order 
to have a finite magnitude ()188p . since we have 



k=i 



) (t V) 



dy 



(254) 



dv 



p,i 



dx , 



.p.i 



dy < oo. 



if G . We choose the space since we have first derivatives of (jl88p 
involved in our estimation. 

Recall that method a) of the introduction applies in the case n = 3. Note 
that a converging series 



^rP''''\T, .) G [H' (R")]" for some s > a + -n 



(uniformly in r) implies that 



v^''(r, .) G [H-' (M")]" for some s>a + 
Hence, in case n = 3 we have 



(255) 

(256) 
(257) 
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for a G (0,0.5) uniformly in r. Hence, if we can show that ()255p holds 
(which implies (j256p . then we have bounded Holder continuous first order 
coefficients in the equation 

^ - PI- Y.U '-^ + P^ = 0' (258) 

and this means that a fundamental solution of ()258p exists. Then we 
have the representation 

v1'\t,x) = vf~^{l - l,y)T[,{T,x]l - l,y)dy+ 

Ii-i PI /m^^ {^^niy - z)) [^%) (r, .)) X (259) 

xr^(T, X, s, y)dzdyds 
and from this representation we immediately get 



G 



n 



C'/iDJ) . (260) 



Now let us finish the proof in the case n = 3 using method a) of the intro- 
duction. We have to show that 

yP'^'^ (r, .) G (M") , (261) 

where we show that the latter relation holds uniformly in r. First, induc- 
tively with respect to I we have 

v^'^-\l - 1, .) G Ci and ^;f''-^ (/ - 1, .) G (M") . (262) 

First we show that 

vf''''\T, .) e (263) 

uniformly in r. The reasoning for the first and second derivatives below 
requires only a little more work. We start with k = 0. Note that the equation 
which defines v^'*^'' has first order coefficient functions ff''~^(^ ~ 1) •) which 
are independent of time. Hence the fundamental solution Fq of the equation 

^ - PI- Ej=i -a^ + Pi Ei=i as^ = (264) 

exists. For this fundamental solution we have a priori estimates. Especially, 
we have 

|r[,(r,x;.,y)| <^-^^exp(-A^^) (265) 
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for some constants C and A (the constants A and C are used generically 
here). Next we have the representation 

uf'°''(r,x) = vf~^{l - l,y)r[,(r, j;;/ - l,y)dy+ 

Ili-i) PI /m" (KnAy - -)) (e • .=1 (^^) (r, ^i) X (266) 

xrQ(r, s, y)dzdyds 
Note that for r = / — 1 we have v^'^~^{l — 1, .) G //^ inductively. We get 

\vf\T,x)\ < 4„ \vP'^"\l - l,j/)l (^_(;5i))„/2 exp (^-aJ^^) 



(^^^)n/2 exp (^-AL;-^ J dzdyds 
=: (/) + (//). 

(267) 

Now (I) and (II) are convolutions. We may apply Young's inequality, i.e., 

|/*5|r < l/Mfflp, (268) 

where 

1 + - = - + -, (269) 
r p q 

and 1 < p,q,r < oo. Note that the family of Gaussian functions 

^ exp (-X^^) (270) 



(r - s)"/2 

(family with respect to time parameters) is in for all 1 < p < oo (although 
not in L°° as s t r). We know that the first term (I) is in because 
vf' ~ G inductively and the Gaussian is in L^. Mote that we may 
set the function in ()270p to zero if r = s. Furthermore, applying Young's 
inequality for fixed t > s first we see that the integrand for the second 
term is in where v^'^"^ G inductively. Since the L^-bound of the 
integrand is uniform the time integral (//) is also L^. Here again, for each 
X G M" you may use decompose the integral in two parts [0, r] x Bx and 
[0, r] X M" \ ([0, r] x Bx), where Bx is a n-dimensional ball around x of finite 
radius. Then on the ball we may use the estimate 

' .expf-A^)< , ,271) 



(r-s)"/2 V T-s y ~ (t - s)^'|x - 
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rUr,x;s,y)\ < ^ ^ exp f-A^^ ) . (272) 



for /X G (0, 1). Outside the ball we may use the estimate 

C ( Ax-y 

Hence we have 

vf'^ (r, .) G l2 (M") (273) 

uniformly in r. Next note that we have representations of the first and sec- 
ond spatial derivatives of -y^'*^'' which involve only the first spatial derivative 
of the fundamental solution (for the second spatial derivative we use the ad- 
joint and partial integration). Then an analogous argument as above which 
uses the first derivative estimates 

« 1 expf-A^')< , ..,„..^.,..,_,„ (274, 



9xi(r-s)"/2 'T - s y ~ (t - s)/^|x - 

for fi G (0.5, 1), and 

|r^(r, x; s,y)\< _ exp ("A^;^) • (275) 

leads to the conclusion that 

vf'^ (r, .) G H'^ (M") (276) 
uniformly in r. Next assuming that 

^p,k-l,l ^ ^1,2 ^^^^ _^ ^ ^2 (j^n^ (277) 

uniformly in r G [/ — 1,1] we know that v'^'^'^ G C^'^ and we want to show 
that 

i;f'^'' (r, .) G ^2 (R") (278) 
uniformly in r. It is convenient to use the series (v''''^'')^ along with 

k 

^p,k,i ^ ^p,o,i ^ 5^p^rn,i. (279) 

m=l 

Then it suffices to show a a contraction property with respect to the H^- 
norm in the sense that 

\6vP'^\t,.)\h^<\\6vP^^-'\t,.)\h^ (280) 



for all /c > 1. This holds if we use a constant pi with 

1 

4C*Cp(C{ 2^ + IC[ 2 



PI < -^^,-^2.^z-i^,^z-u - (281) 
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We consider the main steps in order ro prove this relation. Note that the 
equation which defines v^'^'^ has first order coefficient functions v'^'^~^'^ g 

12 / 

' . Hence the fundamental solution of the equation 

-a^-Pi'' Ei=i -a^ + Pi Ej=i v^' &7 = (282) 

exists. Since we have \vj'^~^'^\i^2 ^ 2C( 2 independent of the number k. For 
this fundamental solution we have a priori estimates. Especially, we have 

\ri{T,x;s,y)\ < _^^)„/2 (-^^^) 

for some A, C > independent of the number k. Next we have the repre- 
sentation 

5v^'''\t,x) = -ft /i8„ (E"=i'^^f^~^'' ^"fe,"''' ) (s,?/)r^(T,x;s,y)ds(iy+ 
xr^(T, X, s, y)dzdyds 

- /(I-l) /m^^ /'^ /m« {-tKn{y - Z)) (e • ™=1 -)) X 

(284) 

This leads to 
{{Elm=i {v?if''is,y) + v?i[''{s,y)))x 

(285) 

Hence, using a priori estimates analogously as above in case k = we get 
(|28U|) . This finishes the proof in the case n = 3. The method a) as it 
is presented above applies in the case n = 3. In case n > 3 we have to 
ensure that the convergence of the series v^''^'\t, .) is in H^'^ for p large 
enough. This can be done using the same a priori estimates. Method b) 
provides an alternative which may be even more elementary. This direct 



{r.z) 
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(more elementary) way shows that the hmit of the functional series exists in a 
Banach space by transforming on a compact domain. Such a transformation 
can be applied if for each k > we have 

|5>^''=''(t,x)|<^^^^ifH< 2 . (286) 

for |a| < 2 and A; < 5. Then we may consider the transformation of spatial 
variables 

Xi — )• Vi ■= arctan (xj) (287) 
from M" to ] — f j f ['^ and consider the transformed function 



w- 



4''\T,y):=vl'P\T,x) (288) 
for 1 < i < n. We have 

Qyf^P^'- dw{'^''' dym dwl''''^ 1 dw{'^''' 1 



dxm dym dxm dym l + a;2, dym l + tan^{ym)' 

and 



(289) 



(290) 



dy'L Vl+tan^{j/m); ^ dym (l+tan2(y™^))'* 



Since 



n2 /iP:^ 

(l + lxl^) 10 as b|too (291) 

^ ' oxi 



m 



Hence we get a series of functions (w''''^''')^ defined on transformed domain 
Dl''" := [I -l,l]x ]-7r/2,7r/2[" along with 

w'''''^(r,y)=v''''''=(r,x) (292) 

and the limit is in the Banach space [Ci^a/2,2+a {-^I'^)]^ with the norm 
|.|" restricted to the domain Dj'^ . Note that we have tt;f'^''(T, s) = for 
all 1 < i < n if for some sj in s = (si, • • • s„) we have sj G {— 7r/2,7r/2} 
according to the polynomial decay of v''''^'' at infinity. Hence we get a 
limit w'''' with w''''(r,y) = v^''(r,x) for all (t,x) G DJ and (r, y) e Dj'"" 
where satisfies the transformed Navier-Stokes equation (in r, y-coordinates 
on D;"'"). 



2.2 step 2: Extension of local existence to extended equa- 
tions for v'' '''' 

In step 1 we did not consider the growth of the solution. We just proved 
that the iterations i^™(v''''~^), Z > 1 of the functional Fi are in a function 
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set 5*0 such that first order terms of the linear subproblems are uniformly 
bounded. In order to get a globally bounded scheme we introduce a control 
function r = (r^)/ and set up a globally bounded scheme for v''''' = v^ + r. In 
this second step we observe that the extension of the local iteration of step 
1 to equations with bounded regular functions r 7^ does not change the 
argument of step 1 essentially. Note that the Leray projection from equation 
for w'^'P includes integral terms for r\ as described in the introduction. This 
means that we have to estimate magnitudes of the from 



Again this reduces to estimates. The similarity of the convergence ar- 
gument of the local scheme for ^^'^'P'^ to the argument for convergence of the 
local scheme of step 1 is due to the fact that the additional terms symbol- 
ized by the term L^'' (r, v) in the equation for w''''^'^ consist of a source term 
dependent on the function r (which we called ^'(r)), and a linear operator 
dependent on the function v''''''' which we called L^'''" (r, v) (resp. v^'P'^''- 
for each step of the local iteration). Note that both additional terms have a 
multiplier pi with the exception of the time derivative of ^ . If is regular 
for given time step number i.e., if the first time derivative and the sec- 
ond time derivatives exist and are integrable on the domain {1 — 1,1] x M", 
then we get a regular solution v'''''''''' of the linearized Navier-Stokes equa- 
tion at the first step of the local time iteration. Furthermore, the additional 
source term 5'(r) depending on the function r appears only in the equa- 
tion for v'''''''^'', because these source terms cancel out in the subtraction 
^^r,p,k,i _ ^r,p,k,i _ ^r,p,k-i,i (^gfj^ing the higher order corrections of the lo- 
cal iterative scheme. It is sufficient to find a lower bound pi ^ j for the 
number pi in order to make the scheme global. 

Remark 2.17. For numerical purposes we shall consider a different solution 
scheme below. If we know by analytical means that v'' and r are bounded 
then we know that v is bounded, and this implies that we can set up a 
scheme where we do not use the function r at all. A second feature of a 
numerical scheme is that diffusion may imply that increasing time step sizes 
Pl may be chosen as the time step number / increases. This makes the scheme 
more efficient. 

Next we construct 



In this section our goal is make the extension of the argument in step 1 to the 
case where r 7^ as easy as possible. Therefore we write down an iterative 
scheme for (|294p which we would not choose for numerical purposes. More 




(293) 



00 




(294) 
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precisely, for each A; > in the computation of the series ()294p we have 
the additional term L^'''' with arguments of the previous time step. This 
means that in step A; = we add a term L^''''' (r, v''''"^) on the right side 
in order to compute the approximation tjp^'"'' (note that the argument r 
is considered to be externally given in this second step of the proof; we 
shall determine the function r in the third step if this proof). Similarly, in 
step A; > we add a term L[''''' (r, v'''*^"^'') on the right side in order to 
compute the approximation y'^'f'^'^ The first term of this functional series 
^r,p,o,i _ ^^r,p,o,i ^ _ _ _ ^^r,p,o,z^ gQ^^j^^Q]^ |^]^g equation (for 1 < i < n) 

+pi /r" {£-Kn{x - y)) Elk=i ( ^""al, ) i^^y)dy + <r> 

^^'''P'^^\l - 1, .) = y'-'P'^-^{1 - 1, .), 

(295) 

and the other terms Sv^'^'^'^ of this functional series are the respective solu- 
tions of the equations (for 1 < i < n) 

dSv'.'''-'^'''' d^Sv'^'''''''' r,p,k-l,l dSv^'''''''^ t fulfil r o k- 

, 1 I n r,p,h — l,l 

-O, V ^ r,p,k-l,ldvJ^ • 

PI /m« [I-Kn{x - y)) E • _i k;;— J (r, y)dy 

Q^^Knix - y)j Elk=i [^-^, ^S^J (^^y)dy, 

S^r;p,k,l^l _!,.)= 0. 

(296) 

Here in case k = 1 the expression ig defined by 

r,p,—l,l r,p,l—l /nr\^\ 

Vf' := Vf' . (297) 

We assume that from the previous time step Z — 1 we have the estimates 

\r'r\ < C°, < C), \r\-\,2 < C,, (298) 

and that 



■'lo < \v\r\,x < C[-\ l^!"'li,2 < CS- (299) 



4" 10 ^ '--O ' l ^i 10,1 ^ '--I . 

1=1 
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Furthermore we assume that 



r,p,«-l 



2 ; ' 



and that 



J-i 



< C: (Cr + (/ - l)Cr) 



(300) 
(301) 



for ah 1 < A; < n and generic C*. Note that this implies 



j,k,m,,p=l 



dv 



r,p,l—l 



dxj 



dv 



r,p,l~l 



dXr. 



2 ; ' 



and that 



j,k,m,p=l 



I dr 



i-i 



dr: 



dxr 



iT,y)dy<c:{Cr + {l-l)Cr 



(302) 



(303) 



If we estabhsh the latter estimates for all / then we have a linear growth of 
the integral terms of Leray projection in the equations for v^'f'K 

Remark 2.18. In the next subsection we shall define and determine the 
constant Cr- We shall show that there exists a constant Cr > (independent 
of I) such that for all 1 < i < n the functions r- satisfy 



k'|l,2 < Cr 



(304) 



uniformly, i.e., independent of the time step number /. Here the norm is 
defined on the local domain Df = (/ — !,/] x R*^, i.e. we have 



k-|i,2 := V sup \r[^{T,x)\+ sup \r[^{T,x) 



(305) 



for multiindices a related to the spatial variables (especially, we have |r' ^| < 
Cr for all 1 < i < ?i). 

In step 3 of this proof below we construct such that the relation ()304p 
and the relation ()303p are satisfied for all / > 1. Next recall that 



-Pi /r" (£-^nix - y)) Elk=l 
where 



Lf (r; v) = +pi J2j=i Tj^^ + Pi E,=i ^ 



dxj 



(T,y)dy + Lf''°(r;v'-.^.'-i) 



n r,p. 
j=l '"j 



j 



-'^Pl /r" [wi^nix - y) ) Elk=i 



(306) 



(307) 
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We may use the rough estimate 

\LP'\r; v'-'^.'-i)| < piuCr + pmCf. + PiCkC^. + |Lf'''°(r; v-'^''-i)|. (308) 

which we get from (j304p . This is not needed essentially for the proof of local 
convergence of this step. For this purpose we use ()304p again, and we get 
from (j307p the estimate (generic use of C*) 

|if'°(r;v-'''''-i)| <pzQ|r%,i|v^'''''-i|o,i (309) 

For I > 1 given the functional series (f['''''^'')fc is determined by the 
representations 

?;[''''°''(r, j;) = u['''''^-^(/ - l,y)T[){T,x;l - l,y)dy+ 

/iRn Pi J^n {£-Kn{y - ^)) (Elk=i ( ^"al, ) ^)) X 

xrQ'(r, X, s, y)dzdyds 

+ /(l-i) /r" ^f''(''' v'"''~^)(s, y)T\f{T, X, s, y)dzdyds 
+ /(l-i) /r" 2/)ro'(r, X, s, y)dzdyds 

(310) 

and (for A; > 1) 

(5u['^'''''(r,x) = /(l.i) /ijn (Ei=i ^"'''a^,""''' ) (s,y)r^''(T,x;s,y)a!sdy+ 

/{«-!) Pi /r" /r" ^n,i{z - y) X 

(Ej,m=i {'"^j',m~^'\s,y) - v2^''~^'\s,y)^^ ^T'f:\T,x;s,z)dydzds, 
+ '° (r, 5v^.A'=-i'') y)dzdyds 

(311) 

where the functions F^' , A: > are fundamental solutions of 

dr'?' v^n d'^r^'' v^n r,p,k—l,l dV^!^ „ /oioN 

- Pii^ J2j=i + Pi J2j=i = (312) 

(recall that v^'''' = -y^'''''"^ in case A; = 0). Analogously as in the first 
step of this proof for some time step / let v''*''''"^ with v^'P'1~^ for some let 
Fr^i : Cl''^{Dl) such that v'^'^''^-' = F/;,(v'^'^''-1). We assume that 
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and start the iteration with 

^r,p,l-l ^ gr ._ jfllfin^ < 4C[-;^^^ . (314) 

Hence, similar as in the first step we get for some generic constant C* 
dependent on the dimension n 



|5^r,p,fc,Z|n^ < piC*Cr + Cr) Cp | -JV-'^-^'' |? 2 

(315) 



for all A; > if we choose 



Pi < ^— r (316) 



In order to have similar estimates with respect to the |.|j|^2-norm we use 
constants pi which decrease linearly with respect to the time step number I, 
i.e. we use 

^' " 4C*Cr ( [Ci;^' + Cr)+l (^{72^ + Cr) ) Cr ^^^^^ 

Moreover with analogous arguments as in the first step for all r € [Z — 1, /] 
we have 

\5v^-'P''^\t,.)\h2 < \\5v'^''^'-'Hr,.)\H^. (318) 

uniformly in r. Then we use either method a) or b) completely analogously 
as in the first step of this proof. Hence locally and for pi as above we have 

/ 12 

got a function v''' G ' which is constructed from a limit of the functional 
series 

oo 

r,p,l r,p,0,l I \ ^ r r,p,k,l /oin\ 

=Vi'^^ + l^Svi'^' (319) 

k=l 

where v*"'^'' is a classical solution of 



-w- - Pi'' z2j=i + Pi Ej=i 



(320) 



^ v'-'P'^{l-l,.) = v'^'^''-i(/-l,.). 
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Finally, we note that by adjusting C* in our choice of pi we can ensure that 

oo ^ 

J^|5v^.P.M|^2<-. (321) 

k=l 

2.3 step 3: Construction of the function r and of a globally 
convergent iterative scheme for v''^ ' 

We have proved the existence of local solutions for v^''''' = v'"'' + r' for 

II 12 

some fixed class of functions r and v'''''' with components in ' and such 
that the components are in the Hilbert space or in some Banach spaces 
H'^'P in order to have convergence with some suitable regularity of the limit 
function. A 'suitable regularity' of the limit function is Holder continuity 
which is uniform in time. This implies another simple reason for choosing 
the classical space: it allows us to have classical representations in terms of 
fundamental solutions not only for some approximating linear equations but 
also for the limit function once it is known. The classical representation in 
terms of fundamental solutions then allows us to get more regularity. Note 
that we used regularity of coefficients when we applied the adjoint equation 
in the first and second step of this proof. Furthermore the existence of 
classical solutions is useful because some forms of the maximum principle 
require the existence of classical solutions. These observations are useful for 
our construction of the function r which controls the global growth of the 
function v''"'' and is itself bounded (and has an integral magnitude which has 
linear growth with respect to time). The time step size at time step number 
I (measured in original time coordinates t) is pi and we concluded that the 
local scheme converges with a time step size of order 



c*n ( {c[;,' + Cr)+i [ci;,' + a) ) ^c} ' 

where C|~^^ is an upper bound for v^'''''"^, where Cr a certain upper bound 
bound for r' with respect to the |.|i^2 norm. Indeed we determine the con- 
stant Cr in terms of the initial data function hi in this section. Furthermore, 
at the end of this section we shall show that this choice of Cr is also an uni- 
form upper bound for all C[ 2 if Pi is chosen as in (j322p . Concerning this 
constant C* we note that it depends only on dimension and and may be 
computed explicitly. The constant Cr in (j322p is a uniform bound of con- 
stants related to the fundamental solutions which appear at the kth substep 
of time step I (cf. step 1 of this proof). Step 1 and step 2 of this proof show: 
if the constants C[~2 and Cr are given then Cr can be computed and local 
convergence is given along with the choice of pi as in (|322p . (Note that we 
have shown in step 1 and step 2 of this proof that the other constants in 
()322p depend only on the dimension and the initial data h and, hence, are 
independent of the time step number I.) 
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We have to show that the numbers C{ 2 ) I > can be estimated in terms 
of a constant Ci^2 which is independent of / in order to get a global scheme. 
We choose a time step size of order ~ j in order to deal with the growth 
of the integral magnitudes related to the integral term of the Navier-Stokes 
equation in Leray projection form. Furthermore we have to show that the 
function r can be chosen together with a finite constant Cr independently 
of the time step number I. In this step 3 of the proof we determine the 
constant Ci^2 and the constant Cr which is a bound for the global function 
r. In order to reduce the number of constants we shall determine Ci^2 such 
that 

|r^|i,2 <Cr = Ci,2 (323) 

for all / > 1 and all 1 < i < n. Since we have a step size with respect to 
original time coordinates of order ~ j we ensure that for the growth of 
the integral magnitude we have 




dri dr\ \ 



(324) 



i.e., we have linear growth of the integral magnitude with respect to the 
time step number I. This ensures that related integral magnitudes of the 
equation for v^''''' in Leray projection form have a linear bound, too. Note 
that the first derivative of the kernel K is bounded especially in case n = 3. 
In order to estimate the quantity ()324p we may estimate 



E 



j,k,m,p=l 



I dx^ 



dXr, 



ir,y)dy<C:iCi,2 + lCi,2) 



(325) 



inductively. However, it is useful to reduce the estimate (j327|) to an 
estimate of the form 



E 

j,k,m,p=l 



dxi 



+ 



dr 



dXr. 



iT,y)dy<Cl{Ci,2 + lCi,2)- (326) 



Note that our local schemes of step 1 and step 2 of this proof shows that we 
need also estimates of derivatives of the integral terms, i.e., estimates of the 
form 



E 



d'^r[. 



j,k,m,p,q=l 



I dXidXn 



dXr 



[T,y)dy<C:iCi,2 + lCi, 



(327) 



However, this can be treated similarly, i.e., all this reduces to an estimate 
in terms of a norm, i.e., to an estimate of form 



|ri(r,.)|^,2 <Q(Ci,2 + /Ci,2) 



(328) 
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where we note that C* is a generic constant. 

Depending on the representation of the function r- we may iterate the 
reduction of mixed products to sums of squares. Therefore we shall establish 
the estimate 

'\T,y)dy<Ci,2 + lCi,2, (329) 



and the estimate 



/ 






dxj 


/ 






dxj 



iT,y)dy<Ci,2 + lCi, 



(330) 



inductively when we have constructed the function r'. This involves related 
estimates for functions up''''"^ and which occur in the construction of r-. 
We shall show that r can be constructed via certain source term functions 
4>{ of linear parabolic equations such that ()323p holds, and such that for all 
^ > 1 we have 

Cr, (331) 



\ r,p,l\n < ^ 
|v \i2 2i ^-1,2 



where Cr depends only on the dimension n and the initial data h of the 
Navier-Stokes Cauchy problem. We shall also show that 



v^'^'-'l^^ <C:(C7i,2 + ZCi,2) 



and 



r''%2<C:iCi,2 + lCi,2), 



(332) 



(333) 



and that similar estimates are available with respect to other i^^'^-Banach 
spaces. This means that our method can be applied to general dimension 
n. This is outlined in the introduction of this section, and we shall fill in 
the details now and in the next step of this proof. In order to have an the 
equivalent problem (i.e., a problem equivalent to the original Navier-Stokes 
Cauchy problem for v) for v^'P''- we know that r has to be bounded (well, a 
bound linear with respect to the time step number / would be sufficient as 
well, but we show that we can construct a globally bounded function r such 
that the problem for v^'^ is equivalent to the problem for v). The problem 
for v''''' has the advantage that we can solve it step by step where we control 
the growth by the functions r'. Note that ()33ip implies that for the original 
equation we have the bound 



|V^''|5^2 < Ci,2+nCr 



(n + l)Ci,2. 



(334) 



The construction of the function r and of the global scheme for v'^'f is done 
inductively with respect to the time step number I — 1. At each time step 
I > 1 assume that the functions v^'^''~^(Z — 1, .) and r''^^{l — 1, .) have been 
computed. In case / = 1 we set v^''f'''~^{l — 1,.) = h and r*^ = 0. We use 
these data and proceed at each time step I in three substeps as we outline 
next. 
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3i) At each time step Z > 1 we construct first the functions which 
serve as 'consumption source terms' and are constructed such that 
the growth of the functions v''"'' and is controlled on the domain 
[/ — X M". Control of growth is both in the sense of supremum 
norms such as |.|o and |.|i^2 and in the sense of integral norms such as 
|.|x,2 and |.|//2. The functions (f)\, 1 < i < n are constructed as a sum 
of two functions, i.e., 

= + '/'I^ (335) 
where (j)^'^ is mainly constructed in order to control the growth of v^'P'^ 
and is mainly constructed in order to control the growth of r' . Each 
of these functions is defined in terms of the data v'''^''~^(/ — 1, .) and 
r'~^(/ — 1,.) respectively. Concerning the growth in the sense of the 
supremum norm |.|o in the regions of the domain [/ — 1, Z] x where 
the value of ^{l — 1,.) is above a certain threshold y we shall 

define 0^'' = — 1 and in the regions of the domain [/ — 1, Z] x M" where 
the value of vl'^'^~^{l — 1, .) is below a certain threshold — ^ we shall 
define (^^'' = 1. There are two possibilities here: for given 1 <i < n the 
modulus of the function r\ may be itself below or above the threshold 
C/2. We shall define a scheme such that in either case the growth of 
the function f is controlled. Moreover, locally and for each i, if 
both function values r'~^(r, x) and v\'^'^~^{t,x) have opposite sign at 
some point (r, x), then the construction of r\ and (f)\ is such that the 
growth of all functions r\ , t^p^'' , and u^'' is essentially controlled by the 
behavior of the function r\ if the latter function exceeds some level. 
In this way in special situations we control the growth of the function 
r\ and of the original Navier-Stokes solution v^'^ and this leads to the 
growth control of the functions on [Z — 1, Z] x R" (in r-coordinates) 
is controlled. Similarly, in the regions of the domain [/ — 1, /] xM" where 
the value of the function f is below a certain level, and the value of 
the function r^j^^{l — 1, ■) is above a certain threshold y we shall define 
= —1. Furthermore in the regions of the domain [/ — 1,/] x M" 
where the value of the function fp'''' is below a certain level, and the 
value of the function — 1,.) is below a certain threshold — y 

we shall define (j)'/ = 1. The functions are defined via partitions of 
unity in order to have regular functions. The sum (f)\ will appear on 
the right side of a linearized equation for r' first. Furthermore, for 
each 1 < i < n we shall ensure that in the complementary regions 
where the modulus of both functions T^j~^{l — 1, .) or v^''^'^~^{l — 1, .) 
is below C/2 the definition of the functions is extended such that 
the integral magnitudes of the functions have a linear upper bound. 
Since products of functions can be pointwise estimated by sums of 
squares it is sufficient to have linear growth of if ^-estimates. The 
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estimation of the solution expression for involves (first derivatives 
of) fundamental solutions of approximating linear equations integrated 
over time (one time step) and space. These integrals are absolutely 
bounded. The functions are defined in the second substep which we 
sketch in advance next. 

Once we have defined the functions 0' we define the function r' on the 
domain [Z — 1, /] x M" via linearized equations of the form 

+Pi /r" {£-^n{x - y)) Elk=i (%- %r ) " ^' 

En l-l _ sr^n r,p,l-ldr^2_ 

j=l '^j ^j=l dxj 

< 

+^pi /m" {£-Knix - y)) E • .=1 [^'-^) (1 - ^,y)dy 

-PI /m" {I-Knix - y)) E • ^=1 ( "l^r'ilT' ) - 1' y)dy + 4, 
, r'(/-l,.)=r'-i(^-l,-)- 

(336) 

Note that the functions on the right side with time index I — 1 are 
defined on the domain [/ — 2,/ — 1] x M". Hence these functions are 
evaluated at r = / — 1. Note that the function (^,- is defined on the 
domain [Z — 1, /] x W\ However, we shall define (/)[ as functions which 
depend only on the spatial variables, i.e., it they are constant with 
respect to time. We outlined the reasons for this in the introduction. 
Moreover, we assume that r'~^(/ — 1, .) has been defined at the previous 
time step. At time step / = 1 we have defined = Note that all 
functions on the right side of the first equation in ()336|) have a factor 
pi except for the function (/)[ . Choosing pi small enough we ensure that 
the maximum value of each component of the function r-, 1 < i < n 
of the function r' is bounded in terms of the maximum of r-(/ — 1, .) = 
r'~^(/ — 1, .), 1 < i < n respectively for each i. Furthermore we shall 
see that we have a bound to the integral magnitude ()324p where it is 
essential at time step I that (/)^ is defined in substep 1 such that the 
l^p norm of 

/ r^k^i} -^^y)^iiT^x]s,y)dy + I / (j)i{s,y)Tl.{T,x; s,y)dy, 
Jr" Jo jm" 

(337) 
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is finite. Similarly for derivatives. The definition of cj)\ in terms of the 
functions r'~^ and v\''''^~^ makes this possible. 

3iii) Next we plug v^''''' = v'''' + r' into the Navier-Stokes equation on 
[/ - 1,Z] X M". Then we get 



dr 



En d V- ■ ' s-^n r 



n r,p,l (9f['^'' _ ,/ 



j dxj 



, ^r''^P^\l - 1, .) =y'-^p^^-\1 -1,.), 



(338) 



where 



-Pl 



r,p,l dr^ 



(339) 



dvl'"'' dv, 
dxj dxf^ 



{T,y)dy. 



In order to control the growth of v'^''''' we consider the functional series 



rr,P,l 



^r,p,o,i _j_ ^s^i'^p^k,i ^-^^ consider a) the growth of v' 



r,p,0,l 



and b) the fact the sum of correction terms X^^Li (5v^''''*^'' is small for 
small Pl- The equation for 



En d v. - v-^n r,p,L—L/, 



n r,p,l—l I 



, dv. 



r,p,0,l 



dxj 



^ ^r,pfi,l(j^ - 1, .) = v'''/'''-l(Z - 1, .), 



(340) 
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where 



-Pi k« {m-Kn{x - y)) Elk=i y)dy 

I v^n. ; ati^'^''^"*" , ^-^n r,p,l-l dr\ 
+Pl l^j=l '^j dxj + Pi l^j=l '"j dxj 

-'^Pi /r" {£-^n{x - y)) Elk=i (Sj ^^fe ) {r,y)dy 
+P«/r" {m-Kn{x - y)) Zlk=i ( ^"al, ) - 

(341) 

Then we analyze the growth of the functions v^'^' ' , 1 < i < n, where 
we estimate the "01'" using p36p . and then estimate the difference 
V'i ~ i^i^- Note that the estimate for the growth of v^''^'' involves and 
estimate for the Navier-Stokes solution v'''' and the control function 
r' if the modulus of components of the latter function exceed a certain 
level in critical regions where the modulus functions wp''''"^ exceed a 
certain critical level. 

ad 3i) As we described in the introduction of this section the function 
(pl is determined as a sum of two functions 4>l' and 4>^' . Let us consider 
the first induction step / = 1 first. At time step number / = 1 we do 
not have a function r^~^ (all we have are the initial data h). We may set 
r*^ = = 0. Hence there is no need to control the growth of on the 
domain [0, 1] x an we may set (p^^^^i = 0. This means that we do not 
involve rl in order to control the growth of the Navier-Stokes solution u^'^ 
(in r-coordinates) at the first time step / = 1, i.e. on the domain [0, 1] x M". 
Hence, is defined solely in terms of a function (f)^'^'^ : — )• W^. Next 
we define i;^^'*'^ to be a function which satisfies the following lemma with 
B = D% . = \x\vI'P'\q,x) = hi{x) G [f and 

A = Dt} = {x\vl'''\0,x) = h,{x) G [-C,-f ]}. 

Lemma 2.19. Let A,B CI be two closed disjoint sets. There is a function 
/ : M" ^ M such that 

1) / G C72(]R"), where 

l/lo,2 < C+, (342) 
where is a constant which depends only on the dimension n and 
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the distance of the sets A and B. (recall that we define 

|/|o,2 = sup V \d^f{x)\ (343) 

-eK" |a|<2 

along with multivariate derivatives da = with multiindices a = 
(ai, • • • , On) and \a\ = ai, and where > for all 1 < i < n.) 

2) we have sup^^j^n < 1 by construction, and 

1 ifxGA, 

fix) = { (344) 
-1 ifxeB. 

The lemma (|2.19p can be proved using classical partitions of unity. 
Proof. Let 

where dist(^, S) := sup^^y^ y^^ \x — y\2 denotes the Euclidean distance. For 
integer tuples p = (pi, • • • ,p„) G Z" define a differentiable partition of unity 

</>p^(x) = ntiV'f--Pil (346) 



where -i/; : R — t- M is defined by 



^<y>=m ''''' 



along with 



Hy) ■■= { 



(348) 



and 



exp (-i^p-) ,if \y\ < 1> 
else, 

if(y) = ^/i(y-m). (349) 

Note that 4>p^ has as its support cubes of side length /i of form 

supp((/)p^) = {x £ M'^llxj — PiiJ,\ < /i, 1 < z < n} . (350) 
Furthermore, note that ((/)p^)pg2n is indeed a partition of unity, i.e. 

4)p^{x) = 1 for aU X G M". (351) 
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Next define the function / : 
/(^) = 



M" — 7- M of the lemma via 



where 
and 



Pa := {p G Z"|supp(</.p^) n .4 ^ 0} 



(352) 



(353) 



:= {p G Z"|supp(</.p^) n P / 0} . (354) 

Computing the derivatives of / involves only finitely many terms (atmost 
2"- either of the first sum indexed by Pa or of the second sum indexed by 
Pb- Especially we get 

|/|o,2<C+:=2"C^, (355) 

where = max2;ggupp(0p^) S|a|<2 ^a</'p/x(a^)- Note that the latter constant 
depends only on /x and n. □ 

Note that a fortiori the first part of the preceding lemma especially im- 
plies that 

l/lo = sup \f{x)\ < C+dist {A, B) . (356) 

This implies that the distance of the sets A and B has a lower bound in 
terms of the norm |/|o,i, i-e. 



dist(Ai?)>i|/|o, 



(357) 



where > is some constant depending only on dimension n. Such a 
lower bound can be established independently of the time step I. It is clear 
that this fact is closely related to the boundedness of the functions v'^''' (and, 
hence, v''). It is also related to the linear bound of the integral magnitude 
in p2|) and similar bounds for the functions v''''' and v''. Clearly, we want to 
have a lower bound which is independent of the time step number /. We have 
now defined 1 < i < n and = 0. The growth of the latter function 
is controlled by definition. Next we define the functions for I > 2. We 
assume that the functions r'"^(Z — 1, .) and v'''''''~^(Z — 1, .) are given. For 
each 1 < i < n the functions are constructed as a sum = (f>^'^ + 
where both functions are determined in terms of the data v^'f'''~^{l — 1,.) 
and r'~^(Z — 1, .). These functions are obtained at the previous time step 



I — 1. In order to construct the function 



we define 



D 



vA-l 



r,p,l- 



and 



D 



XV, 



r,p,l—l 



{I -l,x)G 



^ c 

2'^ 



c 9_ 



(358) 



(359) 
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Since "^(/ — 1, .)|o,2 is bounded there is a distance 



> 6 



(360) 



for some constant 5 > 0. We shall show that this distance can be preserved 
at time step /. We apply lemma [2.191 with A = D^^~^ and B = D^''~^ for 

each 1 < i < n and get the functions 0^'''^'' : M" — )■ M. We define 

1 V,S,l I 



C':(/-l,Z]x 



I, 4>l {t,x) := (x), 



(361) 



and according to lemma [2. 191 we have c^^''^'' G C't'^, where 



4>i [x) 



' -1 ifxez5i;ir\ 



a;- fx) else 



1 ifx€ D, 



(362) 



for some function a G Cl'"^ with a(x) G [0, 1] for all x G IR"\ ^-Dj'^ ^ U Dj'l ^ 
This function a(a;) is not arbitrary. Especially it is useful to ensure that 



a£H^ on its domain M*^ \ [^D^;+^ U D";!"' j • Since 
as a source term of a linear equation for r • we define 



v,l I jr,/ 



appears 



r,p,l—l 



{l-l,x). 



(363) 



where x G M" \ [of-^UDfj^ 



" \ Dj"'' ^. On this domain this 
function is clearly in (by induction). Furthermore the function is 
globally Lipschitz, hence Holder continuous. The functions are defined 
similarly. Recall that we defined subsets of D^^~^ and D^2^~^ of 

'C 



D^2^f ^ which are mutually disjoint with the sets Dj'^, ^, Z)}''' ^ . We defined 



D'_;[-':=\x\ri{l-l,x)e 



-,c 



(364) 



and 



D 



rA-l 



:={x\rl{l-l,x)G }, 



(365) 



For our analysis it is important to consider certain subsets of the sets D 
and D^'^~^: the set 



D 



r,0,i-l 



x|r-(/ - 1, x) G 



2'^ 



r,p,l-l 



{l-l,x)<^D\ 



vl-l 



(366) 
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and the set 

(367) 

where we may define (f)l''' rather independently of cl)J'K On the complemen- 
tary sets 



and 

we ensure that we can control the growth of the functions v^'''''' and r • simul- 
taneously. How to do this is considered in detail in the third substep of this 
third substep. Again, since \rl~^{l — 1, .)|o,2 is bounded there is a distance 

dist (^D!;[;\ D^'Ji^^ >S (370) 

for some constant 6 > 0. Note that this constant 6 is not a small constant 
but bounded from below uniformly as a consequence of the upper bound 
with respect to the |.|o.2 norm of the functions r^~^ and up''''""^. We shall 
show that this distance can be preserved at time step I in the third substep 
below. Note that this implies that the constant 6 in (j370p can be chosen 
independently of the time step number I (cf. substep (iii) below) . We apply 
lemma [2.191 with A = D^^~^ and B = D^^~^ for each 1 < i < n and get the 

functions : M" ^ M. We define 

f/ : (/ - 1, Z] X M" ^ M, f/{T, x) := <P-'''\x), (371) 
and according to lemma [2. 191 we have (/>['^'' G C^'^j where 

' -1 ifxGZ)[;|-\ 



(3{x) else (372) 



1 ifrrGA^'i"', 



for some function /3 G Cl'"^ with /3(x) G [0, 1] for all x G ]R"\ U L>[;i~^ 

Again the function /3 is not arbitrary. We define 



f3l{x) = -^rl-Hl-l,x). (373) 

on the domain M" \ ^-D['^^ U -D,-'i~^^ . Now we have constructed the func- 
tions = (f)1'^ + This finishes substep i) of step 3, i.e., the defintion of 
the functions (b\. 
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ad 3 ii): Next we construct the functions r' for 1 < i < n. We want to 
define constants C^,C^,Cr which are independent of the time step number 
I and such that for ah / > 1 

Ir^lo < C^, \rl\o,i < C}, |r^|i,2 < (374) 

holds, where the norms are defined as above with respect to the domain 
DJ = (/ — 1,/] X M". Furthermore we want to ensure hnear growth of the 
intergal magnitude for the function (see below). At the first time step 
I = 1 we have defined = 0. For I > 2 the functions r- are defined by the 
Cauchy problem ()336p . Note that the right side is evaluated at r = Z — 1. 
We rewrite (|336|) in the form 

< Si{l-l,.) + <Pi{l-l,.), (375) 

^ r'(/-l,.) = r'-i(/-l,.)- 

where 

j=l fj dxj Pi 2^j=l dxj 

+^pi {irKnix - y)) E • .=1 [^'-^) ii - ^,y)dy 

-PI {£-Knix - y)) Elm=i [^^'-^) (1 - hy)dy. 

(376) 

This is a linear parabolic equation for which can be solved. Furthermore, 
the solution has a representation in terms of the fundamental solution Pj, of 
the equation 

rl, - piuArl + PI - 1' = 0, (377) 

which allows us to analyze the growth of r'. We have 

= 4„ rl-\l - l,y)Tl{T,x;0,y)dy 

+ Ii-i /r" - 1' y)ri(T, x; s, y)dsdy (378) 
+ /i-i /r" 4>\[s, y)r{.(r, X] s, y)dsdy, 
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We shall exploit the following ideas related to the representation (|378p : con- 
cerning the the relation of the second and the third term on the right side 
of ()378p we note that all terms of the function S'' have the factor pi while 
the third term (p^ does not have this factor. If the factor p; > is small, 
then the diffusive effect of the fundamental solution is small in the sense 
that the values of the second and third term on the right side of ()378p de- 
pend largely on the source terms Sl{l — l,y) and (/)[{s,y) (we shall make 
this remark precise in a moment). This means that for pi small enough for 
each 1 < i < n the term in (|378p involving dominates the term involving 
S"'. Furthermore a bound for the first term of the right side in ()378p can be 
obtained using the maximum principle. Considering derivatives and shifting 
of derivatives using the adjoint of the fundamental solution as we have done 
in the first step of this proof leads to an estimate of the functions in the 
2 norm. Furthermore, we have the estimates 



\rl{r,.)\m < |E"=i/Mn fRr.rl-\l - l,y) (\r[{T,x;0,y)dy\ + |r^^^.(r, 0, 
+ /«-! /k" 'Pii^'V) (\^UTiX'^s,y)dy\ + |r' (r, s, y)dy| ) dsdy 



< C* 



and 



J-i 



+ C* 



r,p,l—l 



(379) 



\rl{T,.)\H2 < Ej,m=i/iR" " (l^^r* (t, x;0, y)dy\ + \T^;*{T,x;0,y)dy\ 

+ /«-! /e" (l^i,rnis,y) {\^r*iT,x;s,y)dy\ + \r^;*j{T,x;s,y)dy\] dsdy 



< C* 



+ c* 



r,p,l—l 



(380) 

with our generic constant C* > independent of the time step number I. In 
this step we use the definition of (p\ and the adjoint as in step 1. The latter 
fact then ensures that we have a linear growth of the intergal magnitude for 
the functions r-, i.e. the magnitude 



dy. 



(381) 



is bounded by Ci^2 + '^^1,2 for r G [/ — 1,Z]. Indeed we shall establish this 
linear bound for an upper bound of ()38ip . i.e., we shall find that 



n 

j,k,m,p=l 



K 

dxi 



dy<Ci,2 + lCi,2. 



(382) 
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Since we estimate the functions (r, .) and r'(r, .) with respect to the 
norm |.|o,2 (this means spatial derivatives up to second order) we construct 
similar bounds for the first derivatives of ()382p . Again, since products of 
functions can be pointwise estimated by sums of squares estimates the func- 
tions v[''^'\t, .) and r,-(r, .) with respect to the norm are essential. Note 
that the choice of ~ j compensates for this linear growth in order to esure 
convergence of the local scheme with respect to intergal norms of step 2 of 
this proof. Next we have to determine the upper bounds Cr and Ci^2 of the 
functions r and v'''''. We define the constants Ci^2 = Cr just in terms of the 
initial data function h (dependence on the viscosity umay be put into the 
step size). This determines pi, and we check that r' is bounded indepen- 
dently of the time step number I and that the integral magnitude of r' is 
bounded, too. We define 



a = ci,2 = 2 + 2|h|?,2 + /m. E • M,m=i wt(y) 



dy 



(383) 



+ /r" Y^j,k,i,m,p=i dxjdxp (y) dxL 



dh. 



dy. 



The last term in the definition (j383p is due to the fact that we do the 
estimates for local |.|i,2-iiorms, i.e., on the domains {1 — 1,1] x M" and global 
norms for the functions f (r, .) and r-(r, .). Next we prove that the 
constants C^, and Cr are bounds for the functions r- with respect to the 
norms |.|o, Mo,i) and |.|i^2 which are preserved at the Ith time step. Assume 
inductively that 



J-i 



< \r--\,i < a\ \r]-\,2 < Cr. 



(384) 



We have not defined and yet but it turns out that the constant Cr is 
essential. Hence we just define 



Cr '. — Cr , and Cr '. — Cr . 



(385) 



Clearly a closer analysis leads to some bounds , < Cr , but this is 
not essential. However we define r', equation for r' includes an the integral 
term of the form which we know from the Navier-Stokes equation in its 
Leray projection form. At time step / this term may be given in terms of 



the function ^ of the previous time step. Hence, it may seem natural to 
assume inductively that 

n 

, ^ I dxi 

j,k,m,p=l ■' 

Preservation of these bounds at the Ith time step means that 
Ir^lo < k^lo.i < C}, |r^|i,2 < Cr, 



-^{l - 1, y) dy<Cr + {l- l)Cr (386) 



(387) 
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and that 



j,k,m,p 



■{■Ti y) 



dy <Cr + ICr 



(388) 



for all r G [/ — 1,/]. Since we have not estimated v^'''''"^ yet we can prove 
this only in terms of estimates of the latter function which are inductively 
assumed. Our local method described in the first and second substep of 
this proof indicates that we should require a similar requirement for the firs 
spatial derivative, i.e., that we should require inductively that 



E 

j,k,rn,p,q=l 



dXjdXq 



dr: 



dxr, 



-{I - dy<Cr + {l- l)Cr (389) 



Nite that by definition the functions r\ depend on the source terms (f)\ of 
their defining equations, and these source terms depend on f['''''~^(/ — 1, .) 
and r'^^ in turn. Hence, we assume inductively that 

\Vi lo S , |Uj |o,l < , |o,2 < <-'o,2 > ll,2 S <^i^2 

(390) 



for some finite constants Cq ^ , C| ^ , C| 2^ . In the end we shall have 



(391) 



and we may set Cq~^ := Ci^2, and C[~^ = C|^^, and Cq~2 ■= C[~2i but we 
shall show this later in substep (iii) below. Hence we keep the index I for 
this time. Furthermore, we assume inductively that 



j,k,m,p=l 



dxj 



(i-hy) 



UUm 



il-l,y)dy<C['-2' + {l-l)C[-2 



dXr, 



and substep (iii) below shows that this linear growth and that of the 



(392) 



/ 



E 

j,k,m,p,q= 



I -il-l,y) 



^ I dXjdXq 



(/-l,y) dy<C[-2^ + {l-l)C[-2\ 



dXr, 



(393) 

is preserved as well. Well all the inductive assumptions on integral terms 
d vl''^'^^^ can be expressed by ont 



involving r\ ^ and ^ can be expressed by one condition which is 



\H2 



< c: ( c'r^' + (/ - i)c[~2' + Cr + {i- m 



1-1 



(394) 



Next we do the induction step. Note that the first term on the right side 
of (j378|) can estimated using the maximum principle for a corresponding 
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parabolic equation without source term. Next we use the fact that the 
fundamental solution is positive (nonnegative would suffice at this point). 
Indeed, since rj. solves a uniform parabolic problem with constant second 
order and bounded first order coefficients r'~^ which are Holder continuous 
with respect to the spatial variables uniformly in time, we have 



ri > 



(395) 



(cf. Thm. 11 of chapter 2 of [3] for a proof). Next we estimate the functions 
Sl- We consider the right side of equation (|376|) term by term. The essential 
estimate is with respect to the supremum norm |.|o which we do first. We 
have 



PiIr 



dxi 



Knix-y))E 



j,k=l \ dxj dxj. 



{l-l,y)dy ^<PlCK{Cr + lCr), 

(396) 



where we use ()389p . Next we have for the second and third term on the 
right side of (|376p the estimate 



En 
7 = 1^ 



r,pA — l 



+ 



En 



r,p.l~l dr- 



dxi 



< 2npiC[~^Cr 



(397) 



For the fourth term on the right side of ()376p we have the estimate 



+ 



d 



Kn{x-y)) Elk=i 



dxi 



dx,, 



{I - l,y)dy 



(398) 



< piCk (C, + ICr) + piCk {C[-' + IC[ 
where we use 



i-i 



ax I, 



(399) 



< 



dri 



dv 



r,p,l — l 



dxv. 



(^-i,y). 



Finally, we have the estimate 



i-i , iW-i 



j,m=l 



dxi 



dXr, 



{I - l,y)dy 



(400) 



<PiCk[C[-' +IC[- 
Summing up the four estimates (1396]) . (13971) . p98ll - (HOOj) we get 
\Sl{l-l,.)\o 



<PiCk{ 2&f^ + 12C[~' + 2Cr + l2Cr ) + 2npiC[^'C, 



(401) 



84 



Considering first derivatives we get 
I^K^-l,-)lo,i 

, , (402) 

< pi2Ck (2Co~2^ + /2C^-2^ + 2Cr + l2Cr) + AnpiC'o'^Cr, 



where we note the generosity that we used of the constant Cr in the estimate 
(|40ip . Hence some choice 

- —( — r^i TT^ \ 

4 (2Ck f 2Co72 + ^2C^72 + 2Cr + l2Cr \ + AnCly-^Cr 
leads to 

\S\{l-l,.)\,,i<\. (404) 
Note that (j403|) can be put into the form (|322p above, i.e. ino the form 

Pi < 77 \ — 7\ 5 (405) 

Q ((^172' + Cr)+l [C[-^^ + Cr)) 4C2 

where Cr > 1 without loss of generality and and the generic constant C* ab- 
sobs all the constants which depend on dimension n. Note that the constant 
u is implicitly in Cr- Next we apply 

Lemma 2.20. Assume that pi is sufficiently small. Then for 1 < i < n, 
and for (r, x) € D^'^l we have 

4>\{s,y)T[{T,x-s,y)dsdy < -- [t - {I - I)) , (406) 
l-l 4 

and for (r, x) € [/ — 1, /] x we have 

4>\{s,y)T[{T,x-s,y)dsdy>-(T-{l-l)). (407) 

l-l 4 

Proof. Note that for (r, x) G D''^} means that 4)\{s,y) You check first that 
this follows for the heat kernel instead of F^. Then you use the Levy expan- 
sion representation for FJ, and show that for pi it is a perturbation of of the 
heat kernel. □ 

Similarly we have 

Lemma 2.21. Assume that pi is sufficiently small. Then for 1 < i < n, 
and for (r, x) G D''~1 we have 

Sl{s,y)Ti{T,x;s,y)dsdy\ < 1 (r - (/ - 1)) . (408) 
l-l ^ 
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Next we assume that pi is small enough. Note that we assume inductively 
that the coefficients r'~^ are bounded by Cr (even with repsect to the |.|i^2 
norm on (/ — 1, /] x M"- hence saying that pi is small enough means that piCr 
is small enough such that the proof of the preceding lemmas can be applied. 



Next for (r, x) G D^j^ I we have 



|r^(r,x)| < sup(^_^)g[,_^,]xMn 



c2~ 



Ii-i Ir" {^(^ - + (PiisiV)) ^[{T,x]s,y)dsdy 
Similarly, if (r, x) G D^Z] we have 

r\{T,x) > 4„ r^r^{l - l,y)T[{T,x-{),y)dy+ 

/r" ('^'(^ - + 4>\{s,y)) T[{T,x]s,y)dsdy 
> -C^r- 

Furthermore, for (r, x) D^z} U D^^} we have 

|r:-^(r,x)|<^. 



and, hence, (since > 2) we have 

\r\{T,x)\ < SUP(^^^)g[;_i_;]xRn 



C9 



+ 



ji-i /iRn {S{1 - 1,2/) + 0^(5,2/)) K{T,x;s,y)dsdy 
Hence, we get for I > 2 

sup \r\{T,x)\ < sup \rl{l - l,x)\ < C° 



<C0. 



(409) 



(410) 



(411) 



(412) 



(413) 



for all r G [/—I, /]. Next we look at derivatives. Note that the functions ip^ are 
designed in order to control the growth of the functions r- themselves, not the 
derivatives of that functions. Recall that the functions r'f^(l — 1, .) = r'(/ — 
1, .) are the first order coefficients of the linear parabolic equations which 
determine the fundamental solutions F^. Hence, in order to have classical 
solutions r- it is sufficient to control the |.|o,i of the functions 'r'f^{l — 1, .) 
for alH > 1. Nevertheless, we can control the growth of the derivatives up 
to second order since second order derivative estimates can be reduced to 
first order derivative estimates and the supremum norm of the derivatives of 
r-(/ — 1, .) can be estimated in terms of the supremum norm of itself. For 
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the first spatial derivatives we use properties of the fundamental solution. 
For the second spatial derivatives we use the adjoint in addition. Recall 
that we use the constant C* -which depends on n essentially- generically. 
We show 

Lemma 2.22. There is a constant C* > such that the estimates 

sup \rlj{l,x)\ < C: sup \ri{l - l,x)\, (414) 

and 

sup \rl,4l,x)\ < C: sup \ri{l - l,x)\ (415) 

hold. Here, recall the notation 

d d"^ 
r\,j{r,x) := -Q^r\{T,x), r[^.fc(r,x) := -^-^j\{t,x) (416) 

Proof. For the first derivatives we have the expression 

£-j\ir,x) = f^n rl-\l - l,y)Ti^^{T,x;0,y)dy 

+ /m" i^K^ - 1, y) + y)) ^i,ki^^ ^; y)dsdy 

Next we consider the fundamental solution F^. It has the Levy expansion: 
F^T, x; s, y) := N\t, x-s,y)+ [ [ N^{t, x; a, C)M^, s, y)dadi, (418) 



(417) 



where 

N\T,x-s,y) = ] ^exp f-J^^J^) , (419) 

and (pr is a recursively defined function which is Holder continuous in x, i.e., 

oo 

Mr, X] s,y) = Y, mN%{T, X- s, y), (420) 

m=l 

along with the recursion 

mN%{T,x;s,y) = LJN'{T,x;s,y) 

(421) 

— Pll^j=l'j dXj ^ 

(L'-iVOrn+i(r,x) := /^7^ {UN^{r,x;a,0)^L^'N'i^,C;s,y)dadC 
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The lemma may be proved for Ni instead of in the representation ()417p 
first and then the proof may be extended to the correction terms involv- 
ing fj N''{T,x;a,(,)(l)r{o;^;s,y)dad^ using the classical Levy expansion 
estimates. The proof may be obtained in sub-timesteps, i.e. the first the 
norm 



sup HM - l,x)\ < C* sup \rl{l,x)\ < C* sup \ri{l,x)\. 



Then we may use 



(422) 



Afl (n- -T- « 7/^1 — i^~y)k 1 „„„ / \^-y? 



(423) 



and 



Apii>{T-s) 



exp 



\x-y\^ 



mplv{T-s) J (x-J/)fe 



1 



■ exp 



WPiv{t~s) 



< C 



(424) 



for some C > 0. This leads to the first statement of the lemma. Finally we 
use the argument above that 

sup \rl{l,x)\ < sup \rl{l-l,x)\. (425) 

For the second derivatives we use partial integration and the adjoint 

^) = kn r\jHl - l,y)Tr'Ur, x; 0, y)dy 
+ //-I /r" - 1' y)^r!m{T, x; s, y)dsdy (426) 

+ Imp ^i,k(^^ y)'^r'!m{T, x; s, y)dsdy. 

It follows that the estimate can be reduced to the estimate for first deriva- 
tives. Since C* is generic we may replace (C*)^ by C*. 

□ 

Next we observe that we have an upper bound for the integral magnitudes 
for r' and its derivatives is preserved, i.e., we want to show that 



Y 1^ 

^ I dxj 

j,k,m,p=l 



dr 



dXr. 



{t, y)dy <Cr + ICr 



and 



E 

j,k,m.,p,q=l 



d'^rl dri 



dxjdxq 



dxr. 



{t,y)dy < Cr + lCr 



(427) 



(428) 
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for 1 < k < n, where we have assumed this upper bound for I — 1 inductively. 
The estimates ()429p and ()430p . which are of a form which is closely related 
to the integral term of the Navier-Stokes equation in Leray projection form 
- follow from the H^- and H^-type estimates 



E 

j,k,m,p=l 



dxj 



+ 



dXr. 



{t, y)dy < Cr + ICr 



(429) 



and 



j,k,m,p,q=l 



+ 



dXr. 



(r, y)dy < Cr + ICr 



(430) 



Therefore we assume this inductively. Similar for the function v, 
assume inductively that 



r,p,l 



j,k,m,p=l 



dx^ 



uum 



dXr 



T,y)dy<C["i + {l-l)C[-^'. (431) 



We 



and 



/ 



E 

j,k,m,p,q=l 



1 




1 dXjdXq 


dXp 



{T,y)dy<C[-,' + {l-l)C[-, 



i-i 



Again this follows from estimates of type 



(432) 



, E 

j,k,m,p=l 



dXj 



+ 



UUm 



and 



/ 



E 

j,k,m,p,q=l 



dx j dx Q 



+ 



dXr. 



UUm 



{T,y)dy<C[~^^ + {l-l)C[-^\ 

(433) 



dXr 



{T,y)dy<C[-i + {l-l)C[-^ 



i-i 



(434) 

In the next substep we shall see that there is an upper bound Ci^2 for all 
C\ 2 which is independent of the time step number /. 

In order to find these upper bounds (|429p and (|430p we go back to the 
representation of r\ in ()378p . The first term in ()378p is an integral of the 



initial data r^~^ times the fundamntal solution F',. The second order coeffi- 
cients piv clearly have a uniform upper bound and the first order coefficients 
r^~^ (which are also data from the previous time step) have a uniform upper 
bound by inductive assumption. We also know that the first spatial deriva- 
tive of the fundamental solution is locally and globally integrable for fixed 
X G (especially, the integral with respect to y is in L°°). We may use 
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the adjoint for the first and second derivatives in order to estimate the first 
term on the right side of the representation in (j378p in terms of expressions 
of the for ()429p and ()430p . Involvement of the fundamental solution FJ, leads 
to a factor C* = 2Cp, where we may use our generic constant C* since Cr 
depends on dimension n essentially and is defined similar as in the first step 
of this proof, i.e., 



sup 



{T,x)eDl Jl~l JR" Irl^ jSa 



sup 



r' I + ir' • 

r I ' \ r.t 



dyds =: Cr- 



(435) 



Again local integrability and Gaussian-type a priori estimates lead to the 
conclusion that Cr is finite. Note that we can estimate products by sums 
of squares in order to estimate mixed terms. There is one difference when 
we estimate supremum norms and integral norms of r\. Note that we have 



defined 



+ 



in terms of the functions r\ ^ and t)'"''''^ ^ such that 



the i?i-and H'^ -norms of cj)\ are bounded by the sum of the ff^-and H^- 
norms of r'^"*^ and Moreover, by the very definition of 0' this holds 

also for the function (note that pi is small) 



= /jg„r( ^{l - l,y)Tl{T,x]0,y)dy 

+ fj^-i /iRn </>'(«, y)r'.(7", s, y)dsdy. 



(436) 



More precisely, since pi is small we have 



dRi 



dxi 



(r, y)dy < 



dxi 



+ 



dv 



r,p,Z — 1 



dxi 



(437) 



{T,y)dy 



Hence, we get for generic C* > 



p=l 



dri 



dx-j 



n 
dr'- 



dxr, 



Q El 



j,k=l 



dv 



r,p,Z — 1 



dxi 



(r, y)dy < 



j,k=i 



+ h 



(438) 



< CI ( Cr + 1)C, + C[f + l)Cr/ ) + /: 
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where 



h = lun Ej,fc,m,p=i /ffin Siil - 1, y)Tl j{T, x; s, y)dsdy 
^ //-I /r" ^rai^ ~ 1' J')rr,p(T, a:; s, y)dsdydx 



+ kn y Elk,m,p=i /r- ^1(1 - l,y)T[ j{T, X- s, y)dsdy 

+ ( /z-i /r" ^mil - 1. a;; s, y)dsdy^ ^ 

and, similarly, for the first spatial derivatives we have 



En 
j,k,rn,p,q=l 

' Z-^j,k,m,p,q=l 
' Z—/j,k,m,p,q=l 



{t,y)dy < 



,2,,'-,P,i-l 



+ 



dv 



r,p,l~l 



+ 



dx„ 



dx„ 



< C* {Cr + {1- 1)C,) + h 



where 



(439) 



(440) 



-^2 — /iRn Z^j,fc,m,p,g=l /ri 'S'fc,q(^ y)^r,ji'^^ ^' y)dsdy 

^ //-I /r" 'S'm,g(^ - 1, y)rr,p(T, x; s, y)dsdydx 



p (r, x ; s , y)dsdyj j dx 

Note that squared intergals produce mixed terms such as products (of 
derivatives) of and r' which can estimated by sums of squared terms 
again. This leads to estimates of Ii and I2 involving L^-norms and L^-norms 
in terms of functions of the previous time step. Note that all terms of 
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have the factor pi. For Sl we have 

m - 1, .)lo,i < \pi /m. (£-Kn{x - y)) E • .=1 ( ^ t ) - 1, y)dyl^^ 



I I Tr^n l-l dv.-' I I r 

+ \PI Ei=i rj 1 0,1 + k« Ej=i 



i9a;, lo,l 



+\pi [m-Mx - y)) Elm=i ( ^^^^ ) a - 1, 2/)dy|o,i 



< p/Q ( a + iCr + q;/ + ic^. 

(442) 

Similarly, using L^- estimates -which we have from the previous time step 
inductively- we get 



\si{i - 1, < PiC* (Cr + iCr + c[-;,' + ic[-^ 

Hence we have (generic C*) 

h + h< PlCl {Cr + ICr + C[-^ + IC{-^) . 

Hence, we get for generic C* > 



(443) 



(444) 



fc,m,p=l 



dXri 



[r, y)dy 



(445) 



The estimates 



/r" Xyj,fc,m,p,ij=l 



(r, y)dy 



(446) 



and the related H^- and ff^-estimates 



j,k,m,p=l 



dxi 



+ 



dxji 



ir,y)dy 



(447) 



^ ( Cr + ICr ~l~ C*]^ 2 ~l~ IC^ 2 
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and 



k,m,p,q=l 



+ 



dr' 



dx„ 



{t, y)dy 



(448) 



are obtained similarly. 

ad iii) We have defined the functions 



• and Ti in terms of the functions 



^'"''-^ and r'-\ 
data and with r\ 



In the case / = 1 we start with v 







hi, the initial 

0. We have to estimate f ,['''' ""^ first in this special 
situation. Then we proceed with the case I > 2. In order to finish the first 
induction step / = 1 we have to consider the growth of the local solution 
of v^'^'^. Note that for this first step this equals the local solution of the 
Navier-Stokes equation v'''' because we have = 0. Therefore at time step 
^ = 1 we may apply the machinery of the first step of this proof. Hence for 
^ = 1 we consider the functional series (v^''^'"'^)^, where 



p,k,l 

vy 



m=l 



p,m,l 



(449) 



This series converges to the local solution v^'' of the Navier-Stokes equa- 
tion. Note that our symbolism allows for two alternative expressions for the 
members of this series. It can be denoted by tip''''^'^ 



p,k,i 



p,k,l 

vy , or 



by the values of the map Fi of the first step of this proof. Recall that the 
series (jl80p can be generated by an iterative application of the map 



Fi :f 



(450) 



as described in step 1 above in case of general time step number I starting 



with the initial data h. Then the function -Fi(h) 
the equation 



T^h,p,i _ .^p,o,i satisfies 



j=l "'j dXj 



PI 4. [irKr^ix - y) ) Elk=l ( Hfe ) y)dy 



(451) 



where 1 < i < n. Let Fg be the fundamental solution of the equation 



dv. 



h.p.l 



En 
7=1 



Q2h.p,l 



I V^n I, '^'"i 
+ PlY. -,=1^3-3, 



h,p,l 



Let C|^, C^, and be some constants such that 

l^ilo < C'hi l^ili < C\ , \hi\2 < C12 



(452) 



(453) 
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for all 1 < i < n. In terms of the fundamental solution Fq of the equation 
(j452p the components v^'^'^ of the solution of (|45ip have the representation 

vf'°'^(r,x) = hi{y)rl{T,x;0,y)dy 

+ Io PiJ^n [£-Kniy-z)^ Elk=i {wjSk) is,z)Tl{T,x-s,y)dydzds. 

(454) 

From the maximum principle (of. also Corollary 8.1.3. of |TT]) we observe 
that 

hi{y)Tl{T,x;0,y)dy\ < \h\o (455) 



for all (r,x) G [0,1] x M". Hence, the solution v'^''''^ = v^'O'^ = v'^'P'^'^ of 
equation ()459p satisfies 

< CO + \p, {£-K^{x - y)) (Hit) {r,y)dy\Cr 

< CO + piCrCi,2 < Ci,2 + piCrCi,2 

(456) 

where the definition of Ci^2 in terms of the function h is used. We are quite 
generous with the use of the bound C1.2. As we described in the first step of 
this proof we may use the adjoint of the fundamental solution Fq and shift 
derivatives, we get the estimate 

\vf'\2<Cl{Cl, + piCi,2) (457) 

Choosing 
we get 

\vf\,2<C^{Cl + l) (459) 

Since = we have convergence of the local scheme on [0, 1] x where 
we choose C^ = and Z = 1 in the definition of pi of the second step (or first 
step) of this proof. We may use pi with 

in order to get an estimate similar to ()459p for the vector-valued function 
V. Adding the estimate of the correction v'''^ — v'''^'^ = l^t^f |i,2 < 1/4 

from the first step of this proof we get 

|v^'^|?,2 = Er=i \<\,2 < {Cl, + 1) (461) 

where the constant Cr > 1 is from step 1 of this proof. Now we have to 
show that this constant Ci 2 = C,. is preserved for all I > 2. It is clear that 
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our definition of tlie constant Cr fits for the first step Z = 1, since we have 
chosen = in this first step. Indeed the constant Cr is determined in 
the induction step with respect to the time step number I. Furthermore we 
have determined p\ for the first time step according to the convergence rule 
of the local scheme established in step 1 of this proof. Let us assume that 
'^r,p,m g^j^^ j.m j-^g^yg ^gg^ Constructed for m = 1, — 1 as classical local 
solutions on the domains [m — l,m] x M". Furthermore assume that we have 

\vl''no<C and |<'^''"|i,2<Ci,2 (462) 

for all 1 < m < / — 1 and 1 < i < n. Furthermore, let us assume that 

br™(r,.)|i2 <C*(Ci,2 + mCi,2) 

(463) 

|<'''''"(r,.)|H2 <C*(Ci,2+mCi,2) 
for all r G [m — l,m]. 

Remark 2.23. Since we consider C* in the definition of the constant Ci^2 to 
be generic we note that 

|^r,p,m|n^ < Ci,2. (464) 

Let us go back to the local Navier-Stokes equation on the domain [I — 
1,1] X which we wrote in the form ()151|) . Let us repeat this equation here 
for the convenience of the reader. We have 



- PIT' Ej=i -i^ + pi Ej=i -t^ = 

, v^'/'''(/- 1,.) = v'"'/'''-i(/ - 1,.), 



(465) 



where 



4 = r[r - Pl^^rl + pi Y:]=i r^jo^ 

~pi Ir- {ik^nix - y)) E",fe=i (£^£t) i^^y)dy 

+Pi Ej=i + Pi Ej=i vf' Q^. (466) 

+pi /r" {£-Kn{x - y)) Elk=i i^^y^'^y- 

Note that the right side of this equation involves the function v'^''''' which 
is not known. However, the function r' is known by our construction in 
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the previous substep, and determines together with the function v'''''''~^(/ — 
1, .) the right side in the equation for the first approximation v''''''^'' of the 
function v'''''''''' . We may estimate the growth of the function v'''^'"'' first and 
then estimate the growth of the function v^''''' by estimating the difference 



5V^P'^'^ (467) 



k=l 



which we know from the estimates of the previous steps of our proof. We 
have 

r,n I rjO ^:P^O,l , n r,n,0 I 

OT/ V^n d^v-"^' r,p,l-lrj ^ n Of , ,Z,0 

- PI'' Ej=i + PI Ej=i Vj - h = 'I'i 

^r,p,0,l^l _!,.)= V'P'^-^{1 - 1,.), 

where 



(468) 



-Pi /r" (m-Kn{x - y)) E",fc=i (S^St) (J,y)dy 

+Pi Ej=i + PI Ej=i 9^ (469) 

+pi /m" {m-^n{x - y)) Ei,fc=i (^ ^"C, ) - 1. 

We have determined r' in the previous time step such that the first three 
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terms on the right side of ()468p can be replaced. Indeed recall that we have 

En l-l dv^'''''^'^ _ r,p,l-l dr'~^ 

j=l dxj Pi 2^3=1 dxj 

-PI /m" {irK^x - y)) Elm=i [^^'-^) (1 - ^^y)dy + 4- 

(470) 

We may rewrite the function -i/;-'*^ of (j469p . i.e., the right side of the first 
equation in ()468p in the form 

= ^\ + (v.;-" - <Ai) , (471) 

where 

^f-4 = Pii:%M--Y')l^. 

-PI /r" {-ir^nix - y)) Elk=i (Sjfe) iT,y)dy 

+pi /r" {mi^n{x - y)) Elk=i " 

+Pi T.%M - -r^)^ + PI (si - '4^) 

-2pz /m. (^i^n(x - y)) Elk=i ( y)- '4^(1 -I, y?j '-^{l - 1, y?j dy. 

(472) 

Indeed, this difference of ^''"^ and 0' may be obtained by subtracting equa- 
tion (|157p from equation (|156p which we computed in (|158p above. Note 
that all terms on the right side have the factor pi . Estimating this difference 
boils down to estimating the difference r' — r^~^. In order to control the 
growth of the functions tip'''^'' the estimate should be such that the func- 
tions 4>l dominate the differences ipl'^ — in critcal areas where the moduli 
of the functions f ^ exceed a certain level. Summarizing we first estimate 
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v^'^' ' , where we rewrite ()468p in the form 

^ vr.P.O.'(/ - 1,.) = v'^'''''-l(/- 1,.), 

(473) 

and then we estimate f — v'j^''^'^'\ where we use the second step of this 
proof, and the fact that the difference ip\ — is given by 

I sr^n I ( dv'^''''' Si)'''''''"^ \ ( r,p,l r,p,l-l\ dr'. 

+pi E,=i i (^^r fe— j+Pi Ej=i - ^/ ) w- 

-'^pikr. {£-Kn{x - y)) Elk=i (Sj (%r " ^^as^)) 
+ft /r" {£-Kn{x - 2/)) E",fc=i (^^f ^5^) i^^y)dy 

-PI /r" (^^"(^ - y)) E",fc=i ( ^''k', ) {T,y)dy. 

(474) 

Indeed this is essentially estimated in the second step of this proof once the 
first difference is estimated. 
Recall that 

Ir^lo < C^-, \ri\i,2 < C], and |r^|i,2 < (475) 

We see that we have to estimate the difference br\ = rj — r^'j^^i^ — !> •) in 
order to get the estimate of that difference. Well, we an afford to do a rough 
estimate for this difference in the form 

\Srl\i,2 = |rj - r'r\l - 1, .)|i,2 < 2^- (476) 
Next we estimate the difference 'i/'i''' — <Pi- From equation ()472p we get 



IV'-'" - 4\o,i < Pi [^nC^ + ^u^CkCI + 4naCo,2 + SCi^n^aCo.a) • 

(477) 

We get 

- <^'lo,i < \ ■ (478) 
where we choose ^ 



98 



Note the additional factor Cr and recall that C* is a generic constant 
dependent on dimension n. Recall also that we have determined Cr in 
terms of the initial data function h above. Hence the scheme defined is 
global. We can now finish the proof. In order to control the growth of the 
functions v''''''' and the function at time step / we fir control 1) the 
growth of a linearized equation for v^''''*^'', and 2) ensure that the correction 
■^rp,i _ ^r,p,o,i _ §yr^,p,k,o gj^all enough by choosing pi appropriately. 

Let us look at the linear approximation v'''^''^'' of the function v''''''' first. 
From equation (|17ip we get the representation 



r,p,0,l 



(r,x) = 4„ v''P''~\l - l,y)rJ, o(r,x;0,y)dy 



(480) 

+ Iru i'i'^is, yWioir, x; s, y)dyds 
where FJ, g is the fundamental solution of 

r\ r, 0,0,1 i-iO r, 0,0,1 ,_, r-t r,p,0,l 

-^-pi^Ej=i^^ + PiEj=iVf' {l-l,.)-^ = 0, (481) 

and ^1''^ is as in (|17ip . We may rewrite 

^''''O'V, = v-'P''-\l - 1, y)Tl^o{r, x; 0, y)dy 

(482) 

+ /«-i/r" ((V'l'°('S,y) - (^^(s,y)j +4(8, y)j TI q{t,x; s,y)dyds 

where is constructed above and the difference ■0'''^ — is as defined in 
(j472p . Applying the maximum principle we may estimate the first term on 
the right side of equation ()482p by the supremum (indeed, maximum) of the 
initial data. We have for all (r, x) G [Z — 1, /] x M" 

K''''°''(r,x)| < sup(^^^)g[;_i^;]xKn J^„v''f''^~^{l - l,y)TlQ{T,x;0,y)dy 

+ /«~i/r" (i + '/'i('S,2/)) Tlo{T,x;s,y)dyds < sup^gjjn K'''''"^(/ - l,x)|. 

(483) 

For this part the reasoning is the same as before in the case of the function 
We use the maximum principle in order to estimate the term 

v'^P^'-\l-l,y)ri^o{T,x;0,y)dy (484) 

and then consider the different cases making up the definition of 4>\ together 
with the a priori estimate of the fundamental solution F^ q involved. Next 
we apply the second step of this proof which gives 

m=l 
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(486) 
(487) 
(488) 

(489) 

estimate for r' in the 

2.4 step 4: Global existence of classical solutions v'' and v 
resp. 

The functions 0' are bounded functions with supremum less or equal to 1 
which vary with the time step number in general. At each time step / > 1 we 
defined them on the domain (1 — 1,1] x W^. The functions v^''' and r' are only 
weakly differentiable with respect to time across the points / = 0, 1, 2, • • • 
in general. They are also Lipschitz. Furthermore, in the construction of 
each time step the functions (pl are only Lipschitz with respect to the spatial 
variables in general. However, this is sufficient in order to find that the 
solution = v^'^ + r constructed is classical. Consider the function V'^*'^ 
constructed at time step / > 1. The components 1 < i < n have the 

representation 

oo 

^r,p,l ^^r,pAl ^^^^r,P,k,l^ (490) 
k=l 

1 r r,p,k,l r,p.kl r,p,k—l,l j , • i • i i- 

where dv^ — '^i determined successively by linear 

Cauchy problems with zero initial conditions. These linear Cauchy problems 
have bounded classical solutions with 

Q^yr,p,k,l 

lim (r,2;)=0. (491) 

rii-i dr 

This implies that the regular behavior of the function v*"'^'' with respect to 
time is determined by the behavior of the function v^^p^^^^ as t ], l — l and the 
initial data (resp. the final data of the previous time step) ^/^'P'^'^{1 — 1^ .) = 
v^'P''-i(Z — 1, .). We conclude that the function v'''^ is uniformly bounded 
continuous with respect to time. Furthermore it is Holder continuous with 
respect to the spatial variables uniformly with respect to time. Similarly, 
the function r — t- r(r, x) is only weakly differentiable at the integer values 



Note that 

\vr lo < \vr I - 2' 

We conclude that 

|l,2 S (-/1,2- 

The reasoning that 

|<'^''|h2 <Ci,2+/Ci,2 

reduces to the reasoning that 

The estimate ()489p is obtained analogously to the 
second substep of this third substep. 
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I € N. Especially it is uniformly continuous with respect to the time variable 
r and it is Holder continuous with respect to the spatial variables uniformly 
in T. Hence we conclude that the solution v'"''''^ on [/ — 1,/] x M" of the 
Navier-Stokes equation in transformed time coordinates r, i.e., the function 

shares these properties. Hence, if we consider the family of fundamental 
solutions r(J'' of the equations 

"'^^ pz.AC' + pz^<f^ = 0, (493) 



7 = 1 



then we observe 1) that all exist in their Levy expansion form since the 
coefficient functions f^'' are uniformly continuous with respect to time and 
Holder continuous with respect to the spatial variables 2) we can build a 
global bounded continuous coefficient functions Vj : [0, oo) x M" out of the 

local coefficient functions f , where the local restrictions of Vj to the domain 

[I — 1,1] X M" are equal to the functions f . The different coefficients pi 
in the family of equations ()493|) are artefacts of the time transformations, 
of course. We can get rid of them by transforming back to original time 
coordinates 

r ^ t(r) := piT, if r G [/ - 1, /] (494) 

for all Z > 1. Note that we can consider the function Vj as given because 
we have constructed it- the turn to the fundamental solution then gives us 
global classical solutions. Furthermore, the fundamental solutions of 

j=l 1 

exists on the transformed domains and in original coordinates, and we can 
build global bounded continuous coefficient functions Vj : [0, oo) x out of 
the local coefficient functions Vj in original time coordinates, where the local 
restrictions of Vj to the domain [/ — 1 , /] x are equal to the functions Vj . 
The related global fundamental solution can be obtained also by successive 
use of the the relation 

T{t,x;s,y)= T{t,x;si,z)T{si,z;s,y)dz (496) 

for t > si > s, and where si will run through the time step sizes Ylm=i Pm 
in orignal time coordinates. We conclude that the fundamental solution F^, 
of 

vAF. + ^.,§i = (497) 



dt ^ dxj 
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exists on [0, T] x M" for arbitrary T > 0, where Vi is the function that equals 
v\ on the domain [/ - 1, Z] x M" for all I > 1. Now for all t G [0, T] and x G M" 
we have the representation 

Vi{t,x) = J^„hiiy)T,it,x;0,y)dy + fe^) (^'^)x 

(498) 

xK^i{z — y)r^(t, x; s, y)dsdydz 

for 1 < i < n. Hence we have that Vi G C^^'^ ([0,r] x W) for all T > 0. 
Furthermore the latter representation and the standard a priori estimates 
for classical fundamental solutions together with our assumptions on the 
initial data hi show us that we have 

Vi{t,.)£H^ (499) 

for all t > 0, where we use Young's inequality. □ 



3 Regularity, uniqueness, and extensions with ex- 
ternal forces 

The existence of global bounded classical solutions is essential in order to 
prove further regularity and uniqueness. Furthermore this is essential in 
order to extend the proof to equations with external forces, and to study 
asymptotic behavior. Maybe it is well known that bounded classical solu- 
tions lead to full regularity. However, we shall show this for the sake of 
completeness. 

Next we prove regularity, i.e. we prove that for 1 < i < n we have 

G C°°([0,oo) X M"), (500) 

and for all t G [0, oo) and s G M 

Vi{t,.) e H' {R'') . (501) 

One way to do this starts from the representation of the solution in (j498|) 
Note that we may write 

Vi{t,x) = v'^it,x)+p,iit,x), (502) 

where 

^iit, x) = /ig„ hi{y)T^{t, x; 0, y)dy (503) 
solves the Cauchy problem for the Burgers equation, i.e., 

\ (504) 
I v(0,.) = h, 
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and p^i describes the gradient of the pressure, i.e., the solutions to the Poisson 
equations 

-Ap. = E • .=1 £- (afc^.) i^^k) (505) 

for 1 < i < n. Note that in terms of the Navier-Stokes solution function 
V the first equation in ()504p is a linear parabolic equation for for each 
1 < i < n. Hence a simple strategy for regularity is the following: start 
with some known regularity of v and consider ()502p . Then use regularity 
theory for linear parabolic equations in order to prove more regularity for 
the summand v^. Then starting with some regularity for the right side of 
(|505p use elliptic regularity theory for Poission equations in order to get more 
regularity for p^i. This gives more regularity for v. Iteration of the process 
leads to full regularity. In a first step we look at the Poisson equation for p 
itself. Our main theorem tells us that we have 

(Well, it tells us that the right side is even in H^). Hence the convolution 
with the fundamental solution n of the Poisson equation is a well defined 
locally integrable function, and we have 



(507) 



Here recall that an estimate for the term in ()506p can be reduced to an 
estimate. From our main result we know that for each time t G [0, oo) 
the right side of ()507p is in C". We get more regularity of this distributive 
solution for p from the standard result 

Lemma 3.1. Assume that k > 0, and Q is an open set in M". Assume u is 
a distribution solution of 

Au = f (508) 
where the data f G (7^+" (j^) along with a G (0, 1) and k > 0. Then u G 

Hence we have p{t^ .) G C^^" in the first step. Furthermore, we have 

Lemma 3.2. Suppose that for some a > we have u G C^"*"", where u 
satisfies the Poisson equation 

Au = f (509) 
. Then for all s > and p G (1, oo) we have 

f G H''P ^ue H'^+''P. (510) 
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This may be combined with Sobolev type embeddings for Zygmund 
spaces CI, i.e., 



Lemma 3.3. For s > and p G (l,oo) we have 

11 

H'^^CCl, ifr = s--. (511) 
P 

This means that p{t, .) E H"^'"^ C C^"^" where p may be large. From 
our main result we have Vi{t, .) G H'^'P for p > I. Hence Vi{t, .) G C2+" n 
i?^'*', and these are the first order coefficients of the equation for v^. Then 
differentiating the equation (|504p and using classical theory we get 

v^iit, .) G C^'^ n H^'P (512) 

for all t > Hence, we have 

Vi{t, .) £ C'^'^ n H^'P (513) 

for all t > and p > 1. Iterating this argument we get full spatial regularity, 
i.e., 

Vi{t,.) G C°° for all t > 0. (514) 

The Navier-Stokes equation itself then tells us immediately some regularity 
of the first time derivative. We may then consider the time derivative of the 
Navier Stokes equation and consider this as a Navier-Stokes type equation 
for ^ (with a certain source term and a potential term). We may iterate 
the argument above and get 

^Ui^ {t, .) G for all t > 0. (515) 

Another way to prove regularity is the following which uses both the Leray 



projection form and the original form of the Navier-Stokes equation. We 
shall use the following We have Vi G C^'^ for all 1 < i < n. Hence Vij{t, .) G 
C^, which implies for all t G [0,oo) that 



J2 {vj,kVk,j){t,.)GClcC". (516) 
j,k=i 

Next consider multiindices a = (ai, • • • , ,a„) and the multivariate partial 
derivative of order a, i.e., 

-■■ =: , . (517) 



dx\°'\ 
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The spatial derivative of order a with \a\ > 1 of the Navier Stokes equation 
leads to 



n n y \ 

Vi,a,t - v/\Vi^a + ^ VjVi^j^a = -P,i,a " ^ '>Jj,aVi,j - ^ i j Vi^pVij^a-p- 
j=l j=l 0<l3<a ^"^^ 

(518) 

Next, assume inductively that 

Vi^fS S Cf['^, for all /5 < a, and I < i < n, (519) 

and 

Vi^l3{t, .) G L^, for all P < a and all t G [0, oo) (520) 
for all /3 < a, and 1 < i < n. Then 



\0</3<a ^'^^ 



(t,.)GQcC", (521) 



uniformly in the time variable t. Hence, we have the representation with 
fundamental solution T 

Vi,a{t, x) = /jg„ hi^a{y)^{t, X, 0, y)dy+ 
- !oP,iA^^ y)r(i, x, s, y)dyds 
+ E?=i {'^j,aVi,j) (s, y)r(t, X, s, y)ds(iy 



(522) 



+ (Eo</3<a {s ,y)T {t, X , s ,y)dsdy , 

where F is the fundamental solution of 

--^Ar + ^.,— = 0. (523) 



dt ^ dx. 



Hence, we have 



Vi,a e C^''' (524) 

Regularity with repect to time can be treated similarly using regularity with 
respect to the spatial variables. 

1 2 

Next we have uniqueness. Assume that vi,V2 G C^' are two solutions 
of the Navier-Stokes equation. Then from the basic energy estimate 

sup -?;2|o < |(t^i(0,.) -r;2(0,.))|oexp f / \Vv2\L'^dt] (525) 

0<t<T \Jo J 

we get uniqueness. This is an application of Gronwall's lemma and can be 
found in standard texts. We mention it here for the sake of completeness. 
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Finally we mention extensions to equations with external forces and asymp- 
totic behaviour. The solution of the Navier-Stokes equation wih external 
forces fex can be represented as 



Vi{t, x) = hi{y)T{t, X, 0, y)dy 



lo Ir" p,i(s,y)r(t,x,s, 



y)dyds 



(526) 



+ /o /m" /e2:(s, y)r(t, X, s, y)dsdy. 



Since the gradient of the pressure is bounded in case fex = we see that the 
global solution converges to as 1 1 oo. 

4 The algorithm 

The construction of the solution above can be extended to boundary value 
problems straightforwardly. We consider first initial-boundary value prob- 
lems and second initial-boundary value problems. The observation is that 
the global existence for these problems boils down to the existence of related 
scalar first initial-boundary value problems and second initial-boundary value 
problems of parabolic type. Furthermore, the scheme proposed has linear 
subproblems of parabolic type which can be computed explicitly. Expan- 
sions of this form are considered in [8] . Note that there is a difference to the 
Taylor expansion (operator form) which applies only for affine coefficients 
in general (cf. [T] and 0). Furthermore the subproblems considered here 
have been implemented in the context of weighted Monte-Carlo methods 
in finance (cf. [10], [5], and [3]). The discussion here is rather conceptual. 
Further details of implementation and error estimates will be considered in 
subsequent paper. 

Since we have proved that the solution is bounded we may set up a 
uniform time discretization scheme. The size of the time steps is limited 
by a time step size which ensures the convergence of the local time scheme. 
A lower bound of this time step size can be extracted the global existence 
proof. The step size numbers pi may be increased as time goes by if there 
are smoothing effects due to the strictly parabolic subproblems. 

Since we have proved that the solution fj, 1 < i < n of the incom- 
pressible Navier-Stokes equation is globaly bounded and Holder continuous 
with respect to the spatial variables uniformly in time, we know that the 
fundamental solution r„ of the equation 



exists (given v). Hence, in terms of the solution itself the solution of the 




3 



(527) 
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Cauchy problem has the representation (original time coordinates) 
Vi{t, x) = hi{y)T^{t, x; 0, y)dy+ 

(528) 

//-I /m" Em,«=i {vi,mVm,i) (s, z)K{y - z)T^,{t, x; s, y)dy. 

This representation cannot be used for computation of course, since we do 
not now the solution. However, we have shown that for a time step size p < 1 
which is small enough we may compute the solution in a time -discretized 
scheme where at each time step / we compute successive approximations 
yp,k,i T^iii^ initial values from the previous time step. 

Note that the precise values from the the previous time step are vP'^~^ 
which we do not know except for the case I = 1 where ff''~^ = hi, i.e. equal 
the original initial data (well even these have to be approximated upon 
implementation). Hence, the initial data of the previous time step are some 
given by a function vP'''~^'* which approximates v^' ~^ . The approximation 
of ff'' is then computed by iterative approximations v^'^'^'* for /c > where 



Pkr. {£-Knix-y)) E-.., (^^) ir,y)dy, (^29) 

, v'''°'''*(/ - 1,.) = v^''-i'*(Z- 1,.), 

and for A; > 1 recursively defined functions v^''^'''* are determined by the 
respective solutions of 

av'r o^vr' p,k-l,l,* avr 



^ v^'^'''*(; - 1,.) = v^''-i'*(/-i,.). 

We have to stop after finitely many steps. Hence we shall have 

p,Z,* p,m,l,* /roiN 

<" =<' (531) 

for some m if we perform iterations A; = 0, • • • , m at time step I. Let r^''' be 
the fundamental solution of the equation 

^-^^■^^E^r-^-O- (532) 
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Then the solution t;^''''*^'* of ()529p has the representation 
ff'''°'*(t,x) = vf~^'*{l - l,y)T^f{T,x;l- l,y)dy+ 



Furthermore, let T^^^ be the fundamental solution of the equation 



(533) 



.Arl'^ + j:.f-'''*^ = 0. (534) 

Then the solution ()530p has the representation 
?;f'''''*(t,x) = /ig„ ?;f''"^'*(y)rl'''(T,x;0,y)(iy+ 

/z-i /iRn Em,«=i z)KAy - ^)rl'^T, x; s, y)dy. 

(535) 

This scheme involves the fundamental solutions of equations of type 

J=l J 1=1 

essentially. For computational point of view it is interesting that the solution 
has the locally pointwise valid representation 



1 ( ELiAx2 

pit, X, 0, y) = —f=n exp ' 



fWr \ 4t 



) (^f;^A:(i,^,2/)i'j , (537) 



for j = I,-- - ,n, i.e., the representation is valid on some time interval. 
However, this fits with our scheme since this is defined locally in time anyway. 
Note that we have the coefficients outside the exponential. This implies 
that the first term is damping the polynomial terms (i^ as the moduli of 
the Axj = (xj — yi) become large. The coefficient functions have explict 
representations in terms of the coefficient functions h^. for A; = we have 

(io(t, X, y) = exp ^^(y™ - x^) j bm{t, y + s{x - y))dsj , (538) 



^ k 

dm{t,x,y) = } —dm-k Rk^i{t, 



y + s{x-y),y)s''ds (539) 



with 



Rk-iit,x,y) = ftCk-i + Ack-i + J21^^Ylr-=o (^^V^Cfc-i- 
+ Eibi{oo)£-Ck-i 



(540) 
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If the coefRcients bi are given in terms of bounded analytical expansions 
(finite Fourier series for example), then the functions can be computed 
explicitly (cf. [8]). Note that we cannot use Taylor expansions (operator 
form) as in [1] since complete sets of analytic vectors are difficult to define 
if the coefficient functions are not affine (which is true in case of the Navier- 
Stokes equation) . These analytical expansions have been proved to be quite 
efficient in a different context (cf. [U El [TUl E]). Fluid dynamical models 
in applied sciences have boundaries of course, so let have a look how our 
algorithm can be adapted to these situations. We consider boundary value 
problems which are related to the second initial-boundary value problem for 
scalar parabolic equations - cf. [3] for a classical treatment. We consider 
problems on the domain [0, T] x 0,, where 0, C M" is a bounded domain. 
Let B := {0} x Q, Bt = {0} x n, and let S = dn\{BLI Bt). Then we 
consider the following initial-boundary value problem on [0, T] x $7. Let 
ai : [0, r] X — )• M, and : [0, T] x O — )• M be 2n functions. We consider a 
problem for Vj, 1 < ■« < n, where 

( ^ — j/Y^" _|_ V^n „, . dvj _ 

dt ^ 2^j=l &J l^j=l ^3 dxj - 

I {£-Kn{x - y)) E • ...1 (life) (t, y)dy, ^^^^^ 
-§^Vi{t, x) + ai{t, x)vi{t, x) = gi{t, x) on [0, T] x 5, 
. v(0,.)=h. 

The scheme we proposed for the Cauchy problem can be adapted to this 
situation straightforwardly. Consider a time discretization in transformed 
coordinates and assume thhat v^''~^ has computed for I — 1, where / > If 
I = we set v^'~^ = h which is known. We choose a fixed p independent of 
the time step number /. For Z > 1 we have functions v'''' : [Z — 1, /] x — )• R"' 
on successive domains where the final data of the function v''''"^ are the 
initial dats of the function v''''. Each function v'''' is determined as a limit 
of a functional series (v''''^'')^,. Note that transformation to coordinates 
t = pT leads to a domain [0, -RT] x 17, where R = j^. Having computed the 
k — 1 iteration step of the Ith. time step the problem for v^''^'' is given by 
n linear parabolic scalr problem which are classical initial-boundary value 
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problems of second type. 



^ff'''''(r,x) +ai(t,a;)ff'''''(r,x) = gi{t,x) on [0, i^T] x S", 

vP^^'^{l - 1,.) = v''''-i(/ - 1,.), 

(542) 

where for A: = we set v''''^ = v''' ^'^ = v'''^ ^(Z — 1,.) in order to 
determine the first order coefficients in the first step of a local iteration. We 
could add external forces in ()542p but we leave it out for simplicity. The 
local convergence of the scheme (with the right choice of p is proved similarly 
as in the first step of the proof of the main theorem above, i.e., by proving 
the convergence of the functional series in the form 

oo 

^p,i = ^pfl,i ^Y^^^pAI^ (543) 

k=l 

where 

(5v'''*^'' i as I 0. (544) 

Note that the initial conditions and the boundary conditions for the func- 
tions dvf''^'^ simplify to 

^6v^''''\t,x) + a^{t,x)5vf''''\T,x) = on [0,RT] x S, 

(545) 

wP^^'\l- 1,.) =0. 

This leads to representations of the solution in terms of the fundamental 
solutions rj, of the equations 

<~i P,k,l o2 p.k.l , -, , n P.k.l 

-g^ Z^j=i + P l^j=i Vj - U. i54bj 

Let us start with the representation for v''''^''. The solution is given in the 
form 

^r°''(^>^) = - l,y)r[,(r,y;0,y)dy 

+ LpE-,„=i L (^^n(y - z)) (^^) (/ - l,z)<izr[,(r,x;.,y)d.dy 
+ Is Ms, y)ro(r, x; s, y)dsdy. 

(547) 
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where (pi is the solution of the integral equation 

^(f)i{T,x)=l I Kr{r,x;s,y)(l)i{s,y)dsdy + fi{T,x) (548) 
Ji-i Js 



2 

along with the kernel 



d 

Kr{T, x; s, y) = ^r[)(r, x; s, y) + Qi(r, x)r[,(r, x; s, y), (549) 
and the functions fi which satisfy 

/i(r,x) = /j^Kr(r,x;0,y)uf'~^(/ - l,y)dy - gi{T,x) 
-^/n/z-i/n^r(r,2;;s,2/) 



(550) 



X E • ^=1 L (^^n(y - z)) (^^) (/ - 1, z)dzdy 



Well for the corrections 6v^' ' these expressions become simplified. We 

(5uf'''''(r, x) = 



+ L P E • ™=i In {£-Kn{y - z)) ( ^^'-^ ) (., .)rf.ri(r, x; y)dsdy 



where (j)^ is the solution of the integral equation 



(551) 



^4{t,x)= I I K^{T,x;s,y)4>Hs,y)dsdy + f,Hr,x) (552) 



2 

along with the kernel 



i-i Js 



d 

KriT,x;s,y) = —T[{T,x;s,y) + ai{T,x)T[{T,x- s,y), (553) 

and the functions which satisfy 
f^Hr,x) = 

-pL!i-iLKr{T,x-s,y) ^^^^^ 
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Well, the functions and (f>'^ have an explicit Levy-type expansion. For 
k = we have it in the form 

i(/)i(r,x) = 

(555) 

fi{r, x) + Em=i //-I Is ^ri'^^ y)fi{s, y)dsdy, 

where 

^f(T, s, ?/) = Krl-r, x; s, y), (556) 

and 

K^+\T,x-s,y)= r [ K^{T,x;a,z)K^{a,z-s,y)dyds. (557) 
Similarly, for A: > we have it in the form 

^(p^{T,x) = 

(558) 

fiiT^ ^) + Em=i Ii-i Is x; s, y)ff{s, y)dsdy, 

where 

Kl'\T,x;s,y) = K^{T,x-s,y), (559) 

and 

K^+^'''{T,x;s,y) = / / Kl'''{T,x;a,z)KP'''{a,z;s,y)dyds. (560) 
Ji~i Jn 

It is possible to extend this work to Navier-Stokes on manifolds. Another 
interesting problem is the extension to equations where the coefficients sat- 
isfy the Hormander conditions and may have a stochastic force term. These 
problems will be studied in a subsequent paper. 
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